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y—i The infrared behavior of propagators and vertices is derived for the maximally Abelian 

gauge and the Gribov-Zwanziger action relying on functional equations. The derivation 

fH and analysis of Dyson-Schwinger equations increase considerably in complexity when going 

I beyond the standard Landau gauge fixing and the available tools have to be improved. 

D For the derivation of the equations a computer program (DoDSE) was developed in order 

'— ' to handle the plethora of terms. The process of determining possible infrared solutions 

^—1 is abstracted to obtain further insight into the structure of the so-called infrared scaling 

> 

\jD solutions. It is found that a few simple steps suffice to determine possible infrared scaling 

f — 
j^ relations directly from the action. This makes also complicated actions easily accessible. 

For the maximally Abelian gauge an infrared enhanced diagonal gluon propagator is 

(^ found, while the off-diagonal degrees of freedom are infrared suppressed. This is in agree- 

^ — ^ ment with the idea of Abelian infrared dominance. Furthermore, it is proven that SU{2) 

^ and higher SU{N) have the same infrared behavior, although the corresponding actions 

S^ differ. Under a suitable truncation the Dyson-Schwinger equations are solved in the deep 

Q^ infrared to obtain values for the exponents of the power laws. The maximally Abelian 

gauge constitutes the first instance where a consistent scaling solution was found besides 

the Landau gauge. 

Restricting the integration in field configuration space to the Gribov region of the Lan- 
dau gauge with the Gribov-Zwanziger action leads to two solutions qualitatively equiv- 
alent to the one obtained with the standard Faddeev-Popov gauge fixing. In both cases 
the gluon propagator is infrared suppressed and the infrared enhanced ghost propagator 



dominates the Dyson-Schwinger equations. This result corroborates the conjecture by 
Zwanziger that the functional integration can be cut at the first Gribov horizon and only 
the applied boundary conditions are important. For one solution the Dyson-Schwinger 
equations reduce in the deep infrared exactly to those obtained with the usual Faddeev- 
Popov gauge fixing. 

For further investigations of Yang-Mills theory the possibility to derive Dyson-Schwinger 
equations with the computer will prove useful. The presented method for determining 
possible infrared scaling solutions will also alleviate their analyses. 
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Introduction 



Symmetries play an important role in physics especially in the theories of the elementary 
forces observed in nature: gravitation, the electromagnetic force, the weak force and the 
strong force. They are described by gauge theories, which are invariant under certain local 
symmetry transformations. The force of interest for this thesis is the strong force whose 
theory is quantum chromodjTiamics (QCD) jlHT]. The elementary particles of QCD are 
gluons and quarks, which build up hadrons like protons and neutrons. To be precise this 
thesis is about the gluonic part of QCD, which is called Yang-Mills theory |8]. 

The overall structure of the Lagrangian of QCD is similar to that of quantum electro- 
dynamics (QED) which describes electrons and photons: The electrons/ quarks have mass 
and interact via the massless photons/gluons. But while electrons and photons can be 
observed as free particles, as attested by the human eye, quarks and gluons only appear 
in bound states like protons or neutrons. This phenomenon of QCD is called confinement 
ini [ID] and will be of special interest in this thesis. 

For the investigation of QCD many different tools are available, each with their own 
disadvantages and advantages. The approach chosen for the present work are functional 
equations. They form an infinite tower of coupled equations, which naturally has to be 
truncated for calculations. However, in some asymptotic cases it is possible to obtain 
general results for the whole system. This provides then constraints for the solutions 
of the equations and can determine the usefulness of a given truncation. Hence such 
information constitutes a reliable starting point for numerical calculations. The main 
topic of this thesis are systems where such solutions can be obtained. 

Functional methods, mainly Dyson- Schwinger equations [TT] - [l3] and functional renor- 
malization group equations |14| . have by now a long tradition in the investigation of 
Yang-Mills theory and also QCD. As it is the most convenient choice in this context, the 
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Landau gauge was usually employed [T5HT8] with a few exceptions. During the last decade 
the methods have improved considerably and while at the beginning several assumptions 
had been required, many of them could be put aside and the results got more and more 
rigorous. Unfortunately the situation is not as clear as one would wish, because two dif- 
ferent solutions are known, which only differ in the deep infrared |19| . The connection 
between these two solutions has been elucidated only partially and it is not settled if both 
or only one solutions are correct. It has even been conjectured that in the former case they 
correspond to two different non-perturbative Landau gauges that coincide perturbatively 

m- 

In Chapter[2]l shortly review these two solutions for the Landau gauge. In the remainder 
of the thesis I focus then on one of them, the so-called scaling solution. Chapter |2] contains 
also an overview over confinement scenarios and aspects of the infrared regime of Yang- 
Mills theory as well as of the gauge fixing procedure. 

In the Landau gauge the methods for the investigation of its scaling solution have 
proven to be very reliable and we have by now a fairly good understanding of many 
aspects of QCD like dynamical chiral symmetry breaking |2n |22] and confinement in this 
gauge. Transferring this knowledge to other gauges seems like a natural next step to 
understand the complete picture of QCD. So it is promising to focus on a gauge with 
different features than Landau gauge. The maximally Abelian gauge j23] provides a 
perfect first candidate for such an endeavour: It is a covariant and renormalizable gauge 
that is closely related to the dual superconductor picture of confinement [2^ [25] . The 
meaning of this confinement scenario in Landau gauge is somewhat unclear, as there are 
no chromomagnetic monopoles which seem to play a decisive role. The findings of this 
thesis about the maximally Abelian gauge, presented in Chapter [5| further elucidate the 
connection to Landau gauge confinement scenarios. 

The second issue investigated is the influence of the Gribov horizon on the Landau gauge 
scaling solution. It arises when attempting a complete gauge fixing, which is required for 
non-perturbative calculations, as the standard Landau gauge definition is insufficient. As 
a consequence the integration in the path integral should be restricted to the interior of the 
Gribov horizon |26| . This is implemented in the Gribov-Zwanziger action by additional 
terms [26| [27] . Interestingly the gluons are then confined already at tree-level as their 
propagator vanishes at zero momentum. The qualitative picture of an infrared vanishing 
gluon propagator and an infrared enhanced ghost propagator is also found by Dyson- 
Schwinger and functional renormalization group equations |T5l - [l8| |28] . However, these 
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functional equations did up to now not explicitly take into account the restriction to the 
interior of the Gribov horizon as it was argued that one can formally cut off the integration 
at the Gribov horizon without changing the equations [29]. This conjecture is explicitly 
tested in Chapter [6] and confirmed by analyzing the Dyson-Schwinger equations of the 
Gribov-Zwanziger action. 

As both the maximally Abelian gauge and the Landau gauge with the Gribov-Zwan- 
ziger action are considerably more complicated than the standard Landau gauge, the 
available methods had to be improved and extended for their investigations. Chapter 
|4] describes this improved method as general as possible so that it can be employed for 
a broad range of different actions and is not restricted to the cases investigated in this 
thesis. To demonstrate its usefulness I give the examples of the Landau gauge, linear 
covariant gauges and ghost anti-ghost symmetric gauges. Under certain assumptions this 
method provides the means to derive all possible scaling solutions for a theory from the 
structure of its interactions. 

Additionally, the basic derivation of the Dyson-Schwinger equations proved to be quite 
time consuming so that I developed a computer program dedicated to this task. The un- 
derlying algorithm, which is suited for manual derivations as well, is explained in Chapter 
[3} where also an introduction to Dyson-Schwinger and functional renormalization group 
equations can be found. The program is named DoDSE ("Derivation of Dyson-Schwinger 
equations") and freely available on the Computer Physics Communications Program Li- 
hrar'^ Given the fields and interactions of a theory it derives the Dyson-Schwinger 
equations in symbolic notation and even the algebraic expressions when Feynman rules 
are provided. 

Finally, I present in several appendices details of some calculations. The derivations of 
formulae used in Section |4] can be found in Appendix \K\ together with the generalization to 
d dimensions. The employed conventions for Grassmann fields are described in Appendix 
[B| Appendix [C] contains the Feynman rules of the maximally Abelian gauge and Appendix 
[D] a formula for the calculation of the sunset diagram. An overview over the use of 
DoDSE is given in Appendix |E] including the calculation of the two-point function Dyson- 
Schwinger equations in the maximally Abelian gauge and for the Gribov-Zwanziger action. 



^ http://cpc.cs.qub.ac.uk/summaries/AECT_vl_0.html 



Chapter 2 

Yang-Mills theory and its infrared 
behavior 



Yang-Mills theory j2] describes the gluonic part of the strong interaction. It is a gauge 
theory and consequently the choice of a gauge is required in functional approaches. This 
amounts to choosing one representative among all the physically equivalent gauge field 



configurations. In Section |2.1| I will describe the standard procedure for gauge fixing, 
as employed in perturbation theory, as well as the problems arising when going to the 
non-perturbative regime. An improved gauge fixing taking these into account will be 
presented in Chapter |6j 

As in this thesis the aspects of the infrared (IR) behavior of correlation functions is 



investigated, I will describe in Section 2.2 two well known ways of connecting the IR 



behavior of propagators to confinement: First I explain the basic idea behind the Kugo- 
Ojima scenario |5U| EI] and then I introduce the Gribov-Zwanziger scenario fIE\, I52H51] . 
The idea behind the latter will be important for Chapter [6j 

In the final section of this chapter I will shortly review the development of our un- 
derstanding of the IR regime in the Landau gauge. This gauge is widely preferred in 
functional approaches because of its simplicity. Indeed calculations in all other gauges 
lead to one or more additional obstacles that have to be overcome. But although the 
Landau has been intensively investigated for many years by now, there is no consensus 
yet about the solution for the Green functions in the IR. Two possible candidates have 
emerged, and even the variant that they correspond only to different non-perturbative 
extensions of the Landau gauge is discussed. I will describe the features of both solutions 
in this section. Furthermore, I will also explain the notion of positivity violation, which 
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is a mathematical criterion for confinement. 

2.1 The action of Yang- Mi I Is theory 

In this section the Lagrangian of Yang-Mills theory is introduced as the gluonic part of 
QCD. How to fix the gauge, the BRST symmetry [5314H7] and properties of the Gribov 
region are described. 

2.1.1 The gauge invariant action 

The Lagrangian density of QCD is jH |H] 

Cqcd = qi-Jp + m)q + Cym, (2.1) 

Cym = ^tr{F^,F^,}. (2.2) 

The first term in Cqcd describes the propagation of the massive quarks and their inter- 
action with gluons via the covariant derivative D^: 

D^ = d^ + igA^. (2.3) 

The second term, Cym, constitutes the purely gluonic part of QCD and thus contains the 
propagation of gluons and their self-interactions. While it is not possible to find a gauge 
invariant Lagrangian for quarks without gluons, this is not true in the opposite direction, 
i.e., Cym is gauge invariant and can be considered on its own without quarks. It is this 
term that is investigated in this thesis. 

The quantity F^j, is called the field strength tensor and is given by 

Ff,u = d^A, - d,A^ + I g [A^, A,] . (2.4) 

The gauge field A^ lives in an algebra defined by the hermitian generators T"^ of a generic 
gauge group. They obey the relations 

[T'^T'] = ir'T\ (2.5) 

tr{T'T'} = TfS'', (2.6) 
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with Tf = ^ for SU{N). The decomposition of the gauge field is 

A, = A^T (2.7) 

and similarly for the field strength tensor: 

F,, = f;x\ (2.8) 

F;, = d,A: - d^A; - g r'A^^Al (2.9) 



In components the YM Lagrangian reads 



1 

4" 



^YM - -FJ^^FJ^^. (2.10) 



The transformation of the gluon field under which the QCD Lagrangian is invariant is 
given by 

A^^ix) = U{x) A^{x) Uix)'' + -{d.Uix)) U{xr\ (2.11) 

where U{x) is 

U{x) = e*3"^(^) (2.12) 

with uj{x) the Lie algebra valued gauge parameter: 

co{x) = uj'{x)T\ (2.13) 



Eq. (2.11) is called a gauge transformation. In infinitesimal form it reads 

a;; ^ a; + m; = a; - d.oo^ - g r'u-Al = a;- d;^u^, (2.14) 

where the covariant derivative in the adjoint representation D[^^ is defined as 

D;' = 5^^d, + gr'A% (2.15) 

As mentioned above the gauge group for the strong interaction is the compact simple 
Lie group SU{3). In this thesis, however, I will for a large part adopt the gauge group 
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SU{N). The reason is that not only SU{3) per se is of interest, but also general SU{N) 
and especially SU{2). In the latter case calculations can become easier and especially 
in the maximally Abelian gauge SU{2) is often employed as the Lagrangian simplifies 
considerably, see Chapterls} General SU{N) or also other groups like G2, see, for example, 
refs. [38| 139] . are investigated to determine the dependence of certain aspects of Yang- 
Mills theories on the gauge group. Also investigations in the so-called large A^ limit are 
performed 



2.1.2 Fixing the gauge 



To work with functional methods we cannot use the action eq. (2.10) as problems appear 
at several points. For example, in the canonical quantization, see, e.g., ref. [1], the 
commutator relations cannot be fulfilled, as the conjugate momentum of A^, given by 

has no zero component IIo. But this contradicts the equal-time commutation relation 

[A^{x),U,{y)]=tg^,5ix-y). (2.17) 

If, however, a certain gauge is chosen, Hq does no longer vanish and one can proceed as 
usual. 

In the path integral the source of the problems is that one integrates over all gauge 
equivalent gauge field configurations, although the idea is only to take into account gauge 
non-equivalent, i.e., physically different configurations. In other words the defining sym- 
metry of the theory makes the integration in the path integral overcomplete. Furthermore, 
one cannot calculate the propagator of the theory, as the gluon two-point function has 
zero eigenvalues. Consequently it is singular and cannot be inverted to get the propagator. 

It is useful to introduce here the notion of a gauge orbit, which is the set of field 
configurations connected by a gauge transformation: 

0[A]:={A'^\A'^ = A^^}. (2.18) 



A^ is given in eq. (2.11). Ideally one takes only one representative per gauge orbit. The 



idea by Faddeev and Popov to achieve this was to restrict the integration in the path 
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integral to a hyperplane [H]. This can be done by inserting unity, given by 

1 = /\[A]JVU5U[A'']), (2.19) 

into the path integraL The delta functional defines the hyperplane f[A\ =0 and T)!! is an 
integration over group space. The factor A [A] is the Jacobian arising from the transition 
from field variables A to gauge transformation variables U . The integration over U can be 
absorbed into the normalization of the path integral and the factor A [A] can be calculated 

as 

A"[A] = det (^£1^) =■ detM^^{x,y), (2.20) 

where the color indices have been made explicit. M{x, y) is known as the Faddeev-Popov 
operator. Its determinant can be localized by the introduction of a pair of Lie algebra 
valued Grassmann fields, the so-called Faddeev-Popov ghosts c and c: 

.e*M"(..„^/^Ne/-.— <--.). (2.21, 

Although they have zero spin, they obey Fermi statistics, which underlines their status 
as auxiliary and not physical fields. 

As an alternative to the restriction to a hyperplane one can relax this condition to a 
Gaussian distribution over the gauge orbit with the mean value defined as the original 
gauge fixing condition. One example are linear covariant gauges with a gauge fixing 
parameter ^ that gives the width of the distribution around the Landau gauge fixing 
condition d^^A^ = 0. In the limit of ^ — )■ the Landau gauge is recovered. One should 
note that such a Gaussian distribution includes all gauge copies. Since it is normalized it 
does not lead to the same problems as the unnormalized integration before. 

The gauge fixing part which is added to the Lagrangian density is then given by 



^9f = ^(5.^;(^))' - J dyc^ix)M^%x,y)c%y). (2.22) 

Also in the Landau gauge the full expression has to be added, 

CyM -^ ^YM + ^gf, (2.23) 
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and only after inverting the gluon two-point function one can set ^ = 0. In the path 
integral the factor exp{— ^{d^A"^^)^) reduces to a delta functional in the limit ^ — )■ 0. 

2.1.3 The Faddeev- Popov operator and the Gribov region 

In the remainder of this section we will confine ourselves mostly to the Landau gauge in 
order to discuss some specific properties of the Faddeev-Popov operator and the Gribov 
region. In Landau gauge it is given by 

M'^ix, y) = M''{x)5{x - y) = -d,D;'6{x - y). (2.24) 

It is easy to demonstrate that the Landau gauge fixing is not complete as there are still 
gauge copies left. However, the fact that the gauge is not fixed completely by a local 
gauge fixing condition is not a specific property of the Landau gauge, but valid for all 
local gauge fixing conditions |12]. Starting with a gauge configuration fulfilling the Landau 



gauge condition, d^j,A^ = 0, we can perform a gauge transformation according to eq. (2.11 ) 
demanding that the result again fulfills the Landau gauge condition: 



d^A^ = 0^d^A^-d^D^uj = => Muj = 0. (2.25) 

Hence if the Faddeev-Popov operator has zero modes, there are still gauge equivalent 
configurations left. 

The first to treat this issue was Gribov [SB]. He proposed to restrict the configuration 
in the path integral to the region Q, where the Faddeev-Popov operator M{A) is strictly 
positive: 

n := {A\d^A^ = 0, M{A) > 0}. (2.26) 

Positivity entails that the operator is invertible, which is important as its inverse corre- 
sponds to the propagator of the Faddeev-Popov ghosts. Today the region Q is known as 
the Gribov region, which has the following properties in Landau gauge: 

• The vacuum configuration, A^ = 0, lies in the Gribov region [43j. This can easily 
be verified by observing that it fulfills 9^A^ = and the Faddeev-Popov operator 
reduces to the Laplacian, M(0) = — D, which is a positive operator. Since pertur- 
bation theory is an expansion around A = 0, one can understand that it yields good 



13 



2. 1 The action of Yang-Mills theory 



results, as long as the quantum fluctuations do not become large enough to feel the 
presence of the Gribov horizon. 



• The Gribov region is bounded in all directions jl3]. In this it differs decisively from 
the Gribov region of the maximally Abelian gauge |^ H5] . which will be discussed 
in Section [5.3.31 

• It is convex, i.e., two arbitrary conflgurations Ai and A2 within the Gribov region 
can be combined to a new conflguration ^3 that lies again within the Gribov region 
as follows [i3] : 

A3 = aAi + (1 - a)A2, < a < 1. (2.27) 

• Every gauge orbit passes at least once through the Gribov region |l6] . This property 
is important so that one can restrict the integration to the Gribov region without 
missing any gauge orbits. 

• Unfortunately there are still gauge copies within the Gribov region so that the gauge 
flxing is not complete ^7\ HH]. However, expectation values are not influenced by 
these additional copies [43]. 

The statement that every gauge orbit passes through the Gribov region is related to 
the fact that the Landau gauge condition can be derived by minimizing the functional 

R[A]:=^fdxA^{x)A^{x) (2.28) 



with respect to infinitesimal gauge transformations, eq. (2.14): 

5R[A]= f dxA'^{x)5A'^{x) = j dxAl{x){-Dl'{x)uj'{x)) = 
dxA;{x){-d,u^{x)-gr'Al{x)u^{x)) = 
dx{uj^{x)df,A'^ix)). (2.29) 

In order to have an extremum (9^A^ must vanish. The type of the extremum can be 
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determined by the second derivative: 

5^R[A]= I dx{uj'{x){-df,Dl'{x)uj'^{x))= j dxuj'{x)M'\x)uj'{x). (2.30) 

Here we recognize the Faddeev- Popov operator M^*, i.e., if it is positive, we have a 
minimum of the functional R[A\. Thus the restriction to the Gribov region corresponds, 
as asserted above, to the minimization of the functional R[A\. 

As the Gribov region is still plagued by gauge copies ^1\ HH], one can think about 
better ways of fixing the gauge. A natural choice for a unique gauge fixing is the set 
of gauge field configurations that corresponds to the absolute minimum of R[A\. This 
region is known as the fundamental modular region. It possesses a topologically non- 
trivial boundary where degenerate global minima exist that have to be identified. Taking 
the global minimum is also known as absolute Landau gauge and taking one arbitrary 
minimum as minimal Landau gauge. In lattice calculations one can employ algorithms 
that correspond to the latter or approximate the former and indeed finds an influence on 
the correlation functions |1^ 150] . More on different Landau gauges on the lattice can be 
found in Section [6. L3I 

2.1.4 Replacing gauge symmetry by the BRST symmetry 

In the action fixed to a hyperplane in field configuration space gauge invariance is ex- 
plicitly broken. However, there is another symmetry that takes its place. It is named 
BRST symmetry after Becchi, Rouet, Stora [36| [37] and Tyutin [35j. It is very useful in 
proving renormalizability and unitarity of a theory, see, for example, refs. [HTl IST) . The 
corresponding transformations can be derived from the standard gauge transformation by 
replacing the gauge parameter a; by a ghost field c. Due to this choice the pure Yang-Mills 
part is trivially invariant. The invariance of the gauge fixing part is explained below. 

It is convenient to introduce the gauge fixing condition via an auxiliary Lie-algebra 
valued field h that takes the role of a Lagrangian multiplier: 



j^lf-^-JdxwnA]) ^ N5{r[A]). (2.31) 

It is not dynamical and called Nakanishi-Lautrup field [52|l53]. Relaxing the gauge fixing 
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condition f[A] = into f[A] =i^b/2 leads to 

'j^^r^-Jd.{ib^rlA] + lbrb^)^^^~i-^if[A])^ (2.32) 



/ 



^^; 


= 


-D-c% 




sc^ 


= 


-\gr'c' 


'c* 


S(f 


= 


ih\ 




slf 


^ 


0. 





and corresponds to the Gaussian averaging over the gauge orbit in hnear covariant gauges 
as described above. A^ is some normahzation factor and ^ is called a gauge fixing pa- 
rameter. It determines the width of the Gaussian distribution. Specifically .^ = is the 
original gauge fixing condition. However, one should keep in mind that this is only well 
defined if the solution to f[A\ = is unique, what is not the case non-perturbatively. 
In Landau gauge the off-shell BRST transformation reads: 

(2.33) 

(2.34) 

(2.35) 
(2.36) 

Integrating out the Nakanishi-Lautrup field yields the on-shell form, where we just have to 
replace h^ by —i (9^Ap/^. A substantial property of the BRST symmetry is its nilpotency, 
s^ = 0, in the off-shell case. Thus BRST can also be defined as follows: The gauge 
parameter of the gauge transformation is taken as the anti-commuting field c. Requiring 
nilpotency of this transformation, i.e., s'^A^ = 0, determines sc. Finally the fields c and 
b are introduced as a BRST doublet, i.e., they have the trivial BRST transformations 
sc = ib and sb = 0. 

The nilpotency property allows an easy way to fix the gauge without the need for 
a path integral |54j- This is based on the observation that one can add any quantity 
that is the result of a BRST transformation, a so-called BRST exact quantity, to the 
Lagrangian without spoiling its BRST invariance. As the gauge fixing condition f[A] has 
ghost number zero and the BRST transformation itself raises the ghost number by one, 
we introduce the factor c to get ghost number zero and add s(c/[y4]). In the case of the 
Landau gauge, where f[A] = 9^A^, this prescription to fix the gauge directly leads to the 
known gauge fixing terms: 

C,f = s{e-f[AY) = sic-dM = 
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= I h' {d^AD - c'- d^i-Dl'c') = I ¥{d^Al) - <f ikf^ c\ (2.37) 

One minus sign stems from the anti-commutativity property of s and c. This method can 
be used also for other gauges and allows the use of gauge fixing conditions depending on 
ghost fields [54J. It will be employed for the maximally Abelian gauge in Chapter pi where 
it significantly simplifies the gauge fixing procedure. As the expectation value of any gauge 
invariant quantity remains unaffected by adding such a BRST exact form, one can very 
nicely see in this way of quantization that all physical observables are independent of the 
chosen gauge. 

2.2 Aspects of the asymptotic infrared regime 

Solving a Dyson-Schwinger equation (DSE) numerically for the complete momentum re- 
gion is a challenging task and it proves useful to know the asymptotic behavior of the 
correlation functions [T6| |55] . For the ultraviolet (UV) perturbation theory naturally 
provides a good guideline. For the IR regime it is more complicated to determine the 
qualitative behavior and thus it is worthwhile to develop analytic methods for its inves- 
tigation. But the low momentum region is not only of interest for providing input for 
numerical calculations. One can also learn about the realization of different scenarios of 
confinement [9]. In this section I will give a short overview over confinement scenarios 
directly related to correlation functions. The description of an additional scenario, the 



dual superconductor picture, is deferred to Section 5.1 as it does not translate directly into 



conditions for the correlation functions. A purely mathematical criterion for confinement. 



violation of positivity, is described later in Section 2.3 



2.2.1 The Gribov-Zwanziger confinement scenario 



As mentioned in Section |2.1.2[ the usual gauge fixing is not sufficient for the non-per- 
turbative regime. An improvement can be achieved by restricting the integration in field 
configuration space to the first Gribov region. This is explicitly realized by the so-called 
Gribov-Zwanziger action |261I27] . Its derivation and IR analysis are described in Chapter 
[6j Here I only mention the qualitative consequences of this improved gauge fixing for the 
theory. 

The main statements of the Gribov-Zwanziger confinement scenario in the Landau 
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gauge are that the gluon propagator vanishes at zero nionientum |34| and that the ghost 
propagator is IR enhanced, i.e., it diverges stronger than a simple pole [32| l33] . The 
interesting point is that these statements are realized in the Gribov-Zwanziger action 
already at the perturbative level. The reason is that the existence of the Gribov horizon, 
which is a non-perturbative object, is taken into account. The tree-level gluon propagator 
is 

where A^ is the number of colors and 7 a mass parameter. 7 is not free but determined 
by the horizon condition, see Section 6.1.2[ This condition has to be enforced in order to 



make the theory well-defined in the first place |26) . 

Naively the bare ghost propagator goes like l/p"^, but it turns out that the horizon 
condition leads to a cancelation of this term at one-loop level and the IR leading part of 
the ghost propagator goes like 1/p^. One can show this cancelation diagrammatically for 
the exact ghost propagator |56| and perturbatively it was checked up to two loops in ref. 

E 



2.2.2 The Kugo-Ojima confinement scenario 

Although being based on completely different considerations as the Gribov-Zwanziger 
scenario, the mechanism proposed by Kugo and Ojima leads qualitatively to the same 
predictions for the IR behavior of the gluon and ghost propagators [HUj [HT] . The basis 
of the Kugo-Ojima construction is the existence of a global, non-perturbatively defined 



BRST symmetry. It was long not clear if the standard BRST, as given in eq. (2.33), could 
be extended in this sense to the non-perturbative regime. Recent progress in this respect, 
however, leads exactly to this conclusion [581461] . 

The importance of the BRST symmetry rests on the fact that the corresponding charge 
defines the physical state space of the underlying quantum field theory by its cohomology. 
In contrast to the total state space, which possesses an indefinite metric, the physical 
state space must have a positive metric in order to allow a probabilistic interpretation 
of expectation values. According to Kugo and Ojima physical states G Vphys obey the 
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condition 

QB\(f)) = 0, (2.39) 



where Qb is the Noether charge corresponding to the BRST symmetry |jT| |62] . Further- 
more, one can distinguish between states with positive and zero norm, \1/ and x ^ ^o, 
respectively: 

(^1^) > 0, (xlx) = 0. (2.40) 

The zero norm states x are also called daughter states, as they are the BRST variation 
of unphysical states, the parent states \a): 

\x) = QB\a), \a)^VpHys. (2.41) 

Consequently the subspace Vq € Vphys is orthogonal to Vphys and does not influence any 
amplitude in Vphys- The physical Hilbert space equipped with a positive metric is then 
given by the completed quotient space 



Hphys — Vphys/Vo. (2.42) 

In summary states fall in one of the following three classes: 

• BRST non-invariant states \a): These are unphysical states, as the BRST charge 
does not annihilate them, Qb\c() y^ 0. 

• BRST exact states |x): They are the BRST variation of the states |a). 

• BRST invariant but not BRST exact states |\I^): They define real physical states, 
i.e., they are BRST invariant, QbI"^) = 0, and have a positive norm, (^I/|^E') > 0. 

Kugo and Ojima could show that the states \a) and \x) always appear in so-called 
quartets, which is given by a "Faddeev- Popov conjugated" pair of BRST doublets |31] . 
Using the notation \k,N), where k represents all quantum numbers except for the ghost 
number N, the four particles of a quartet can be denoted as follows: 

|0i) = |A:, iV), \^,) = Q^\cf>,) = \k,N-l), 
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103) = \k, -N), 104) = gB|03) = \k, -{N + 1)). (2.43) 

|0i) and 102) are one BRST doublet, |03) and |04) the other. Furthermore, the states 
|0i) and 103) belong to the class of BRST non-invariant states, and |02) and |04) to the 
BRST exact states. The contributions of these four states always cancel and they will 
never appear as asymptotic observable states. 

The only other possibility of states are states ^ that are BRST invariant but not BRST 
exact, i.e., they obey QbI"^) = and cannot be written as a BRST variation of another 
state. These states are called singlet states and should represent the physical states of the 
theory. For QCD this means that the fundamental, confined fields of quarks, gluons and 
ghosts should belong to some quartet, whereas physical states such as mesons, baryons 
and glueballs should be BRST singlets. 

Based on the global color charge Kugo and Ojima derived a simple criterion for color 
confinement in Yang-Mills theory. Using the equation of motion of the gluon field the 
conserved global color current can be written as 

J: := 9,F;, + {Qb, Dfc^}, d,J: = 0. (2.44) 

The related charge is the global color charge 

Q- := G" + N'':= f d^'x [d^F^, + {Qb, D^c'}) , (2.45) 

where G"' and N"- correspond to the first and second terms in the integral. The first 
criterion of Kugo and Ojima is that diF^j^ contains no discrete massless pole. The charge 
G" then vanishes, as the integral is over a total derivative. If there were massless con- 
tributions, G"' would be ill-defined. The second criterion requires that N"- is zero and 
consequently the total color charge Q"- is also zero. 

In the charge A^" the anti-ghost field appears. It belongs to the so-called elementary 
quartet and has a massless contribution from the asymptotic field 7". The asymptotic 
contributions due to the composite operator gf°'^^A'i(f are given by u'^'^d^'^'^ characterized 
by the dynamical parameter m"^, which can be determined from the correlation function 

dxe'^^--y\{D;^c^) {x)gf^'At{y)c\y)) = (^g,^ - ^^ u'^\p'). (2.46) 
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The condition for the charge A^" to be well-defined can be inferred from the asymptotic 
behavior of {D'^''c^){x), given by 

(D»'c')(x) "" =^°° (5^^' + u^'')d^f{x), (2.47) 

as 

^ab ^ ^afe(Q) _ _^ab_ ^2.48) 

In summary the two confinement criteria due to Kugo and Ojima are: 

1. There is no discrete massless pole in diF^^. 

2. The parameter m"^ has to be —6""^. 

The last point was taken up again by Kugo sixteen years after his article with Ojima, 
ref. ini]. In ref. [SU] he derived how the Landau gauge ghost propagator, parametrized 

by 

Df (p2) = _s-b?^ (2.49) 

is related to the parameter u: 

Cc(0) = -^. (2.50) 

1 + u 

Hence the ghost propagator is more IR singular than a simple pole, if the criterion u = —1 
by Kugo and Ojima is met. 

2.3 Solutions of Landau gauge Yang-Mills theory in 
the infrared 

The non-perturbative calculation of propagators and vertices is involved and in the course 
of history the emerging picture underwent some changes and improvements. The best 
investigated gauge in this context is Landau gauge for which I will give a short review in 
this section. 
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2.3.1 From infrared slavery to ghost dominance in the Landau 
gauge 

The Landau gauge is a preferred gauge in functional approaches for several reasons. First 
of all, it is the gauge with the simplest form of the action. There are only two fields which 
interact via three vertices. This simplicity is supported by additional information such 
as, for example, the fact that the system of the transverse components of Green functions 
is closed, see, e.g., ref. [19j, and the longitudinal terms can be discarded. Furthermore, 
there have been arguments that the ghost-gluon vertex stays bare in the IR [UH j IM ] . This 
was used as input at the beginning, but having understood the gauge better and better 
it could be proven directly from functional equations [651 - I67] . 

A well-known example of early non-perturbative calculations of the gluon propagator 
in the Landau gauge is due to Mandelstam |68]. Motivated by perturbation theory he ne- 
glected the ghost loop, which only contributes a small amount at high momenta, and the 
four-gluon vertex, which only occurs at two-loop level. With an approximated three-gluon 
vertex he obtained an IR divergent gluon propagator and also an IR singular running cou- 
pling. As such a gluon propagator apparently is perfectly suited to yield a linear rising 
potential by single gluon exchange, this picture was widely accepted and became known 
under the name IR slavery. The results of the Mandelstam approximation were subse- 
quently confirmed in refs. [691 - 171] . However, this is a perfect example of how assumptions 
based on perturbation theory can lead in the wrong direction in the non-perturbative 
regime. In the late nineties it was shown in refs. [T^ FIB] that the ghost propagator yields 
the most important contribution to the IR sector of Landau gauge Yang-Mills theory. In 
fact, it was found that the ghost contributions were dominating all DSEs |66 t 1721 173] sup- 
porting the scenarios of ghost dominance as proposed by Gribov and Zwanziger f2E\ IB5] 
and Kugo and Ojima [301 131] . 

This solution is characterized by power laws for the dressing functions ca(p^) and Ccip"^) 
of the gluon and ghost propagators, respectively: 

DA{p')=dA-{py^ D,{p') = d,-{py^. (2.51) 

dA and dc are momentum independent coefficients and the qualitative behavior is deter- 
mined only by the exponents, called infrared exponents (IREs). An important feature of 
these power laws is that the exponents depend only on one parameter, usually denoted 
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by K, as they are related by a so-called scaling relation: 

K ■= -5c = 5a/2. (2.52) 



The value of k can be calculated analytically. The most reliable value is 0.5953 . . 

It is obtained for a bare ghost-gluon vertex and changes only slightly when an IR constant, 

but momentum dependent dressing is employed |65| . 

The fact that the ghost-gluon vertex is not IR enhanced was originally used as an input 
to solve the DSEs. It relied on an argument by Taylor |63| IM] and was confirmed by 
lattice [7H477] and DSEs studies |78l 179] . The self-consistency of this assumption was 
confirmed later on [72] and finally it was shown that this is indeed the only possible 
IR solution |66| |67] . This proof was possible by combining the two distinct systems of 
functional renormalization group equations and DSEs as suggested in |80j. A different 
approach with the same result was taken in ref. |73| . where as an assumption the existence 
of a stable skeleton expansion in the IR was used. It turned out that such an assumption 
necessarily holds as can be derived from functional renormalization group equations [67J. 



Details on the proof can be found in Section 4.2.3 where this emerges merely as a side 



result from the general analysis of the systems of functional equations. 

2.3.2 The decoupling solutions and the scaling solution 

Although the uniqueness of this solution in terms of scaling laws was established, the 
results did not agree with calculations on the lattice, see, for example, refs. [7S| IHU- 
[QT] . For exceptions see refs. [89l [92] . the two-dimensional case [76| l87] and the strong 
coupling limit |93ti95] . Indeed another solution was found with continuum methods that 
possesses quite different characteristics: The ghost propagator is not IR enhanced but 
stays bare in the IR, and the gluon propagator becomes finite instead of going to zero 
[T9l [96VI99) . This behavior is responsible for the name decoupling solution, as the gluons 
decouple below the scale given by their mass. One should note that a massive behavior 
for the gluon propagator is also possible in the formerly found solution, now called scaling 
solution, since the value of the parameter k could be 1/2. But then the ghost propagator 



would still be IR enhanced via the scaling relation eq. (2.52). In the decoupling solution, 
however, the propagators do not obey an IR scaling relation and the ghost stays bare 
in the IR. Nevertheless there is some realization of ghost dominance in the decoupling 
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Figure 2.1: Numerical solutions of the gluon and ghost two-point DSEs. Plotted are the 
gluon (left) and ghost (right) dressing functions for the decoupling and scaling 
solutions. Printed with permission from Elsevier from ref. |19] . 



solution, as the gluon is IR suppressed. Another qualitative difference between the two 
solutions is that the scaling solution is unique, while there exists a family of decoupling 



solutions |19| . In fig. 2.1 both solutions are plotted. 

The connection between the two solutions is the renormalization of the ghost propaga- 
tor DSE [19]. Solving the DSEs one requires renormalization conditions and for the ghost 
propagator it is convenient to specify the renormalization condition at zero momentum. 
Choosing a finite value for the ghost dressing function yields a decoupling solution, where 
the value of the ghost dressing function at zero momentum is connected to the mass of 
the gluon (TU]. For an infinite ghost dressing function at zero momentum the scaling 
solution is recovered [T5l [T6| [T9] . One should note that it is not necessary to fix the value 
at zero momentum but the normalization prescription can be defined at any arbitrary 
momentum, since the values for different momenta are uniquely connected. Mathemati- 
cally the necessary choice of the renormalization condition for the ghost dressing function 
corresponds to a boundary condition for the DSEs. 

This is also related to an argument by Zwanziger that one can in principle cut the 
integration in field configuration space at the first Gribov horizon without modifying the 
form of the equations [29]. What is changed, however, are the boundary conditions. In 
case one cuts at the first Gribov horizon one uses as boundary condition that the ghost 
propagator is more enhanced than a simple pole, as derived from the horizon condition 
|32| |33| [56] . One can also argue for the same boundary condition with the Kugo-Ojima 
confinement scenario [30l l3T] without referring to the Gribov region. Unfortunately it 
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is not understood how these two arguments are related, since the implementation of the 
Kugo-Ojima condition u = —1 into the action leads to the Gribov-Zwanziger action [lOOj. 
However, the Gribov-Zwanziger action breaks the BRST symmetry, while the Kugo-Ojima 
formalism relies on an unbroken BRST charge. 

What speaks in favor that both solutions are valid is the fact that physical quantities 
do not seem to depend on the deep IR, the only region where the two solutions differ from 
each other. In ref. |1U1] the transition temperatures of the confinement /deconfinement 
and the chiral transitions at vanishing chemical potential were calculated with identical 
results for both solutions. 

In the present work I concentrate only on the scaling solution, which is more accessible 
to analytic approaches and offers explanations for many aspects of Yang-Mills theory. The 
most important one is certainly the confinement of gluons which can directly be inferred 
from the IR vanishing gluon propagator, as one can show analytically that the Schwinger 
function A(t) is not positive |34| . This, however, is necessary to interpret a particle as 
physical |102[ 1103] . It boils down to show that 

A(t) := J d'x J (^e^^ V(/) > 0, (2.53) 

where cr(p^) is a scalar function characterizing the propagator. The Fourier transformation 
of the gluon propagator at zero momentum is related to the Schwinger function by 

L)(p2 = 0) = /" d^xD{x -y)= [ dtA{t). (2.54) 

Thus a D{p^ = 0) = 0, A(t) is either zero or it has equally distributed positive and 
negative contributions. In this case one speaks of maximal positivity violation. But also 
for the decoupling solution a violation of positivity is observed both by lattice simulations 
|85| 11041 1105) and functional equations [19j . 
A delicate point is the global BRST symmetry of the solutions. If a global symmetry 



exists that has the same form as the perturbative definition, eq. (2.33), we know from 
the seminal paper of Kugo and Ojinia how the propagators of the ghost and the gluon 
should behave in the deep IR [301 |3T] : The ghost propagator should be IR enhanced and 
the gluon propagator should vanish at zero momentum. Only the scaling solution fulfills 
these criteria. However, one cannot exclude a different realization of a global BRST that 
leads to a different qualitative behavior in agreement with the decoupling solution. 
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Finally I would like to note that both the decoupling and the scaling solutions fulfill 



the criterion for a confining Polyakov loop potential |106) and thus quarks are confined. 
This is a sufficient criterion depending only on the asymptotic IR part, but in actual 
calculations it is found that the region responsible for confinement is the mid-momentum 
regime, where both solutions agree |106] . 
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Chapter 3 
Functional equations 



This chapter serves as a short introduction to Dyson-Schwinger (DSEs) and functional 
renormahzation group equations (FRGEs). Both systems of equations form complete 
sets describing the theory exactly. As such they are suited for the investigation of non- 
perturbative phenomena. 



The derivation of DSEs is described in Section |3.1[ First some general relations among 
the generating functionals are provided in order to allow a self-contained description. Then 
the process of deriving the equations is explained and the computational tool DoDSE 
based on the presented algorithm is introduced. 



FRGEs are treated in Section 3.2 They require the introduction of the so-called effec- 
tive average action which depends on an artificial momentum scale and can be interpreted 
as an action interpolating between the classical and the quantum action. After introducing 
the effective average action, the derivation of FRGEs is sketched. 



3.1 Dyson-Schwinger equations 

Dyson-Schwinger equations (DSEs) are named after F. J. Dyson [llj and J. S. Schwinger 
[T2I IT3] . Put into one sentence, DSEs are the equations of motion of Green functions 
and describe the propagation and interaction of the fields of the theory. The full system 
of DSEs provides a complete description of the theory. Therefore DSEs can be used 
to generate a perturbative expansion in the weak coupling regime, but they show their 
true strength when applied for strong coupling. This makes them a perfect tool for 
the investigation of Yang-Mills theory and also QCD, where non-perturbative methods 
are needed to explain such phenomena as confinement and dynamical chiral symmetry 
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breaking. Both are important properties of these theories and cannot be accounted for 
by standard perturbation theory. 

3.1.1 The generating functionals 

The easiest way to derive DSEs is via the path integral A derivation with the canonical 
formahsm is possible, but less instructive and definitely more tedious. The interested 
reader is referred to ref. [107]. The path integral also allows for an easier algorithmic 
implementation. 

In the following we will consider a theory with fields (pi. The multi-index i contains the 
field type, all indices of the fields, as can be, for example, Lorentz or color indices, and 
also the space-time or momentum dependence. Equal indices in a product are summed 
and integrated over. This avoids cumbersome notation. The action is therQ 

11 1 

'S'M =W-^Srs(pr4>s - ly^Srst4>r4>s4>t " Y\^rstu4>r4>s4>t<t)u- (3.1) 

The statistical factors are chosen such that the coefficients Srs-, Srst and Srstu denote the 
bare two-, three and four-point functions, and the choice of signs is a consequence of 



the definition of the vertices, see eq. (3.6) below. From this action the path integral is 
constructed as 

Z[J] = [ D[(f)]e-^^'^'-^' = e^W, (3.2) 



where Jj is the source of the field 0j. The path integral Z[J] is also called the generating 
functional for full Green functions and W[J] that for connected Green functions. A 
Legendre transform of W[J] with respect to the sources yields the so-called effective 
action that generates the one-particle irreducible (IPI) Green functions, which are those 
Green functions that are still connected after one internal line is cut: 

r[$] = -W[J] + ^,J,. (3.3) 



""^Here we restrict ourselves to three- and four-point functions, but the generahzation is straightforward. 



28 



3.1 Dyson-Schwinger equations 



It depends on the averaged fields $ in the presence of external currents J, 

$, := (0.), = ^^ = Z[J]-' j Dm.e-'^^^'K (3.4) 

The current, on the other hand, can be expressed as the derivative of the effective action: 

J^ = -^- (3.5) 

Derivatives of the effective action with respect to fields are abbreviated as r/^...j^: 

r;'.:^-.:.^^^. (3.6, 



The minus sign is chosen for later convenience such that no additional minus signs in the 
DSEs of vertices appear. The T/^...^^ are not yet the physical n-point functions of the 
theory as the external sources J are still non- vanishing. Setting them to zero physical 
propagators Dij and vertices Ti-^...i^ are obtained: 

A, := Af °, (3.7) 






^h-in '■— ^i-,-i„- (3.^ 



The function D':^, is given by 



6^w[j] /52r[$] 



^j 



-1 



^i ■ 6mj! [ 5$,i$!- ' ■ ^^-^^ 



It is important in the derivation of DSEs that the sources J are set to zero only at the 
end. Otherwise one would miss contributions. 

The propagators Dij are the inverse of the two-point functions Tij. It should be stressed 
that this is a matrix relation, i.e., if the two-point matrix is not diagonal there is a non- 
trivial relationship between propagators and two-point functions. This complicates anal- 
yses of actions with mixed propagators like the Gribov-Zwanziger action. For a diagonal 
matrix in case of bosons or an off-diagonal matrix for fermions the situation is simpler, 
since the propagator can be directly calculated as the inverse of the corresponding two- 
point function. 
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3.1.2 Derivation of Dyson-Schwinger equations 

Using our multi-index field the formal derivation can be done in a few lines. The real 
work is then to expand in the physical fields of the theory. DSEs for full, connected and 
IPI Green functions can be worked out, but since the full theory can be reconstructed 
from IPI correlators we focus on those. 

To start we note that the integral of a total derivative vanishes so that 



/^M 



6S 



+ Ji]e 



-S+<f>jJj 



6S 



+ J, Z[J] = 0. (3.10) 



<P',=5/5J, 



Employing further derivatives with respect to sources yields the DSEs of the full Green 
functions. Substituting Z[J] by e^t*^' and using 



-W[J] 



SJi 



,W[J] 



6W[J] _5_ 
5Ji 5Jj 



(3.11) 



after multiplication of eq. (3.10) from the left with e ^^"^ we find 



6S 



+ J, = 0. 



(3.12) 



SW[J] S 



This is the generating DSE for connected correlation functions. The DSEs of connected 



Green functions are obtained by acting with further source derivatives on eq. (3.12). To 



get the corresponding version for IPI functions we perform the Legendre transformation 
of W[J] with respect to all sources. Thereby 6W[J]/6Ji changes to $i and 6/6 Ji becomes 



6 _ 6^j 6 
6Ji 6Ji 6^j 



6 6W[J] 6 _ 6^W[J] 6 _ j^j 6 



6Ji 6Jj 6^j 



(3.13) 



This yields 



6S 



b^=^^+DJ 5/5<S>j 



6Tm 

H — = 



(3.14) 



which is the basic equation. All DSEs for IPI Green functions can be derived from it by 
further differentiations with respect to the fields. 
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Let us write down this expression in the case of an action with three- and four-point 



interactions as given in eq. (3.1). The derivative is simply 



SS 



Ji: 



2! 



S., 



istVsVt 



3! 



Jir. 



irstVsHJtHJu- 



(3.15) 



Before replacing the field 0j by $, + D/t^ we need to know how the differentiation oper- 



«i5$ 



ator j^ acts on fields and propagators. For completeness we also include the derivative 



of an n-point function: 



$. = A.. 



-£>• 



(3.16a) 



52r[$] 






5^j5^r 






5$ Ji-J" 



5r[$] 



^^if^^ii 



5$, 



-.J 

ijl-jn- 



DiJ^LnDL (3.16b) 



(3.16c) 



Eq. (3.16b) can be derived from the equation d{M M ^] 
Mj~^ = r^-. Using these relations in eq. (3.15) yields 



0, where Mij 



Df^ and 



5r[$] 

5^i 



--s,s<^s-^su^s^t^Di)+ 



1 



- j^s,,tu{<^s<^t<^u + ^<^sDi + z^f^Ai^iriJ. 



(3.17) 



In fig. 3.1 a graphical representation of this equation is shown. Also the differentiation 



rules from eq. (3.16) can be depicted graphically as given in fig. 3.2 Since all required 



ingredients are now available in graphical form, it is very convenient to proceed like this. 
The DSE for a generic two-point function is derived by performing another differen- 



tiation of the generating DSE eq. (3.17) using the diagrammatic replacement rules of 



fig. 3.2 in the corresponding diagrammatic representation fig. 3.1 The result is shown in 



fig. 3^ Proceeding to higher vertex functions the number of diagrams increases rapidly: 
For three-point vertices there are 15 generic diagrams and for four-point functions 60. For 
real applications it is therefore recommendable to exploit possible simplifications. 



First, the final number and form of graphs depend on the first differentiation in eq. (3.17) 
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Figure 3.1: The generating DSE for IPI functions. Crosses in circles denote external fields. 
All internal propagators are dressed and the big blob denotes a dressed IPI 
vertex. The double line represents the generic field $. 

as the corresponding field determines which bare vertices appear in the diagrams. For 
example, the DSE of the ghost-gluon vertex in Landau gauge QCD has only four terms, 
when the first derivative is performed with respect to a ghost field. In this case one can 
drop all diagrams with bare gluonic vertices. On the other hand, if one starts with the 
gluon field, all vertices have to be kept and one ends up with twelve graphs. A detailed 
derivation of these two examples is provided in ref. |108| . Secondly, one can skip some 
diagrams taking into account where one is going. Simple examples are that for a three- 
point function we do not have to drag along the bare four-point vertex or diagrams with 
an external field can be dropped if no further derivatives with respect to this particular 
field follow. 

An important point is keeping the internal indices general and do not set them to specific 
values too early in the derivation process in order to try to simplify the calculation. This 
may discard some diagrams. 

The algorithm described so far did not take into Grassmann fields. As they represent 
only a technical complication and the basic structure remains the same, their inclusion is 
described in Appendix [Bj 



_6_ 



_5_ 



-r 




Figure 3.2: Diagrammatic rules for differentiating an external field, a propagator or a 
vertex. The circle with the cross denotes an external field, small blobs denote 
dressed propagators, and big blobs IPI vertices. The double line represents 
the super-field $. 
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Figure 3.3: The DSE for a generic two-point function. 
3.1.3 DoDSE: Derivation of Dyson-Schwinger equations 

Deriving DSEs can easily become very tedious: First of all, in the case of actions with a 
small number of interactions the number of terms blows up so fast that only the lowest 
n-point functions can be done by hand. Secondly, the interest in using DSEs of actions 
with many interaction terms is increasing as we get DSEs better under control. La- 
grangians with up to eleven vertices are used and necessitate a quick and easy method to 
derive DSEs. The algorithm presented above allows the implementation into a symbolic 
programming language and thereby the automated derivation of DSEs. 

For such a purpose several programming languages are available, e.g., FORM |109| or 
Mathematica |110) . The latter was chosen due to its greater accessibility, its more user- 
friendly interface and the possibility to draw the results directly as Feynman diagrams. 
The result is a package called DoDSE, which is short for "Derivation of Dyson-Schwinger 
Equations". For the beginner it is a useful tool to derive DSEs and represent them 
graphically and for the advanced user it provides further possibilities like step-by-step 
calculations, keeping external fields or using it for actions with mixed two-point func- 
tions. For simple applications, like, for example, summing up the IR exponents of a given 
graph, the primary output of DoDSE can be used for further calculations in Mathemat- 
ica. However, if one would like to obtain the integral in full glory with color and Lorentz 
indices and so on, some more effort is required. 

For this purpose DoDSE can use Feynman rules provided by the user to write down the 
complete expressions for a DSE. They can then be further processed with other functions 
as required by the user. Automating this process proved very useful for the calculation of 
the numeric values of IR exponents. While in normal Landau gauge this is a manageable 
task by hand, it can become arbitrarily complicated for other examples like the maximally 
Abelian gauge or the Gribov-Zwanziger action where there are more diagrams with more 
dressing functions to calculate, see Chapters [5] and |6l respectively, for more details. A 
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short guide to DoDSE and the derivation of the DSEs for the maximally Abelian gauge 
and the Gribov-Zwanziger action can be found in Appendix |E} 

3.2 Functional renormalization group equations 

FRGEs are used in a variety of fields in physics, for example, in the investigation of 
ultracold fermion gases, e.g., |11H 1112) . for supersymmetric models, e.g., |113till5j . in 
gravitation, e.g., |116H120| or for the phase diagram of QCD, e.g., |106l I12H - I124| . This 
list is by far not exhausted. It should only give a glimpse at the vast applicability of 
FRGEs, for reviews see refs. [Ml [T25HI27 . 

This section is devoted to the description of some properties of FRGEs and of their 
derivation. Details on how to use FRGEs in the case of interest for this thesis, the IR 



analysis of Yang-Mills theory, are provided in Section 4.1.3 



3.2.1 The effective average action 

The functional renormalization group is formulated with the help of the effective average 
action r[$]'^, which is closely related to the effective action or generating functional of IPI 



Green functions we encountered in Section |3.1.1[ The defining property of the effective 
average action T[^Y is the dependence on an artificial momentum scale fc, denoted by the 



superscript. The flow equations, derived below in Section 3.2.2 describe the dependence 
on this scale by an integro-differential equation. The quantum fluctuations above the scale 
k are all integrated out. Thus two special limits of the effective average action are A; — )■ 
and k — )■ oo, where it has to correspond to the standard effective action r[$] and the 
classical action 5, respectively: In the former case all quantum fluctuations are included 
and in the latter none. In statistical physics a natural UV cutoff A can be present. In 
this case r[$]'^ is equivalent to S for /c — )■ A. In the following, however, we will focus on 
the use of the functional renormalization group in quantum fleld theory. 

The dependence on k is introduced as an IR cutoff in the deflnition of the path integral: 

W^[J] = lnZ^[J] = In /D0e-^['^]+^'"^'-^*»^^^*^ (3.18) 

We used again the notation from the last section, i.e., writing out the integration in the 
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last term of the exponent gives 

^4'/'. = i / i^KMM-q)Mq), (3.19) 



1,0/^, 1 / d'q 

where R^^^ is an IR cutoff function with the following properties: 

• It has to vanish for A; — )■ to obtain the standard effective action. 

• It has to diverge for /c — )■ oo so that the classical action is recovered in this limit. 

• For small momenta q^ < k'^ it is proportional to k"^, which can be interpreted as an 
effective mass that constitutes an IR cutoff for fluctuations with small momenta. 

• Finally, it has to vanish for large momenta g^ > k'^ so that it does not interfere with 
the high momentum behavior of Green functions. 

The effective average action is then defined via a modified Legendre transformation: 

r'=[$] = ~W''[J] + J,<I>, - ^-K.^-K, (3.20) 

with 

*. ^ ^ ^ (*). (3.21, 

3.2.2 Flow equations 

The dependence on the scale k is used to derive the so-called flow equations or functional 
renormalization group equations. They describe how one gets from the classical action S 
to the effective action r[$] by lowering k, which is equivalent to integrating out more and 
more fluctuations. The trajectory from S to r[$] is called a flow. It can be interpreted as 
moving through theory space, starting at the microscopic theory and heading towards a 
macroscopic description. The trajectory depends on the chosen IR cutoff function R'', but 
the endpoint is always the same, because it is determined by physics. Thus approximations 
and truncations can only be checked by comparing the endpoints, not the flows. 
Differentiating r[$]^ with respect to k yields 



35 



3.2 Functional renormalization group equations 



HdkR'lj(pj 



^^^kR■,^J 



-a,4 ((0,0,) -$,$,) 



---dkR-jGfj. 



(3.22) 



where dk '■= d/dk and eq. (3.21) was used to cancel the second and third terms in the 



first hne. Furthermore, (0i0j) was decomposed as Gfj + {(f)i){(f)j) = Gf, + $i$j, where 






(3.23) 



is the connected propagator in presence of the sources J . Its inverse is the two-point 
function but with an additional contribution from I^: 

Oi^ = ^7^ = ir^ir^ — = :r^^r^ ^^ - ^ = U ,7 ["^J + ^jl G, 



S^j 5^j SJi SJi 



S^jS^i 



SJi SJi 



J 
(3.24) 



with 



and 



Ji 






S^i 



(3.25) 



(3.26) 



Thus eq. (3.22) can also be written as 

1 



a,rn 



2vrS'[*] + 4) ^^4- 



(3.27) 



From this equation functional identities for all Green function can be obtained by further 
differentiation. I will illustrate the procedure for the two-point function of a theory with 
three- and four-point functions: 



a.r 



k,J 



S^iS<^ 



-dkTm" 
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Figure 3.4: The FRGE obtained from eq. (3.28) for a theory with one field. A grey blob 
denotes the regulator insertion d^R^, a black blob a dressed n-point function. 
Internal propagators are all dressed. 
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(3.28) 



The following derivatives have been used: 



^ -pk,J 



sir^f + R^, 






jm mln nji 



■^kJ 



5^i5^j, ■ ■ ■ 5$j„ 



«J1---Jn 



(3.29a) 
(3.29b) 



Again one has to keep in mind that propagators and vertices only correspond to physical 
quantities once the external sources are set to zero at the end of the derivation. Doing so 



in eq. (3.28), the FRGE depicted in fig. 3.4 is obtained for a theory with one field. The 



inclusion of Grassmann fields goes along the same lines as for DSEs, see Appendix [B| 

Performing more differentiations with respect to fields one can derive flow equations 
for all higher n-point functions, which can also be presented graphically like DSEs. The 
appearing integrals are all one-loop and only dressed quantities appear. This is a decisive 
difference to DSEs, which contain also bare n-point functions. Furthermore, every integral 
has an insertion of the regulator R^ . Suitable choices of it for the IR analysis are discussed 
in Section [4.1.31 
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Scaling solutions 



In this chapter a method to assess the possible existence and form of a scahng solution for 
a given theory is derived. The original goal of the underlying work was to investigate the 
IR behavior of Yang-Mills theory in the maximally Abelian gauge (MAG). However, it 
soon became clear that the system of interactions given by its Lagrangian requires a more 
refined treatment than provided by any previous work, the reason being the large number 
of terms in the DSEs. Subsequently a generic method was developed that can handle a 
large set of interactions. It allows to determine for an arbitrary system of interactions 
if a scaling solution is possible and what its qualitative features are. The main proof, 
albeit relying only on simple mathematics like combinatorics and inequalities, is rather 
technical. However, it provides a simple recipe for finding possible scaling solutions. 

The underlying idea is to combine information from DSEs and FRGEs as suggested in 
ref . jHn] • Both systems of equations describe the complete content of a theory and their 
combination allows to extract enough additional information to derive all possible scaling 
relations. With the method of a combined analysis it was even possible to prove that there 
is only one scaling solution in the Landau gauge [66j. This is often called uniqueness of 
the scaling solution and should not be confused with uniqueness of the solution, see the 
discussion on scaling and decoupling solutions in Section |2.3.2 The generic method for 



obtaining the scaling relations of more general cases, originally developed for the study 
of the MAG, has been presented in ref. |67] . 

The first section of this chapter provides an introduction to the method of IR power 



counting, which is employed in Section 4.2 to derive several general statements about 



the IR behavior of dressing functions. Section 4.3 explains the method to extract pos- 



sible scaling solutions from a given Lagrangian, which is exemplified in Section 4.4 by 



employing it for linear covariant and ghost anti-ghost symmetric gauges. Up to this point 
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all considerations are for the special case of a Lagrangian diagonal in the fields at the 



two-point level. Section |4.5| provides the generalization to the case of mixed two-point 
functions. 



4.1 Infrared power counting 

In order to determine the qualitative IR behavior of correlation functions I rely in this 
thesis on a technique called IR power counting. It is based on the assumption of a 
power law behavior of all dressing functions in the IR. The qualitative information is 
then contained in the corresponding exponents and the analysis can be shifted from the 
integral equations to the exponents. In the following I will explain how they are extracted 
from DSEs and FRGEs. The role of possible or required cancelations in the functional 
equations is discussed at the end of this section. 

4.1.1 Integrals of DSEs at low external momenta 

The goal of this subsection is to explain how the IR behavior of an integral can be assessed. 
The starting point are the fully dressed propagators of the theory, parameterized by 

where P^^ is the part containing color and Lorentz structure and c*^*^ is some dressing func- 
tion. The propagators are labeled by the superscript (i). A basic assumption supported 
by renormalization group arguments is that the dressing function c^*-* obeys a simple power 
law in the IR: 



c 



».^«(p2)=rfW.(p2)M»)_ (4 2) 



Here (i*^*^ is a constant coefficient and the exponent S^'^^ is called infrared exponent (IRE). 
It describes the qualitative IR behavior of the propagator. The assumption of power law 
behavior does not take into account the possibility of a logarithmic momentum depen- 
dence if 5*^*^ is zero. One speaks of IR enhancement /suppression if S^'^^ < or 5^^^ > 0, 
respectively, and clearly the propagator vanishes if 5^^' > 1/2. The generalization to 
several dressing functions corresponding to different tensors is straightforward. 
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In the following we investigate DSEs at low external momenta. Such momenta occur 
within the integrals in factors like l/{p — k)"^, where one momentum is a loop momentum 
and the other an external one. Hence the integrand itself is dominated by low internal 
momenta and all dressing functions take their IR form in the limit of vanishing external 
momenta. Here we also assume a power like behavior for vertices in the IR with their own 
IREs. We will only treat the case that all momenta of a vertex vanish simultaneously. 
This is called uniform [73j or global scaling |66) . 

If all dressing functions are given by power laws the integral can in principle be solved 
by methods like Negative Dimensional Integration (NDIM), see, for instance, |128H137] . or 
employing the Mellin-Barnes representation, see, e.g., |138] . It is vital that these methods 
can handle non-integer exponents which is not the case for most standard methods like, 
for example, integration by parts |139l 1140) . 

The simplest case is the one-loop two-point function, where the result can be written 
down as a simple closed expression: 



(27r)''' y vv -. y v . T (-z/i) T (-z/2) T (rf + z/i + z/2) 

(4.3) 

The convergence of the integral gives constraints on the exponents Ui and z/2: UV con- 
vergence dictates d/2 + z/i + z/2 < and IR convergence d/2 + z/i > and d/2 + 1/2 > 0. 
Also for the one-loop three-point integral a solution in closed form is known |132t 1138) . 
The qualitative behavior if all external momenta scale equally can easily be assessed. A 
detailed numerical analysis, however, is complicated by the appearing Appell's function 
F4 that have several singularities in the space of the two variables given by ratios of the 
external momenta. Series representations that converge for Euclidean momenta can be 
found in refs. |79l [Tin 1112) . 

Detailed solutions for higher n-point functions can be derived and allow the assessment 
of uniform scaling. Numerically the closed forms are, however, of no use due to their 
complex form. Fortunately we only need the IREs of given diagrams. To obtain these it 
suffices in the case of uniform scaling to count the powers of all momenta in the integrals 
as upon integration the internal momenta turn into external momenta. 

In the case of higher n-point functions with n > 3 one could wonder what happens if 
only a subset of external momenta vanishes. Is it possible that additional IREs occur that 
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have to be taken into account? In principle the answer is yes. In ref. [73j it was shown 
that taking into account these additional kinematic IREs a self-consistent solution can be 
derived that is an extension of the uniform case. Numerical evidence for the existence of 
additional IREs of three-point functions can be found in ref. [7^1 1143| . where also the in 
general non-trivial dependence on the kinematics is confirmed. Interestingly additional 
divergences only appear in the longitudinal part of the three-point functions |66 | I79 | . which 
is irrelevant for Landau gauge as the set of transverse dressing functions is closed [19j. 
Also higher n-point functions do not show kinematic IR divergences in their transverse 
parts [66]. Hence kinematic divergences do not influence the IR fixed point of Landau 
gauge qualitatively, but the non-trivial kinematic dependence on external momenta can 
have a quantitative impact on quantities like k |79| . 



4.1.2 Dyson-Schwinger equations in the infrared 

Knowing how to obtain the IRE of a given diagram we can continue to the analysis of the 
DSEs themselves. The key observation is rather simple: Every diagram on the right-hand 
side of a DSE can only be as divergent as the one on the left-hand side and it cannot be 
more divergent as this would render the equation inconsistent. On the other hand one 
diagram has to scale exactly as the one on the left-hand side. This allows us to write 
down simple inequalities for the IREs. 

Let me illustrate this with a simple example. The DSE of the gluon propagator in 



Landau gauge is depicted in fig. |4.1[ With the IREs 6a, Sc, i^aaa, i^Acc and k.aaaa of 
the gluon and ghost propagators, the three-gluon, the ghost-gluon and the four-gluon 
vertices, respectively, the integrals can be evaluated for low external momentum p in four 
dimensions as 

j^-l('p2\l-5A _p2 _ / 2\l+KAAA+'iSA J^g-loop _|_ / 2\l+KAcc+'iSc j^gh-loop _ 

_ tJ^Y+i^AAAA+^SAj^sunset _ /2\1+2kaaa+4:Sa J^squint /^ ^\ 

The tadpole was not taken into account as it does not depend on the external momentum 
and can thus be absorbed in the renormalization when employing dimensional regulariza- 
tion. d-A denotes the coefficient in the power law of the gluon propagator dressing function 
and the Ls constant terms for the gluon-loop, the ghost-loop, the sunset and the squint 
diagrams. They depend on the IREs and the Ls contain also the constant coefficients 
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Figure 4.1: The gluon two-point function DSE in the Landau gauge: On its left-hand side 
the dressed gluon two-point function and on its right-hand side the bare gluon 
two-point function and the tadpole, the gluon-loop, the ghost-loop, the sunset 
and the squint diagrams. Small blobs represent bare vertices and large blobs 
dressed IPI vertices. All internal propagators are dressed. 

from the power laws. 

Every diagram on the right-hand side could be leading. Thus we can extract the 
following inequalities: 

-^A <f^AAA + 25a (4.5) 

-Sa <kacc + 24 (4.6) 

-6a <kaaaa + 3(5a (4.7) 

-6a<2kaaa + 4:6a. (4.8) 

It is worth pointing out that at this point we can get a boundary for the gluon propagator 



IRE by combining eq. (4.5) with the inequality kaaa < 0, which just reflects that the 



three-gluon vertex has a contribution from the bare vertex in its DSE. We rewrite eq. (4.5) 
and plug in the boundary on kaaa'- 

< kaaa + 35a < 35a ^ Sa>0. (4.9) 

Thus the gluon propagator cannot be IR enhanced and confinement scenarios based on a 
1/p^ behavior of the gluon propagator can be ruled out in the Landau gauge. A stricter 
argument for the non- negativity of the gluon propagator IRE can be found in ref. |73] . 
where a similar argument based on the three-gluon vertex DSE has been derived and 
possible cancelations have been investigated. 
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One could now proceed by writing down all inequalities obtained from DSEs up to a 
certain order. For DSEs with many terms this will get quite tedious and it is not sure 



how much useful information can be derived from this program. In Section 4.2.1 I will 
describe a better way of proceeding by determining the most restrictive inequalities. It 
will be found that all inequalities required for the analysis can be written down in closed 
form. 

4.1.3 Functional renormalization groups equations in the infrared 

As the IR analysis relies also on FRGEs, I explain here shortly how their IR behavior may 
be assessed. Although at the end it turns out that the IREs can be counted similarly as 
in DSEs, some special care is required as the regulator function R^ appears. 



The choice of the regulator R'^ is free, as long as the constraints listed in Section 3.2.1 
are obeyed. For an IR analysis there are some especially attractive choices which have 
been discussed in ref. [8U]. One possibility is of the form 

4(p') = r?,(p')5e(p'-fc'), (4.10) 

where 6^ is a smeared 5-functional around k'^. Such a regulator projects out modes with 
p^ ^ k"^ and the momentum dependence of physical, i.e., in the limit /c — )■ 0, n-point 
functions in the IR is not altered. 
A second possibility is 

4(p') = r?;(pXp'A')- (4.11) 

The physical two-point functions then become 

Gi/{p') = Tl{p'){l + r{pye)) (4.12) 

and again the asymptotic behavior for small p is not altered qualitatively. 

The IRE of a diagram in an FRGE is determined in the usual way by counting the 
IREs of all quantities. Let us consider as an example a diagram from a two-point FRGE 



as depicted in fig. 4.2 



(^k^ ij\P ) ^mt^ tiu^ur^ rjs^ sn^k^rnn ~^ ^ ^i-'' irjs^ sn"k^rnn' \ '-'} 
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Figure 4.2: A typical two-point FRGE. The grey blob denotes the regulator insertion 
dkR^, the black blob dressed n-point functions. Internal propagators are all 
dressed. 



As the regulator takes the momentum dependence of a two-point function, its IRE cancels 
in the power counting with that of a propagator. Hence we can just ignore regulators in 
the power counting and count one propagator instead of two propagators and a regulator. 
An alternative way of determining the IR behavior of an n-point function is given in 
ref. [66] using integrated flow equations. The flow equations take then a similar form as 
DSEs but with only dressed vertices in the diagrams. As in these equations the integral 
over the scale k is performed, also the classical term of an n-point function, defined by 
the limit k — )■ oo, appears similar to DSEs. For details I refer to ref. |66J. 

4.1.4 The role of cancelations 

Shifting the analysis from the integrals to the level of the IREs may be a delicate step. 
The reason are possible cancelations in the IR leading terms. A detailed analysis is 
necessary to assess if the IR leading terms cancel and the next leading terms take over. If 
cancelations occur they are likely linked to some symmetry and do not appear at random. 
An example is the longitudinal part of the gluon propagator in the linear covariant gauge. 
A naive power counting analysis yields for its IRE the same result as for the transverse 
part, but it is known that it does not acquire a dressing at any momentum. Details can 
be found in Section 14.4. 1[ 

The occurrence of cancelations can often not be completely disproved, but one might 
exclude some types of them. They can happen between different diagrams or between 
different contributions of one and the same diagram. These two types of cancelations can 
be ruled out for the two-point DSEs in Landau gauge, since there is only one IR leading 
diagram in each equation and the ghost-gluon vertex, the only appearing vertex in these 
equations, has only one relevant dressing function. For vertices there are always several 
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IR leading diagrams and cancelations could happen between these. Also other vertices 
except the ghost-gluon vertex with several relevant dressing functions appear. Another 
possibility is a cancelation between parts stemming from different momentum regions of an 
integral. This could happen, for instance, via a change of the sign of a dressing function. 
However, this again can be excluded in the Landau gauge: The dressing functions of both 
propagators are necessarily positive and also the ghost-gluon vertex dressing function does 
not change sign as can be inferred from lattice data [7i] - [77] and other calculations [75 | [7U ] . 
On the other hand there are cases when cancelations are required. The most prominent 
example is probably the ghost propagator of Landau gauge. In its DSE the bare two- 
point function appears that leads to the inequality 6c > and thus an IR suppressed 
propagator. As long as the bare ghost two-point function is in the DSE, one cannot 
get an IR enhanced propagator. However, it is possible to cancel the bare term in the 
renormalization by adopting the renormalization condition 

[Dc{p')pT' ^^ (4.14) 

for the ghost propagator Dc{p^)- This condition is not chosen without reason but due to 



the Kugo-Ojima and Gribov-Zwanziger confinement scenarios as discussed in Section 2.2 
Nevertheless it has to be mentioned that also a different renormalization prescription may 
be adopted which leads then to the decoupling solution, see Section 2.3[ 



The cancelation of the bare two-point function in the DSE of an IR suppressed two- 
point function is indeed crucial, as otherwise the bare term is leading. We will neglect 
in the following the bare two-point functions in the DSEs and determine a posteriori if a 
cancelation is required. In both systems investigated we will find that such cancelations 
are possible: In case of the Gribov-Zwanziger action the same reasoning as in the standard 
Landau gauge can be applied and in the MAG the renormalization of the diagonal gluon 
propagator apparently plays the same role as that of the ghost in the Landau gauge 
connecting the obtained scaling solution with the massive solution found on the lattice 
|144| and in the refined Gribov-Zwanziger scenario [IS f 1145) . 

Having clarified the role of possible cancelations it remains to be stressed that the 
analysis is exact from the point on where we go to the level of the IREs. This is quite 
remarkable as we will be able to make statements about the complete, infinitely large 
tower of integral equations without any truncations or approximations. 
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4.2 Basic relations derived from functional equations 

There is quite much information about the relations between the IREs of Green functions 
that can be extracted from functional equations. This section is devoted to the study of 
these relations. First I will show how to derive two classes of constraints for the IREs. 
Then a formula for the IRE of an arbitrary diagram will be derived. Finally I point out 
the connection between the inequalities provided by FRGEs and the skeleton expansion, 
which was employed in earlier investigations of DSEs. 

In order to keep this part as general as possible I will not refer to any special Lagrangian 
except in some examples for illustration. Therefore, generic, not necessarily different fields 
denoted by A, B, C, . . . are used. The only restriction on them is that they are of mass 
dimension {d — 2)/2. Otherwise the canonical dimensions of correlation functions will 
change. Note that in general no real simplifications occur if the equations are written 
down for a definite Lagrangian. Even worse, in many cases writing out all the equations 
would lead to extremely large expressions the treatment of which is not feasible. The 
MAG, which is the topic of Section [5| is a perfect example for such a case. In this and the 
following sections all expressions are only valid for four dimensions except when denoted 
otherwise. The corresponding expressions valid in d dimensions are given in App. |A} to 
which all detailed calculations of this section are deferred. 

4.2.1 Inequalities from FRGEs and DSEs 

At first sight it seems that DSEs and FRGEs allow to derive a plethora of inequalities as 



described in Section 4.1.2 However, on closer inspection one can see that a part of them 
is contained within others. Therefore we should look for the minimal set of inequalities 
required. It turns out that there are two groups which play a distinguished role. The first 
connects the IRE of a vertex and the IREs corresponding to its legs, and the second relates 
the IREs of propagators among themselves. The first group consists of infinitely many 
inequalities and the second of only a few. In fact we will see that there are at most only 
as many members of this group as there are primitively divergent vertices, i.e., vertices 
appearing in the Lagrangian. This is the decisive aspect that makes the IR analysis a 
manageable task. 

Before we start a few explanations on the notation are required. Throughout this thesis 
I use the letter k for IREs of n-point functions and the letter 6 for IREs of propagators. 
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Figure 4.3: One specific diagram in the FRGE of a generic tliree-point function. Internal 
lines represent dressed propagators, black blobs dressed vertices. The grey 
blob is a regulator insertion. 

The reason for this is the discrimination of two-point functions and propagators as it 
will be required in Section |6j The subscripts of the Greek letters denote the fields of 
the external legs. In the case of propagators and two-point functions the abbreviations 



644 and /€,■ 



are used, respectively. Furthermore, vertices are abbreviated by 



monomials of their fields, e.g., ABCC corresponds to a four-point function with one A, 
one B and two C legs. Such a monomial can also denote a DSE. Then the order of 
fields refers to the order in which the derivatives are applied as this can lead to different 
realizations of the equation. The field (pi corresponds to that field indicated by its index 



(see also Section 3.1.1). 



We start with three-point functions. In their FRGE there appears a triangle diagram 



as given in fig. 43 The corresponding inequality is 



KABc + ^{Sa + 6b + 6c)>0. (4.15) 



This equation already gives a hint at the general form of these inequalities: Given the 
IRE of a vertex, we can add half the sum of the propagator IREs corresponding to the 
legs and we get a non-negative quantity: 



K. 



....v + ^E^"''"'^^^o- 



(4.16) 



The symbol k^' denotes the number of times the field 0j appears in the vertex 
It remains to be shown that this inequality is true for all n-point functions. 
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"m 




b) 





Figure 4.4: Parts of the FRGEs of generic four-point functions. 



For the general four-point inequahty corresponding to eq. (4.15) we have to make a 



detour via a special case. The FRGEs from which we infer these inequalities are given in 



fig. 4.4 They read 



-Sa — 5 b < f^AABB, 

f^AABB < '^f^ABCD + 5c + Sd- 



(4.17) 
(4.18) 



Combining them yields 



kabcd + ^{Sa + Sb + Sc + Sd) > 0, 



(4.19) 



what is the expected inequality for a four-point function. It should be noted that the 
inequalities of this group can only be derived from FRGEs, since the appearance of the 
bare vertex in DSEs leads to a different, less restrictive numerical factor in front of the 



propagator IREs in eq. (4.16) 



In a similar way one can prove by induction the validity of eq. (4.16) for higher n-point 



functions. I give here only the general structure of the proof and refer for details to 
Appendix |X| For the proof one needs two relations which can be inferred from fig. 4.5 
The figure on the right-hand side relates two n-point functions with an m-point function 
where every field appears twice and m = 2n — 4. The connection of this m-point function 
to an (m-2)-point function is given in the left figure. Plugging the one inequality into the 
other establishes an inequality between an n-point and a (2n-6)-point function. For the 
latter one can repeatedly use relations obtained from the type of FRGE depicted on the 



left-hand side of fig. 4.5 until one gets to the equation as given in fig. 4.4 on the left. Thus 



one can establish eq. (4.16) for all n-point functions. 



Eq. (4.16 ) allows an interesting observation concerning the so-called skeleton expansion. 



48 



4.2 Basic relations derived from functional equations 



k 



dk 





A / \ D A / \ D 

a) B C B C 



b) 





Figure 4.5: Parts of the FRGEs of generic n-point functions. 



The corresponding discussion is deferred to Section 4.2.3 

Now we turn to the second group of inequahties. They can be derived directly from 



eq. (4.16) by noting that for all interactions appearing in the Lagrangian, the primitively 



divergent vertices, the corresponding bare vertices appear in their DSEs. Thus their IRE 
is non-positive, 



prim. div. ^ p. 



(4.20) 



and we infer from eq. (4.16) that 



- 2^ kl^"'^''5i > V primitively divergent vertices. 



(4.21) 



This group of inequalities can be derived purely from DSEs, but not from the FRGEs 
without help from the DSEs. 

In the preceding analysis one should always keep in mind that for a realistic Lagrangian 
it is not sufficient to replace the generic fields A, B etc. by the fields of the Lagrangian, 
but one has to check if the corresponding diagrams really exist. The reason for the 
disappearance of a diagram can be as simple as a trivial color factor or as complex as a 
non-manifest symmetry. The former case will be encountered in the analysis of the MAG 
but only for diagrams not relevant in the IR. Thus they do not invalidate the result. 



4.2.2 Infrared exponent of an arbitrary diagram 



We established in Section AA_ how the IRE of any given diagram can be evaluated by 
power counting. This is now formulated in very general terms by considering an /-loop 
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diagram v with n^ internal propagators of the field (pi, rrii external legs of the field 0j, 

■■4>k- 



^t...k tiare vertices of the type (pi- ■ ■ 4>k and nf j, dressed vertices of the type 



K,; 



', 


+ >> 

i 


i(«.- 


1)+ 


vertices,r>3 


<.... 


+ 


>: 


"i 


..iM- 


.ir ^V- 





\'^il...ir ' ^il...ir )' 



(4.22) 



vertices,r>3 



Here Ci,,,j/K,i,,,j denote the canonical dimension/ IRE of the vertex (pi- ■ ■ (pj. Eq. (4.22) is 
valid for d dimensions, but in the following I only give the results for d = 4, since this 
makes the expressions more transparent. The arguments, however, are the same and the 



details for d dimensions can be found in Appendix A. 2 The double subscripts of the n"^ 
and n^ indicate all possible combinations of r fields. E.g., for the Landau gauge the term 
corresponding to r = 3 is 



^AAA ( I^AAA + 2 



(> 



+ ^A-cc \ I^Ac^ 



cc 



1 

H — n 



AAA + 2^^cc- 



(4.23) 



Using topological relations that connect the numbers of propagators, vertices, loops 



and legs eq. (4.22) can be rewritten: 



K,; 



•i vertices,r>3 \ i / 



vertices,r>3 



(4.24) 



Details on the derivation of this formula are given in Appendix A. 2 The reason for 
rewriting it in this way is that now it depends only on the external legs and the numbers 
of vertices and we got rid of the dependence on internal propagators. 



Eq. (4.24) is not only useful in general, but it allows also to determine a lower bound 



for the IRE of a diagram. To this end we resort to the results of Section 4.2.1 which 



are given by eqs. (4.16) and (4.21). Comparing them to eq. (4.24) we see that they 



correspond exactly to the terms in parentheses. Even the fact that eq. (4.21) is only 



valid for primitively divergent diagrams is mirrored as the corresponding coefficient in 
eq. (4.24), ra^^ j^, can only be non-zero for primitively divergent vertices. Hence the first 



term in eq. (4.24) is a lower bound for the IRE of the diagram v due to the non- negativity 
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of the other two expressions. Even more, since the first term only depends on the external 
legs it is the same for all diagrams appearing in a DSE. This establishes a lower bound for 
the IRE of an arbitrary vertex function. It is called the maximally IR divergent solution: 



__lv- . 



(4.25) 



4.2.3 Skeleton expansion of higher vertices 

Before continuing with the IR analysis I want to make a short detour to a tool that 
was used repeatedly in IR analyses of Green functions: the skeleton expansion. It can 
be seen as a loop expansion of a vertex using only dressed propagators and vertices; for 
more details on the original idea see, e.g., ref. |146] . In functional equations one usually 
restricts the set of vertices in the expansion to those appearing in the Lagrangian. One 
can then derive the IREs of all vertices in Landau gauge [T^ I147 J. The assumption that 
the skeleton expansion does not explicitly diverge, meaning that higher orders do not 
feature more IR divergent diagrams, provides sufficient additional information to prove 
the uniqueness of the IR scaling solution of Landau gauge [73]. The key observation 
thereby is that higher orders of the skeleton expansion are constructed by adding further 
loops using certain combinations of correlation functions. These insertions into diagrams 
must not make the diagrams more IR divergent and hence they must have non-negative 
IREs. 

Originally the validity of the skeleton expansion was an assumption, but the functional 
renormalization group clarifies how and why the skeleton expansion indeed works for 



power counting. The explanation is given by eq. (4.16), which was derived from FRGEs. 
The inequalities provided by it correspond exactly to the inequalities derived from the 
insertions used in the skeleton expansion, see ref. [73] for details. Thus the use of the 



skeleton expansion is well justified. However, eq. (4.16) is even more powerful, because 
it covers all possible vertices in contrast to the skeleton expansion, which could - at least 
in the form as it was used - provide inequalities only for the vertices appearing in the 
Lagrangian. 
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4.3 Solution for the system of IREs 

The information collected up to now will be used to determine the scaling solution(s) of 
a given Lagrangian. For this purpose we analyze the DSEs of the two-point functions. 
However, instead of considering single cases we again keep things general. The analysis 
proceeds by assuming that one diagram in a two-point function DSE is leading without 
specifying which one. This will lead to a new inequality that has only to be fulfilled in 
this one case. Taking into account the inequalities obtained earlier, we can discard some 
solutions and come up with a set of possible IR scaling solutions for the given system. 
This set can be empty, consist of a unique solution or of several different ones. 

4.3.1 Analysis of two-point equations 



We consider a generic two-point function DSE. Fig. 4.6 illustrates possible diagrams if 
only three- and four-point interactions appear in the Lagrangian. In the following we 
do not refer to specific diagrams and the analysis covers all possible diagrams. For the 



leading diagram the following equation can be derived from eq. (4.24) 



«i = -^$^"^i'^.+ Yl '^n.-.^J ^J^^r-"''^.] + Yl <...v'^n...v, (4.26) 

i vertices \ i / dressed 



vertices 



where Kj is the IRE of the two-point function under investigation. The only values de- 
pending on the specific leading diagram are njj.,.j^ and nf j^, but we can leave them 
arbitrary for now. Using again some topological relations and the lower bound for the 



IRE of a vertex, eq. (4.25), we arrive at the inequality 



<...^r\lT.'^^^'"'^] ^0- (4.27) 



The details of the calculation can be found in Appendix A. 3 The only information from 
the leading diagram that enters here is the type of its bare vertex as given in n^ ^ , i.e., 
we get as many different possible inequalities as there are bare vertices. 



Interestingly eq. (4.27) looks exactly the same as the inequalities given in eq. (4.21 ), but 



with the opposite sign. Hence, both inequalities can only be fulfilled if they are saturated. 
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{a) 




(b) 





id) 




Figure 4.6: Topologies of diagrams appearing in the DSE of a two-point function for three- 
and four-point interactions only. Conventional names for the diagrams are (a) 
tadpole, (c) sunset, (d) squint. 



I.e., 



2 '^■ 



J25,kf--^ = 0. 



(4.28) 



The hat was added in order to indicate that this equality was only derived for one specific 
interaction, namely the one appearing with a bare vertex in the IR leading diagram. From 



eqs. (4.16), (4.20) and (4.28) it follows directly that the IRE of the corresponding vertex 



IS zero: 



Ml...Jr 



0. 



(4.29) 



Thus for a scaling solution at least one vertex does not get IR enhanced. This vertex is 
called the IR leading vertex for reasons that will become clear soon. 

In order to obtain valid scaling relations one has to determine all possible realizations 



of eq. (4.29 ) and determine if they are compatible with the residual system of inequalities. 



In some cases the scaling relation does not relate all propagator IREs to each other and 
a second scaling relation has to be obtained from the remaining DSEs. Very often this 
leads to the trivial solution where all IREs are zero. Known examples where one sector of 
the theory possesses a scaling relation independent of a second sector are massless QCD 
|148j and a fundamentally charged scalar coupled to Yang-Mills theory |149l 1150) : The 
Yang-Mills sector is not influenced by the quark/scalar so that the usual scaling relation is 
valid. As the quark/scalar does not couple directly to the ghost, which is the IR dominant 
field, the quark/scalar propagator is not IR suppressed and its IRE is zero. 

It has to be stressed that the notion of a not IR enhanced vertex is not tantamount 
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to a bare vertex in the IR. In general a vertex has several tensors with corresponding 
dressing functions and the IR analysis just takes into account the most IR divergent 
one(s). There can always be other dressing functions with a smaller IRE. For the numerical 
determination of an IRE it suffices to use only the IR leading tensors, but for a complete 
numerical solution of the equation also the subdominant tensors have to be taken into 
account. If a vertex has a bare counterpart and is not IR enhanced, then the bare vertex 
will be part of the IR leading terms unless it gets canceled. 

For Landau gauge the issue of choosing the IR leading tensor is easy, as its IR leading 
vertex, the ghost-gluon vertex, only possesses two tensors of which only one is transverse. 
The longitudinal tensor does not contribute in Landau gauge and hence is irrelevant |19j . 
The impact of different dressing functions on the value of the parameter n was investigated 
in detail in ref. |65]. Thereby several constraints, like the unchanged uniform scaling, 
were taken into account with the result that k, only changes slightly. Additionally the IR 
dressing function was calculated in ref. |79) as a function of the external momenta, where 
the solutions for the propagators and a bare ghost-gluon vertex were used as input. Even 
in this simple case the dressing function depends on the kinematics. Thus a more detailed 
calculation of n should take into account the momentum dependence of the ghost-gluon 
vertex. However, the qualitative behavior of the propagators would not be affected by 
this. 

Having determined the set of possible scaling solutions we have to identify the leading 



diagrams in the two-point DSEs. To each solution corresponds via eq. (4.28 ) a bare vertex. 
Since all dressed vertices and propagators dropped out in its derivation, all diagrams with 
the same bare vertex lead to the same scaling relation. However, a word of caution is in 
order here: This does not necessarily mean that all diagrams with the same bare vertex 
scale equally. In the case of the MAG such an instance is encountered and will be discussed 
in detail in Section I5.3.2[ 

We can summarize the procedure to obtain the possible scaling solutions of a theory as 
follows: 

1. Determine the inequalities that are derived from the interactions appearing in the 



Lagrangian from eq. (4.21). 



2. Reduce the number of inequalities, if some of them are contained within others. 



3. Try to saturate one inequality after the other according to eq. (4.28) and see, if you 
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can find any contradictions with the remaining inequahties. 

Let me illustrate the simpleness of the method just described by showing how it works 
for Landau gauge: 

1. There are three interactions in the Lagrangian: the ghost-gluon, the three-gluon and 



the four-gluon vertices. According to eq. (4.21) they lead to the three inequalities 



Isa + 6,>0, ^Sa>0, 25a >0, (4.30) 

where 6a and 6c are the IREs of the gluon and ghost propagators, respectively. They 
are all valid at all times. 

2. We reduce this to the two inequalities 

16a + 4 > 0, 6a > 0. (4.31) 

3. a) Saturating the second inequality yields 6a = 0, which in turn gives a non- 

negative IRE of the ghost propagator: 6c > 0. The equation 5^ = is obtained 
from diagrams with bare three-gluon or four-gluon vertices. As such diagrams 
do not appear in the ghost two-point function DSE, we still have to determine 
its IR leading diagram. Based on the general argument given above we use 
that the only option for a scaling relation involving the ghost propagator IRE 
is 6a + 26c = 0. Consequently we find 6c = 0, i.e., this case corresponds to the 
trivial solution. 

Since the ghost two-point function DSE in the Landau gauge as depicted in 



fig. 4.7 has only two diagrams, one could also investigate it directly: For the 
bare two-point function we find 6c = without any effort. For the loop diagram 
we get with the analysis described above 6a+26c = 0, which also leads to 6c = 
for 6a = 0. 

b) The first inequality corresponds to the case when the ghost-gluon vertex is the 
IR leading vertex. Hence the scaling relation is 6a + 26c = and the ghost-loop 
and the mixed loop in the gluon and ghost two-point DSEs, respectively, are 
IR leading. Conventionally the parameter k := —6c is used. 
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Figure 4.7: The ghost two-point function DSE in the Landau gauge. 



The generahzation to d dimensions is straightforward. The only difference is that the 
canonical dimensions do not factor out so that additional terms proportional to d/2 — 2 
appear. Appendix |A] contains all required formulae in d dimensions. For the scaling 



relation we use eq. (A.25): 



5a + 25c+^-2 



0. 



(4.32) 



Following a simple list of instructions we could derive the known scaling relation of 
Landau gauge [T51IT6]. What remains to be determined is the IR behavior of the vertices. 



4.3.2 Infrared exponents of vertices 

From the IR leading diagrams in the two-point DSEs one can construct the IR leading 
diagrams of higher vertices by adding more legs via insertions of the IR leading vertex. 
For the higher vertices the corresponding IRE is given by the maximally IR divergent 



solution, eq. (4.25), since we only add vertices for which the other two terms, as appearing 



in eq. (4.24), are zero due to eqs. (4.28) and (4.29) 



As an explicit example we consider again the Landau gauge. From the arguments above 
the solution for a vertex with m gluon and 2n ghost legs follows as 



^m,2n 






{n — 'm)K, 



(4.33) 



which was found by a different route in ref. |72) . 



For the generalization to d dimensions one plugs the scaling relation given in eq. (4.32) 



into the maximally IR divergent solution in d dimensions as provided by eq. (A. 19): 



m 



1 



m 
~2 



n 



(n 



m)K + (1 



n] 



d 



- 2 . (4.34) 



This is in agreement with an earlier investigation |147l 1151] . 

The Landau gauge represents the ideal case, where all vertices can be constructed from 
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an IR leading diagram of a two-point DSE. However, it can also happen that this is not 
possible, e.g., if the IR leading vertex has four legs. How would one construct then a 
three-point function from a two-point function by inserting four-point functions? The 
answer is one cannot do so, as can be inferred by simple combinatorial considerations: 
Out of vertices with an even number of legs one cannot construct one with an odd number 
of legs. 

4.4 Examples of infrared scaling relations 

The method described in the previous section can directly be used for the investigation 
of different gauges. Two gauges investigated in the past besides the Landau gauge are 
linear covariant gauges and ghost anti-ghost symmetric gauges [152J. In their analysis a 
truncation scheme with bare vertices was used |152) . For linear covariant gauges a non- 
trivial scaling solution was only found, when the longitudinal part of the gluon propagator 
scales equally as the transverse part and for ghost anti-ghost symmetric gauges only the 
trivial solution was found. I generalize these results here by including dressed vertices, 
but the main outcome of ref. |152| is corroborated. 

4.4.1 Infrared analysis of the linear covariant gauge 

As described in Section [2. 1 . 2| linear covariant gauges can be considered as a generalization 
of the Landau gauge by relaxing the restriction to the hyperplane 9^v4^ into a Gaussian 
distribution on the gauge orbit. The corresponding gauge fixing part of the Lagrangian 



is given in eq. (2.22). An immediate consequence of this gauge choice is that the gluon 
propagator is no longer purely transverse. 

A peculiar feature of the linear covariant gauge is that the longitudinal part of the 
gluon propagator does not acquire a dressing but stays proportional to the gauge fixing 
parameter ^: 



Dt.ip') = 5"' (9,. - ^"j ^ + e-^. (4.35) 



This can be inferred either from a Slavnov- Taylor identity [U [T7] or from the fact that 
the longitudinal part of the gluon propagator, {{d^A^){d^A^)), corresponds to the second 
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moment of the Gaussian distribution 

and is hence proportional to ^ [153]. 

To include the longitudinal part of the gluon into the formalism developed so far we 
interpret it as an additional field. Its propagator has the IRE 6A,iong- This allows to 
distinguish also for the vertices between transverse and longitudinal parts. Naively the 
same inequalities are obtained for the longitudinal gluon field as for the transverse one. 
This could change if some vertices vanished when they are contracted with a longitudinal 
gluon propagator. However, as the Landau gauge result showed the ghost-gluon vertex to 
be the IR leading vertex, we obtain for the longitudinal gluon propagator IRE the same 
result as for the transverse one, as the longitudinal part of the ghost-gluon vertex does not 
vanish. The result 6a = ^A^iong = — 2(5c leads together with the triviality of the longitudinal 
part, 6A,iong = 0, to 6a = —'26c = 0. This allows two possible conclusions: Either there 
is no non-trivial scaling relation for linear covariant gauges, or the naive application of 
the IR analysis is insufficient here. Most probably there are some cancelations related to 
the triviality of the longitudinal part of the gluon propagator. How these can be made 
manifest within DSEs is currently unknown. 

4.4.2 Infrared analysis of the ghost anti-ghost symmetric gauge 

The most general Lagrangian of dimension four constructed from gluons and Faddeev- 
Popov ghosts which is Lorentz invariant, globally gauge invariant, HermitiaiQ BRST 
invariant and anti-BRST invariant without topological terms is |152[ I154[ 1155] 

^,a =\fi + l-^{d,A,f + f (i - f ) \r'^r'^c'c^c'^c^+ 



{D^c)d,c - (l - {d^c)D,c. (4.38) 



^For the hermiticity of the Lagrangian the following transformation properties of the ghosts are chosen 



.t _ ^ ?!■__;; (4,37) 



c = c, c = —c. 
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The gauge fixing parameter p allows to interpolate between the Landau gauge (p = 0, 2) 
and the ghost anti-ghost symmetric case (p = 1). 

The IR analysis of this Lagrangian can be done along the same lines as for the Landau 
gauge. The vital difference between the two gauges is clearly the occurrence of a quartic 



ghost interaction, which leads to a new inequality due to eq. (4.21): 



4 > 0. (4.39) 

It is now straightforward to conclude that only the trivial solution can be obtained. The 
inequalities from the vertices in the Lagrangian are 

5a > 0, 5c > 0, 5a + 25c > 0. (4.40) 

If the first inequality is saturated, i.e., 6a = 0, the IR leading vertex is purely gluonic and 
we still need to determine the IR leading vertex of the ghost two-point function DSE. In 
both possible cases, quartic ghost vertex and ghost-gluon vertex, we obtain 6c = and 
thus the trivial solution. If we start with the inequality 6c = 0, the arguments are the 
same. Finally, for 6a + 26c = we immediately get 6a = 5c = 0, &s both 6a and 6c are 
non-positive. Hence all possibilities lead to the trivial solution. 

This confirms the result from ref. |152) . where bare vertices were employed, also in the 
case of dressed vertices. Again we can only conclude that the existence of a scaling solution 
is not possible or that more effort is required to expose some not manifest symmetries 
which lead to cancelations. 

The purpose of this section was to demonstrate how the method for obtaining scaling 
relations facilitates the analysis compared to earlier investigations. Before we turn to 
systems which have not been investigated yet in Chapters [5] and |6| I will generalize the 
method to the case of Lagrangians where the fields mix at the two-point level. 

4.5 Extension to actions with mixed two-point 
functions 

The method developed in the preceding sections works very nicely for many actions. There 
are, however, also Lagrangians where the fields mix on the level of two-point functions. 
In this case additional complications arise, because a propagator is no longer the inverse 
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of the corresponding two-point function. In fact the connection between propagators and 
two-point functions becomes a matrix relation. This requires a more detailed analysis to 
determine the relation between their IREs. The identity 6i = —Hi, valid before, does only 
hold as a special case and in general the IREs of a propagator and the corresponding two- 
point are related by an equation involving also the IREs of other propagators/two-point 



functions. Thus some formulae of Sections 4.2 and 4.3 have to be generalized. 



The first important observation is that the inequalities (4.16) and (4.21) obtained from 



FRGEs and DSEs do not change, the reason being simply that in the case of mixing two- 
point functions the number of diagrams in the equations grows, but all diagrams from 
before will still be there in general. Furthermore, the new diagrams do not lead to new 
relevant inequalities. 

The formula for the IRE of an arbitrary diagram, however, is altered by the appearance 
of the new propagators. A detailed derivation is given in Appendix A. 3 For the sake of 
simplicity I only give the result if there are two mixing fields, denoted by A and V. The 
final expression is 



K^ 



i=A,V 
vertices,r>3 \ i=A,V / 



vertices,r>3 \ i=A,V / 



(4.41) 



The new term contains uav, which is the number of mixed propagators, and the combi- 
nation Aav '■= "26 Av -Sa-^v of the IREs of the propagators. The value of the latter has 
to be considered for each case. It should be non-negative or the IRE is unbounded from 
below, as one can find diagrams with an arbitrary number of mixed propagators. They 
are explicitly constructed by appropriately replacing dressed quantities by their respective 
DSEs. 

Finally, we have to redo the analysis of the leading diagram in two-point function DSEs. 



The resulting inequality corresponding to eq. (4.27) also contains a new term proportional 

to Aav- 



'^i+oJ2^: 



3^3 



n, 
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1 



/ 



{25av - Sa- S 



V 



nAv+ ^ 



l\...lr 



^il...ir^AV 



>0. 



(4.42) 



V 



dressed 
vertices 



The symbol fc^y *"■ denotes the number of times a mixed propagator is contained in the 
IR leading diagram of the vertex ipi^- ■ ■ 0v ■ 

In the remaining analysis ^av plays a special role. If it can be shown to be zero, all new 
terms vanish and one can proceed as before. If it is not zero, more care is required. Since 



an analysis only makes sense for positive values, it is clear from eq. (4.41 ) that the leading 



diagrams are those with the lowest possible number of mixed propagators, because for 
every one of them the IRE is raised. Hence one determines first the IREs of all diagonal 
two-point functions and subsequently those of the mixed ones. An application of this 



procedure can be found in Section 6.2.2 
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Chapter 5 

The infrared regime of the maximally 
Abelian gauge 



For quantizing Yang-Mills theory naturally a variety of gauges is available, each with 
its own advantages and disadvantages. Depending on the aspect one is interested in 
one chooses one gauge or another: Not only that the complexity of the calculations is 
influenced by such a choice, but also some properties are more accessible in certain gauges. 
In this chapter I discuss the maximally Abelian gauge, which is especially amenable for 
investigating the role of Abelian field configurations. These are of relevance for the dual 
superconductor picture of confinement, see, for example, pll l25l I156H159] . I will first 
give a short overview of this confinement scenario and the related hypothesis of Abelian 
IR dominance |160] . Then I will introduce the maximally Abelian gauge [23l 1161] and 
use the method devised in Chapter |4] to obtain its scaling solution, which supports the 
hypothesis of Abelian IR dominance. These results have been published in refs. p7 tll62] . 
I also present additional explicit numerical calculations yielding values for the IREs. 

5.1 The dual superconductor confinement scenario 

The absence of free chromoelectric charges led to the idea that the QCD vacuum may 
be a sort of superconductor |l^^ [25] . In a conventional superconductor of type I the 
magnetic field is expelled due to the Meifiner-Ochsenfeld effect except for a thin layer 
below the surface. If, however, the magnitude of the external magnetic field is raised, 
superconductivity breaks down and there is a sharp transition to the non-superconducting 
phase. For a type-II superconductor there is an additional phase in between, where there 
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is still zero resistance, but the magnetic flux can penetrate the material by means of flux 
tubes, also called Abrikosov vortices. This mixed phase is characterized by two critical 
field strengths: The first one gives the field strength where the magnetic flux starts to 
penetrate the material and the second gives the value when superconductivity breaks 
down. 

One might think of the QCD vacuum in a similar fashion, namely being a dual su- 
perconductor where not the magnetic but the chromoelectric field lines are confined into 
vortices. This picture provides a natural explanation for the formation of flux tubes be- 
tween quarks. A possible reason for the formation of such vortices is the condensation of 
chromomagnetic monopoles, which can be identified after reducing the gauge symmetry 
from SU{N) to f/(l)^^^ by choice of an Abelian gauge [23]. The most prominent one 
is the MAG as introduced in ref. [23] which makes the gauge field as diagonal as possi- 



ble, because the diagonal part is identified as the Abelian part, see Section [5. 2. 1[ Hence 
one speaks of the diagonal and the off-diagonal parts of the algebra. Indeed in early 
lattice calculations the condensation of magnetic monopoles was observed in this gauge 
|1631I164] . Alternative possibilities for the realization of a dual superconductor exist if no 
magnetic monopoles occur, e.g., in the Landau gauge, where also a dual Meifiner effect 
was observed |165) . 

Directly connected to the dual superconductor scenario is the hypothesis of Abelian IR 
dominance |16U| . Based on the assumption that chromomagnetic monopoles cause con- 
finement Ezawa and Iwazaki argued that the Abelian parts of the gauge fields should be 
dominant at large distances, since classic magnetic monopoles live in the Cartan subalge- 
bra [TBDJ. This is the algebra constructed from the maximal set of commuting generators 
and is thereby directly related to the Abelian part of a gauge field. The IR dominance of 
Abelian configurations in turn means that off-diagonal fields do not propagate over large 
distances. In subsequent years this was usually attributed to the occurrence of a mass 
for off-diagonal fields due to which these fields decouple according to the Appelquist- 
Carazzone theorem |166| . In the present work another another possibility is realized, 
namely an IR suppression of the off-diagonal propagators without a mass. The IR dom- 
inance of the Abelian degrees of freedom is realized by an IR divergent diagonal gluon 
propagator. 

Although never explicitly proven indications of Abelian dominance were found in many 
instances. One is the fact that the string tension between a static pair of quarks has 
almost the same value if calculated from the diagonal part of gauge fields alone as from 
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the complete fields, see, for example ref. |167| . Related to this is the notion of monopole 
dominance, which means that the string tension calculated from the monopole part of the 
diagonal fields is again close to the true value |168l 1169] . 

There are also lattice simulations that calculated the propagators in the MAG. This 
allows direct comparisons to the results of functional methods. Early calculations can be 
found in refs. |170l 1171) . The most recent data with the lowest momenta available so far 
is in refs. |144l I172J . The difference to Landau gauge results is immediately obvious as 
all three propagators (diagonal gluon, off-diagonal gluon, ghost) become finite at vanish- 
ing momentum. Furthermore, there are less statistical errors. The results support the 
hypothesis of Abelian IR dominance as the propagators of the off-diagonal gluon and the 
ghost fields are suppressed at low momenta compared to that of the diagonal gluon field. 

5.2 The maximally Abelian gauge 

This section describes how the Abelian part of the gluon field is identified and how the 
gauge is fixed in the MAG. The gauge fixing procedure is done in such a way to allow a 
continuous interpolation to the Landau gauge. From the action of the MAG the DSEs of 
the two-point functions are derived. 

5.2.1 Identifying the Abelian part of the gauge field 

One way to determine the importance of the Abelian part of Yang-Mills theory is to use 
the freedom to choose a gauge. By minimizing the norm of the off-diagonal gauge fields 
one can emphasize the role of the diagonal part. For this it is necessary to identify the 
diagonal part of the gauge field what is done via the generators of the corresponding 
algebra. The generators of an algebra T*" are defined by 

[r,T']=ifT\ (5.1) 

where /^** are the structure constants. In the following we employ a widely used conven- 
tion where the indices of Abelian generators (and corresponding fields later on) are given 
by i, j, ... and those of the other generators by a, b, . . .. Furthermore r, s, . . . represent 
both indices. The Abelian part of an algebra, the so-called Cartan subalgebra, is defined 
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SU{2) 


SU{N > 2) 


pk 








pa 








r'' 
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fabc 





/ 



Table 5.1: Overview over structure constants with diagonal 



.h J, 



indices (a, 6, . . . ). 



and off- diagonal 



by the maximal set of generators that commute among each other: 



[T\T^] =0. 



(5.2) 



The matrices for the Abelian generators can be chosen diagonal so that one speaks of the 
diagonal part. The diagonal generators are given in SU{N) by 



2^j 



J(J + 1) 



1/2 



X diag{l,...,l,-j,0,...,0), j 

j times 



.N~l. 



(5.3) 



The other generators are off-diagonal matrices and so one denotes the corresponding 
fields as off-diagonal. This nomenclature avoids confusion about an Abelian field in a 
non-Abelian field theory, albeit sometimes in the literature the term non-Abelian field is 
used for the off-diagonal field alone. 

For later it is convenient to determine the non-zero structure constants. Clearly those 



with at least two diagonal indices vanish due to eq. (5.2). If only one index is diagonal 



it is non-zero, but if all three indices are off-diagonal, one has to consider the specific 
algebra. Only if it possesses three off-diagonal generators it can be non-zero due to the 
antisymmetry of the structure functions. For SU{N) this is always the case except for 
A^ = 2, as the number of off-diagonal generators is A^^ — N = 2. For SU{2) the only 
possible combination is then two off-diagonal and one diagonal indices. This leads to a 
simplification of a Lagrangian based on SU{2) compared to higher SU{N), when diagonal 
and off-diagonal parts of the fields are treated separately. Tab. 5.1| gives the vanishing 
structure constants. 

The gauge field A^ is given by its components in the Lie algebra AJ^, so its natural to 
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separate the diagonal and off-diagonal parts as 

A^ = TM; + T^B"^. (5.4) 

In addition to the index convention the off-diagonal gluon field B was given a new symbol 
in order to alleviate the distinction between the two gluon fields in the following. These 
two fields may behave quite differently, what can already be inferred from the Yang-Mills 
Lagrangian, which is split in the same fashion: 

Sym = J dx\F;,F;, = Jdx Qf;,f;, + ^f;,f;,^ (5.5) 

with 

^;. =d,B^u - duB; - g r'^B^B: - g r'^B^Al - g r'A^B^ = 

=DfBl - DfBl - g r'^B^Bl (5.6) 

F;, =d,A', - d.Al - g r'BlBl, (5.7) 

Df:=5'''d, + gr''A'^. (5.8) 

Explicitly the diagonal and off-diagonal terms in the Yang-Mills Lagrangian read as fol- 
lows: 

F;.F;^ = 2{d,Bl){d,B:) - 2{d,Bl){d,B^^) + 

- 2g r'^BlBliid.Bl) - {d.B;)) - Ag r'^ BlAl{{d,B:) - (9.S;)) + 
+ 2/ r'^r'^Bl.AliBiAi - Aj^Bt) + 2/ r'j'^'' B'^B^iB^Al - A^^Bt) + 
+ g'r''r'''BlBlBlBl (5.9) 

F;,F;, = 2{d,A',){d,Al) - 2{d,A',){d^Al)+ 
-2gr'B;Bl{{d,A^,)-{d^A\)) + 
+ g'r'r''B;BlBlBt. (5.10) 

In summary there are BBBB-, ABBB-, AABB-, ABB- and 555-interactions. 

It should be stressed that for later convenience the covariant derivative was defined 



in eq. (5.8) only with the diagonal gluon field and not as its projection from D'^J' 



S^^dfj, + g f^ A , where also an off-diagonal field would appear. At this point it does not 
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make a difference, yet, as it acts on an off-diagonal gluon field and the additional term 
would be zero: f^^^B^^B^ = —f'^^'^B^B^^ = 0. However, if it acts on a ghost field, this term 
has to be taken into account. Accordingly we define the "full" covariant derivative as 

D;' := 5^% + g r'^Bl + g t^^A^. (5.11) 

5.2.2 Gauge fixing and renormalizability 

Next we will determine the gauge fixing part of the action. For this we need to derive 
the BRST transformation, where we discriminate again between the diagonal and off- 
diagonal parts. Furthermore the Jacobi identity splits into several identities. The modified 
expressions can be directly derived from the usual ones by ascribing the free indices to 
either the diagonal or off-diagonal sector. As an example we have a look at the BRST 
transformation of the gluon field, given by 

sA; = -D;V. (5.12) 

Choosing r to be either a diagonal or off-diagonal index, i and a, respectively, we get the 
BRST transformation for the two gluon fields: 



s B; = -Dfc' - b^d = -{D^'c" - g r'-^B^ - g r^B^), (5.13) 

s a; = -d;v' - d;"c" = -{d^d - g r^B''/). (5.14) 

Similarly the transformations for the ghost fields and the Nakanishi-Lautrup fields are 
derived: 



sc^ = --gr"'c'd- 


-9r''c'c\ 


se = --gr'c^c\ 


(5.15) 


sc"" = ib'', 




sd = ib\ 


(5.16) 


s6" = 0. 




sb' = 0. 


(5.17) 



The standard Jacobi identity, 

rrst rruv i rrsu rrvt i ersv rrtu n ( ^ A 9i\ 
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splits up into 

rabc cadi , rabd raic , rabi racd r\ /r on\ 

rrab rrcd , rrac rrdb , rrad rrbc r\ (^ '}^\ 

Now we can proceed to fix the gauge. The basic idea is to minimize the norm of the 
off-diagonal part given by 

Rmag = I I dxBl{x)Bl{x). (5.22) 



It is found that a (local) minimum corresponds to 

DfBl = 0. (5.23) 

This is the gauge fixing condition that defines the MAG. In contrast to the Landau gauge 
condition it is non-linear as the diagonal gluon field appears in the covariant derivative. 
This will have far-reaching consequences for the Faddeev-Popov operator and the structure 
of the interactions in the gauge-fixed Lagrangian. 

Having the gauge fixing condition, we can proceed in the usual way of BRST quanti- 
zation, i.e., we calculate the gauge fixing part of the Lagrangian via a BRST variation: 



Smag = s j dxnbfBl-i'^h'^) 



dx {%%" + ib'^DfBl + ebfD^^c" - g r^'e^b^^B^Y- 

- g f^-d^bfB^c'^ - / C/"'7'*5j5^c"c'^ - (1 - C)^ r'^-d'Bld.d), (5.24) 

where the covariant derivative with the hat, D"^, is defined as 

bf:=5''% + gCr'"Al. (5.25) 

It allows to interpolate between the MAG and the Landau gauge via the parameter C^: 
For the MAG it is set to one and for the Landau gauge to zero. The parameter a is the 
conventional gauge fixing parameter for the MAG. Alternative formulations for gauges 
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interpolating between the MAG and the Landau gauge can be found in refs. |173[ 1174) . 

We do not integrate out the Nakanishi-Lautrup field yet, as we have to face the following 
problem of the MAG first: The Lagrangian as we have it now is not renormalizable, as 
there appear divergences due to a quartic ghost interaction |1611I175] . In order to absorb 
these divergences into a counterterni, this interaction is introduced into the Lagrangian 
via a BRST-exact term so that the BRST symmetry is not spoiled: 

Sr = s j dxi-^^Xg^^'c^c^d - -X g P^^ec^c^) = 

= I dx{-i g A r^'Wc^e -i-g\' pH^c^c" + -g'^\p^'r'^'ec^c''c'^+ 

+ X'jg^r^T^'c^c^d^d + X'-g'^r^T^'^^^d^d') (5.26) 

The new terms have two parameters A and A' that can be set to zero to recover the original 
MAG Lagrangian. Note that the second term in the first line is zero in SU{2). Now the 
Nakanishi-Lautrup field is integrated oulQ 

S'macr = Jdx [ebfDl^c^ - g r'^e{D';B;y - g f^'^c-DfBy- 

- / Cr'^^r'^^BlB^c^c" - (1 - C)^ r'^c'-Bld.d - l-{bfB\f+ 

8a a ^ ^ ^^ >"' 

A' AA' 
gr^%bfBf,)c''c'' + — g'^p^'P^'^c^dc^c'^+ 

+ ^g'^Xr^'T'^'e^d'c'^ + A'^^ V"^7'"^'c^c"c'^c* + )^-g^r^T'^%^d'c'^A . (5.27) 

The MAG corresponds to the formal limit a, A, A' — )■ and (^ = L A final simplification 
can be obtained by identifying a, A and A'. This is not only possible because all three 
parameters are zero for the MAG, but it is also supported by a Ward identity, the so- 
called diagonal ghost equation |175] . It only has the meaning of a Ward identity if 
a = A = A'. Identifying these three parameters the terms containing diagonal ghosts drop 
out completely in 5"^/^^ R- Hence the off-diagonal gauge fixing sector does only depend 
on off-diagonal ghosts. The interactions that appear are AAcc, BBcc, ABcc, cccc, Ace 



^Note that in principle one could also keep the Nakanishi-Lautrup field, but for the coming analysis it 
is preferable to set this field on-shell as otherwise there are mixing terms at the two-point level. 
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and Bcc. 

The definition of tlie MAG does only fix tlie gauge for tlie off-diagonal gluon fields 
and the diagonal Yang-Mills part still has a [/(l)^^^-symnietry. This remnant of the full 
gauge symmetry has to be fixed as well for functional equations. From several possible 
gauge fixing prescriptions we choose the Landau gauge, i.e., d^A^ = 0. The corresponding 
gauge fixing parameter is ^ and the required Nakanishi-Lautrup field is again b with a 
diagonal index: 

Sdiag = s dxc\d^A^^ - i-U) = 

= jdx (ib^d.A'^ - i^-U) + c^d,{d,d - g r'^B^c")] . (5.28) 

Integrating out the auxiliary fields also here yields 

SL, = jdx (^Y^{d,Alf + ed,{d,d - g r'-B^c")^ . (5.29) 

In this Lagrangian the diagonal ghosts appear again. However, they can be integrated 
out after shifting the diagonal ghost field as |l5l 1176] 

d^d + ^gr'^B^c". (5.30) 

The Jacobian of this transformation is trivial |176j and so we can get rid of the diagonal 
ghosts entirely. Their disappearance is not surprising as ghosts also disappear in the pure 
Abelian gauge theory in the Landau gauge. 

Now we can combine the individual parts to the final action: 

'S' = Sym + Sj^jj^Qj^ + 5'^jag, (5.31) 

with 



'^MAG,R — '^MAG,R 



X'=X=a 



I dx [ebfDl^d - g f^^ebfBic"- 



g\r''r''B';^B;c''c' - ^{DfBlf + '^g^ r^'^p'^^c'dd^d- 



2 



g r^^JD^^BDc^d + -g^ar^T'^'ec^'c'^ + a-g^ T^T'^'c^cfc'^A , (5.32) 
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Figure 5.1: DSE of the diagonal gluon two-point function. DoDSE cannot draw the usual 
wiggly lines for gluons, so the convention is that continuous lines are diagonal 
gluons, dashed ones off-diagonal gluons and dotted ones ghosts. 



^diag — / dx —{d^A^) . 



(5.33) 



One should note that this Lagrangian simplifies considerably for the gauge group SU{2) 
as the structure constants with three off-diagonal indices are zero: 



Q// 






SU{2) 



dx eofD^^d' - g^ U f BlBied^- 



n ^ n 



M /« 



^{DfBlf + ^-g^ar'^r^^ec'd^c^ 



(5.34) 



SU{2) 



2{d,Bl){d,Bl) - 2{d,Bl){d,Bl) - Ag r'^^BlAUid.B^) - (9.5^)) + 



+ 2/ r'^r'^BlAliBiAi - Aj^Bt). 



(5.35) 



5.2.3 The Dyson-Schwinger equations of the maximally Abelian 
gauge 



The Lagrangian of the MAG, eq. (5.31), contains many interactions of the three fields. 



Therefore the derivation of the DSEs becomes quite lengthy when done by hand. Espe- 
cially the numerous possibilities for different intermediate propagators lead to many terms. 
For this reason the program DoDSE is very useful here. The derivation of the two-point 
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Figure 5.2: DSE of the off-diagonal gluon two-point function. 



function DSEs is explained in Appendix E.2 The resulting equations are depicted in 



figs. 5A_^ 52 and 5.3[ The total number of terms at the two-point level, 16+23+18=57, 



should be compared with that of the Landau gauge, 6+2=8. This illustrates nicely the 
increased complexity of the system of equations that has to be treated here. 

One issue is that diagrams can vanish due to the color algebra. This has to be checked 
explicitly for every diagram. At the two-point level this is possible without assumptions 
for all diagrams except the sunsets, where a dressed four-point function appears. For the 
three-point function there is, under the assumption that no symmetric tensor d^'** appears, 
only one basis tensor in color space, which is the antisymmetric structure function of 
SU{N), /'"'**. The dressed three-point functions are then proportional to this color tensor. 
The two-point functions are all taken as diagonal, i.e., either proportional to 6°"^ or 6^^ . 
Contracting the color tensors in all diagrams except for the sunsets reveals that indeed two 
diagrams vanish: The squint diagrams with a bare BBcc vertex and an internal diagonal 



gluon line. They are depicted in fig. 5.4 Note that other diagrams with a bare BBcc 
vertex remain. 

At this point is has to be stressed that the gauge fixing parameter of the MAG a is 
still arbitrary as it appears in the denominator of the bare vertices. Therefore the ofF- 
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Figure 5.3: DSE of the ofF-diagonal ghost two-point function. 



diagonal gluon propagator has a longitudinal part. In order to take into account that 
the longitudinal part could scale differently as the transverse part one could introduce 
an extra IRE 5B,iong and split the off-diagonal field in transverse and longitudinal parts. 
However, inserting the bare vertices explicitly and projecting the DSE of the off-diagonal 
propagator transversely and longitudinally one finds that the equations for both parts is 
the same. Thus the IR analysis yields 5b = SB,iong- This result can only be invalidated 
by cancelations between different diagrams or due to contributions from the unknown 
dressed vertices. Based on this relation we adopt in the following only one common 
dressing function for both tensors. 

For the moment this is as far as we go in excluding cancelations. It cannot be ruled 
out, though, that nevertheless some diagrams vanish or cancel each other. A clarification 
of this point would require a more detailed analysis taking into account the detailed 
Lorentz structure. However, at the moment it seems not manageable to do so even at 




V 



a) • • b) 

Figure 5.4: Diagrams that vanish due to the color algebra. 
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the IR leading order, as it will be found that the IR dominant diagrams involve dressed 
four-point functions which possess a plethora of possible Lorentz tensors. 

5.3 Infrared scaling solution for the maximally Abelian 
gauge 

Having derived the DSEs of the MAG we can now to proceed to its IR analysis. It will be 
found that there is a unique scaling relation, but the IREs of vertices with an odd number 



of legs are ambiguous. In the discussion of this solution in Section 5.3.3 differences and 
connections to the Landau gauge are highlighted. Finally a numerical calculation of the 
IREs is presented. 

5.3.1 Obtaining the scaling relation 

As the MAG features three fields and eleven bare interaction vertices, its IR analysis seems 
much more complicated than that of the Landau gauge. However, with the improved 
method developed in Chapter |4] only slightly more effort is needed. 
The dressed propagators in the MAG are 

Diip^) = 6-'-^ (,,. - ^) + e^--^, (5.36) 

D%\p^) = 5-^ (,,. - (1 - af-f-^ , (5.37) 

The reason for choosing only one dressing function for the off-diagonal propagator is given 



in Section |5.2.3[ The dressing functions obey power laws in the IR: 

ca{p') " = ' dA ■ {p'Y^ (5.39) 

CB(/)' = °rfB■(/)'^ (5.40) 

c.(p')' = °4-(/)'^ (5.41) 



Now we follow the procedure outlined in Section 4.3.1 The first step is: 
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1. Determine the inequalities that are derived from the interactions appearing in the 



Lagrangian from eq. (4.21) 



In the MAG this leads to the following inequalities in SU{2): 

6a + 5b> 0, 6a + 6,> 0, Sb + S,> 0, 26b > 0, 

25, >0, 5a + ^5b>0, 5a + ^5,>0. (5.42) 

The additional inequalities in SU{N > 2) are 

Isa + ^Sb + S,>0, ^5A + pB>0, ^5b>0, ^5b + 5,>0. (5.43) 

The second step is: 

2. Reduce the number of inequalities, if some of them are contained within others. 
This yields the following four inequalities for general SU{N): 

5a + 5b> 0, 5a + 5c> 0, 5b > 0, 4 > 0. (5.44) 



The inequalities for SU{N > 2), eq. (5.43), do not provide additional constraints for the 
system. Hence we will find that SU{2) and SU{N) have the same IR behavior, although 
their actions are different. 
The last step is: 



3. Try to saturate one inequality after the other according to eq. (4.28) and see, if you 
can find any contradictions with the remaining inequalities. 



We first investigate the last two inequalities in eq. (5.44). They yield Sa = Sb = S, = 0. 
To illustrate this I explain the line of argument for the case 6b = in more detail. 

The equation 5^ = is obtained from diagrams with a bare B self coupling. As such 
bare vertices do not appear in the DSEs for the diagonal gluon or the ghost two-point 
functions, we still have to determine their IR leading diagrams. In the DSE of the diagonal 
gluon there are two possibilities for IR leading diagrams, those with bare AABB or AAcc 
vertices. If the latter are leading we obtain 6a + S, = 0. However, since both 6a and 
6c are non-negative if 6b = 0, the only possible solution is 6a = 6, = 0. Similarly we 
get (5a = 0, if vertices with a bare AABB vertex are leading. For 6c the ghost two-point 
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Figure 5.5: Possibly IR leading diagrams in the DSE of the diagonal gluon two-point 
function. The sunset diagrams (first row) are definitely at leading order, while 
it is possible that the squint diagrams (second row) are IR subdominant. 



function is analyzed, where we come back to the argument with the AAcc vertex. Hence 
the emerging solution is that all three propagators scale trivially in the IR, i.e., equally 
as in the UV. This option is of no interest to us, as we do not expect that perturbative 
propagators describe the IR regime of Yang-Mills theory. The argument if starting with 
the equation 5c = is completely the same. 

The two remaining possibilities are 6a + Sb = and Sa + Sc = 0, which can be fulfilled 
independently. Thus the scaling relation of the MAG is 



kmag '■= -Sa = 5b = 5c>0. 



(5.45) 



Before discussing its meaning we determine the IR behavior of vertices. 
5.3.2 The infrared exponents of vertices 

From the scaling relation found above the IR behavior of vertices can be obtained. A first 
observation is that the IREs of AABB and AAcc vertices vanish: 



I^AABB — I^AAc 



0. 



(5.46) 



This can be shown by plugging the scaling relation eq. (5.45) into eq. (4.16) and using 



that primitively divergent vertices have a non-positive IRE, see eq. (4.20) 



Next we want to determine the IR leading diagrams so that we can derive the IREs of 
other vertices from them. However, for the MAG there is an ambiguity due to the fact 
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Figure 5.6: Possibly IR leading diagrams in the DSEs of the off-diagonal gluon two-point 
function (upper row) and of the ghost two-point function (lower row). The 
sunset diagrams (left) are definitely at leading order, while it is possible that 
the squint diagrams (right) are IR subdominant. 



that there are two classes of possibly leading diagrams. The origin of the scaling relation 
are diagrams with bare AABB or AAcc vertices, i.e., (neglecting the tadpole) sunset and 
squint diagrams can be used to derive the scaling relation. The corresponding diagrams 



are depicted in figs. |5.5| and |5.6[ As we did not yet specify in our analysis which diagram 
led to the scaling relation, we look for the consequences if only one class is IR leading. 
Let us assume first it is the squint diagram from which we derived the scaling relation 



and eq. (5.46 ). This provides enough information to count the IREs of the sunset diagrams 
which turn out to be at leading order, i.e., if the squints are IR leading also the sunsets 
are. If, however, we assume that the scaling relation is obtained from the sunset diagrams, 
we cannot determine the IREs of the three-point vertices necessary to count the IREs of 
the squint diagrams. So the sunset diagrams are always at leading order, while we cannot 
say if the squint diagrams scale equally. 

The reason why we cannot determine the IREs of three-point functions is purely com- 
binatorial. We found that the IR leading vertices have four legs. Consequently it will 
be easy to construct the IR leading diagrams of an n-point function if n is even: As de- 



picted in fig. 5.7 we plug AABB or AAcc vertices into appropriate IR leading diagrams of 



two-point functions. The IREs of the resulting diagrams can be inferred from eq. (4.24) 



but as we only add quantities for which the last two terms are zero, we directly get to 



eq. (4.25), the maximally IR divergent solution. The IRE of a vertex with n^ diagonal 
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Figure 5.7: Consecutively inserting pairs of A, B and c fields via AABB and AAcc 
vertices shows that for graphs with an even number of legs the maximally IR 
divergent solution is reahzed. 



gluon legs, Ub off-diagonal gluon legs and Uc ghost legs is then 

K{nA,nB,nc) = -{ua - ub - nc)KMAG {nA + ub + ric even). (5.47) 

However, it is not possible to construct an IR leading diagram purely from four-point 
functions for a vertex with an odd number of legs. The open question is the IR behavior 
of the three-point function for which we can only find upper and lower bounds from the 
bare vertices and the maximally IR divergent solution, respectively: 

1 1 

--;^HMAG < I^AB2 < 0, - -KmAG < f^Ac2 < 0, 

3 3 

-^f^MAG < f^BBB < 0, - -KmAG < f^Bcc < 0. (5.48) 

If the IRE of the three-point functions was known, the IREs of the complete tower of 
vertices with an odd number of legs could be derived. Without this knowledge, which 
cannot be obtained from a power counting analysis, we have to introduce a parameter t] 
to take the different possibilities into account. The IRE of an n-point function with n odd 
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reads then 



n{nA,nB,nc) = -{nA-nB-nc + r])KMAG {nA + ub + ric odd). (5.49) 

For SU{2) 1] can be either or 1, while for SU{N > 2) it can also be 2 and 3. 



5.3.3 Aspects of the infrared solution and comparison to Landau 
gauge 

The IR solution for the MAG can be summarized as follows. The scaling relation reads 



f^MAG '■= -Sa = 5b = 5c>0 



(5.50) 



and the IREs of vertices are given by 



K{nA, riB, He) = -{riA - riB - nc)KMAG {ua + nB + nc even) 



(5.51) 



and 



K{nA,nB,nc) = -{ua-Ub -nc + r])KMAG (nA + nB + ric odd). 



(5.52) 



One can see directly from the scaling relation that this solution is qualitatively quite 
different from the one in Landau gauge, where the ghosts are IR enhanced, as here it is the 
diagonal gluon propagator that is the dominant degree of freedom in the IR. This can be 
interpreted as a variant of Abelian IR dominance. Also in agreement with this hypothesis 
is the IR suppression of off-diagonal propagators. One should note, however, that it is 
different from the indication of Abelian IR dominance seen in lattice calculations of the 
MAG propagators where all propagators show a massive behavior. At momenta below the 
momentum scale induced by this mass the propagators decouple. Since the propagator of 
the diagonal gluon is least IR suppressed the diagonal degrees of freedom are dominant in 
the IR. Here, however, we find an IR enhanced diagonal gluon propagator and it is clear 
that in the DSEs/FRGEs the diagrams with the most diagonal propagators dominate. 

There is even a direct connection between ghost dominance in the Landau gauge and 
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Abelian dominance in the MAG. One can show that upon transforming gauge field con- 
figurations from the MAG to the Landau gauge that the diagonal (MAG- dominant) con- 
figurations are transformed to configurations lying on the Landau gauge Gribov horizon 
|177) . which give rise to the IR enhancement of the ghost propagator. This establishes 
not only a link between the dominant configurations of the Landau gauge and the MAG, 
but also between two confinement scenarios of different origin |177) . 

Besides the obvious difference that the dominant fields are different, there is also another 
aspect that distinguishes the MAG and the Landau gauge. The scaling relation in the 
latter is < 2k,lg := 6a = — 2(5c. A lower bound for the IREs of propagators is the value 
when the Fourier transforms of the propagators are no longer well-defined [55| ^U\, i.e., 
here 6c > —1 and 6a > —1 which corresponds to propagators with a 1/p^ behavior. Hence 
we get 

0<hlg< 1- (5.53) 

Note, however, that more divergent propagators, specifically going like 1/p^, have been 
found also and interpreted via an analytic continuation of the exponent beyond —2 [56]. 
These bounds on Klg allow that the gluon propagator vanishes at zero momentum, for 
which we need 6a > 1, i-e-, k,lg > 1/2- In the MAG the upper bound is kmag < 1- The 
condition for an IR vanishing propagator, on the other hand, leads to kmag > 1- Thus 
the IR suppressed propagators in the MAG cannot vanish. On the contrary, they are still 
IR divergent, but less than the tree-level propagator. It should be stressed that also such 
propagators can violate positivity and belong to confined fields. 

The reason for different IR dominant field configurations in the MAG and the Landau 
gauge is likely related to the different shapes of their Gribov regions. While both are 
convex and contain the origin, only the latter is bounded in all directions. The former, 
however, is unbounded in the direction of diagonal fields |i^l45] . In this case the entropy 
argument that in a high-dimensional space most volume is contained at the boundary 
and thus configurations with a small value of the lowest eigenvalue of the Faddeev- Popov 
determinant are dominant does not apply |27| . The unboundedness of the Abelian direc- 
tion in field configuration space may be directly responsible for the enhancement of the 
diagonal propagator and the dominance of Abelian configurations. 

The discrepancy presently observed between the results from this thesis and from lattice 
calculations |144t I172J may have its explanation in the existence of an additional solution 
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for the MAG of the decouphng type which is also found in the refined Gribov-Zwanziger 



scenario |I5| 1145) . Indeed one can see from the DSEs that a decouphng solution in 
the MAG features three massive propagators. The reason is that there are four-point 
interactions between all three fields {AABB, AAcc, BBcc) which lead to tadpole diagrams 
in the DSEs. If the internal propagator of a tadpole diagram is massive, it does not vanish 
in dimensional regularization and can render the solution of the DSE massive. In Landau 
gauge there exists only a gluon four-point interaction and thus a massive gluon propagator 
does not influence the ghost propagator via a tadpole diagram. In the MAG, however, a 
massive propagator automatically leads to the massiveness of the other two propagators 
by the tadpoles stemming from the vertices AABB, AAcc and BBcc. 

The recent suggestion that the decoupling and the scaling solutions in the Landau 
gauge correspond to different Landau gauges that only deviate from each other in the IR 
regime |20] can also be considered in the MAG. Indeed its DSEs offer a similar setting. 
In the Landau gauge the realization of one of the two solutions is directly related to 
the ghost propagator at zero momentum. Only if the boundary condition is chosen such 
in the renormalization process that the ghost dressing function diverges in the IR, the 
scaling solution is obtained. In the MAG the same condition is required for the diagonal 
gluon propagator. Thus the zero momentum value of its dressing function may serve as a 
parameter in analogy to the B-parameter of the Landau-B gauges |20] . 



5.3.4 Calculation of the parameter k,mag 

The first step in establishing the existence of a scaling solution in the MAG was to find 



the scaling relation eq. (5.45). The next step would be a numerical solution of the DSEs 
over the complete momentum range. However, this is complicated mainly by the fact that 
the IR leading terms are two-loop diagrams. For a numerical solution the DSEs have to 
be truncated properly. This can be done in a consistent way if the IR leading diagrams 
only have one loop, as a one-loop truncation then contains the IR and UV leading parts. 
For two loops, however, a direct UV/IR consistent truncation does not exist: As some 
two-loop terms are required for the IR, one has to include all diagrams up to two loops in 
order to be consistent in the UV. Such a calculation would not only be of a tremendous 
complexity due to the high number of terms in the MAG, but is also aggravated by the 
fact that a new method for renormalizing and calculating the two-loop diagrams and 
ansatze for all involved vertices are required. 
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What can be done is the calculation of the IRE Kmag if we only take into account the 
sunset diagrams. It is not possible to include the squint diagrams as in this case we do not 
have enough equations for the determination of the free parameters. The free parameters 
are the coefficients in the power laws of the dressing functions in the IR, (Ia, ds, dc, and 
the IRE Km AG- So in total we have four unknowns and three equations. However, in the 
DSEs truncated to the sunsets the coefficients appear only in the combinations d^d"^^ and 
d^d"^ and we can solve this system. If we included the squints not only would we need 
to know the IR behavior of the three-point functions, but also other combinations of the 
coefficients would appear rendering the system of equations unsolvable. Thus we assume 
here that only the sunset diagrams are IR leading. 

For the calculation it suffices to determine the integrands in the IR and the system of 
DSEs reduces to the following simpler form after suitable projections in Lorentz space: 

dA^ = -X^ABBi.P^^ I^MAG)dAd\ - XJ^^^{p'^, K,MAG)dAdl, (5.54) 

d^ = —Xj^^^^{p ,KMAG)dAdB, (5.55) 

C^ = -XAAccip"^, nMAG)d\dc. (5.56) 

The X(p^, Kmag) denote the sunset integrals without the coefficients from the propagator 
power laws. The superscript gives the corresponding DSE and the subscript the bare 
vertex contained in the diagram. Using the invariant combinations 

h :=44, (5.57) 

h ■■=d\dl (5.58) 

the three equations can be combined to 



-^ ^ XaABbJP iI^MAg) X^j^^^ip , KmAg) /g ggx 

^AABBiP'^^f^MAG) X^^^^ip^ , KmAg)' 

If the X(p^, Kmag) are known, this equation yields the solution(s) for Kmag- 

For the calculation of the X(p^, kmag) "we need to decide what we use for the dressed 
four-point functions. As we cannot solve their DSEs and get exact expressions we can only 
make ansatze. The simplest choice is using bare vertices, but we could use any expression 
that stays constant when all external momenta go to zero. However, the number of 
possibilities for modifying the dressed four-point functions is tremendous. First of all 
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four-point functions have a large number of possible tensors: For two Lorentz indices 
{AAcc vertex) there are ten possible tensors and for four Lorentz indices there are 138 
{AABB vertex). Of these five and 76, respectively, are transverse with respect to the 
diagonal gluon legs. In color space the situation is a little bit better, as the fact that two 



indices are diagonal allows only a subset of all possible color tensors (see Section 6.2.1 for 
details). One could choose, for example, the following basis for SU{N): 

where d^*"' is the totally antisymmetric tensor 

df^ = -Tr (T^T^^TT^^) . (5.61) 

The Jacobi identity tells us that another possible tensor, y«icjfeic^ jg j^gj^^gfj ^q f^ ^ g^j^fj 
hence does not appear in this basis. For the calculation of Kmag we restrict ourselves to 
bare vertices. 

Another thing we need is an analytic solution for the scalar integral 

Iss{a,b,c,eJ;p'):= I 0^^^/q'nrm{p-q-rm2pqn2pry (5.62) 

for {e, /} G N. As I am not aware of any existing solution in the literature except for 



e = I provide it in Appendix D 

The integrals Xlp"^, kmag) can be calculated with Mathematica. First the DSEs are 



derived with the package DoDSE as outlined in Appendix E.2 Using the Feynman rules 



given in Appendix O the detailed expressions of the integrals are derived. They are then 



processed by contracting Lorentz and color indices and calculated using eq. (D.4). The 
contraction of color indices was done in Mathematica by simple replacement rules, as no 
available program could deal directly with the splitting of the color algebra. The whole 
process is necessarily automated, because the resulting expressions cover several pages 
even in the simple case adopted here. This also allows to test the effect of different 
tensors or employing different dressings. 

Having all analytic expressions for the sunset integrals X{p'^.,HiMAG)-, we can plot the 



right-hand side of eq. (5.59) as a function of kmag- There is still one unfixed value which 



is the gauge fixing parameter a. The obvious choice is a = 0, which corresponds to the 
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Figure 5.8: Left: Both sides of eq. (5.59) for several values of a. Right: Dependence of 



f^MAG on a. The dots and triangles represent the first and second solution 
branches, respectively. 



pure MAG, but it is also interesting to check the dependence on a. In fig. 5.8 on the left 



the two sides of eq. (5.59) are plotted for several values of a. Solutions for Kmag can be 



read off where the straight line (left-hand side of eq. (5.59)) crosses the curves (right-hand 
side of eq. (5.59)). The singularity at kmag = 0.739908 seems to be universal for all 
values of a; this was checked up to a = 100. Consequently there is always one solution 
i^MAG ~ 0.74. However, as is depicted in eq. (5.59) on the right, there exists a second 



branch of solutions. For a between approximately 0.5 and 0.9 both coincide and when at 
a ~ 2.2 the first branch starts to go away from kmag ~ 0.74, the first branch takes over. 
The first branch crosses Kmag = 1 at a = 3.359 and becomes thus unphysical. 

As there is a solution for the parameter Kmag we have one more indication that the ob- 
tained IR solution really exists. The final proof would be a complete numerical calculation 
for all momenta. 
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Chapter 6 

Dyson-Schwinger equations in the 
Gribov region 



In Chapter |2] I explained how to fix the gauge according to the conventional Faddeev- 
Popov procedure and that for non-perturbative calculations this is not sufficient due to 
the appearance of further gauge copies. Now I will go into the details of the improved 
gauge fixing suggested by Gribov and how this leads to a local Lagrangian amenable to the 
usual tools of quantum field theory. The resulting action is called Gribov-Zwanziger action 
[26| [27] and constitutes the best option available so far to achieve a complete gauge fixing 
within the path integral. For lattice simulations the situation is slightly different and I 



will comment on this in more detail in Section 6.1.3 Of course the ultimate goal would 
be a complete gauge ffxing as, for example, provided by restriction to the fundamental 
modular region. This, however, does not seem realizable directly within the path integral 
due to the non-trivial topology of this region. 

One should note that the underlying idea of Gribov |26| works for all gauges which are 
deffned by a minimization problem like Landau gauge. Coulomb gauge or the maximally 
Abelian gauge. The ffrst two were already considered in Gribov's original paper |26| . 
while the last one was mainly investigated by Sorella and collaborators, see, e.g., refs. 
|15| 11451 I176i 1178] . In this chapter, however, I only consider the Landau gauge. 



The Gribov-Zwanziger action is described in Section 6.1 and will be investigated in 



Section 6.2 with the tools developed in Chapter |4] In this action ffelds mix at the level 



of two-point functions what leads to several complications that require a refinement of 



the method as described in Section 4.5 When the smoke clears two possible IR solutions 



remain, which, however, are very similar up to details. In particular, both feature the same 
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qualitative behavior as the result obtained from the Faddeev-Popov theory. Interestingly 
one solution reduces in the IR exactly to the Faddeev-Popov theory. The results of this 
chapter have been published in ref. |179] . 

6.1 The Gribov-Zwanziger action 

This section contains a short overview of the derivation of the local formulation of the 
Gribov-Zwanziger action. It is named after N. Gribov, who was the first to suggest a 
solution to overcome the incomplete gauge fixing, and D. Zwanziger, who subsequently 
brought this idea into its final form in terms of a local action and derived several properties 
of the first Gribov region. 

6.1.1 Restriction to the first Gribov region 

As mentioned in Section [2] the first Gribov region Q is defined as the set of gauge field 
configurations for which the Faddeev-Popov operator is positive and which fulfill the gauge 
fixing condition, i.e., for Landau gauge 

n:={A\d^Af, = 0, M >0} (6.1) 

with M = —d^D^ the Faddeev-Popov operator. The boundary of this region is called the 
(first) Gribov horizon. It is defined by those gauge field configurations for which the lowest 
non-trivial eigenvalue of the Faddeev-Popov operator vanishes. Consequently the second 
Gribov horizon is where the second eigenvalue vanishes and so on. The conventional 
gauge fixing already restricts the field configurations to fields obeying the Landau gauge 
condition d^A^^ = and only the second condition, the positivity of the Faddeev-Popov 
operator, has to be implemented additionally. 

The Faddeev-Popov operator is related to the ghost propagator D'^{k) by 

Dt{k) = -^c-Uk) = -{M-T- (6.2) 

Hence in order to be within the first Gribov region the ghost propagator dressing function 
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Ccc{k) must be positive. Gribov parametrized the ghost propagator by [26] 

As a{k, A) increases with k decreasing, it is sufficient to demand that o"(0, A) < 1. This is 
known as the no-pole condition |26]: When the ghost form factor cr(0, A) becomes one, the 
Faddeev-Popov operator has a zero eigenvalue and one is directly on a horizon. o"(0. A) 
is calculated as a series in A, which can be regarded as an external gluon field here. The 
no-pole condition is enforced in the path integral via a Heaviside functional: 

0(l-a(O,A))= / -^ ' . . (6.4) 

zvr z p — le 



oo 



The integral over /3 can be evaluated with the saddle-point method. It yields a value for 
/3 — )■ 7^ which takes the role of a new mass parameter. 7 is called the Gribov parameter. 
The positivity of the Faddeev-Popov operator means that all its eigenvalues A„ are 
positive. Relaxing this condition to tr(A„) > allowed Zwanziger to derive a closed 
expression, the so-called horizon function, that restricts the integration to the Gribov 
region [27] : 

h{x)= hm f dy{Dl\x)^\x)){M~y\x,y){Dl\y)^\y)). (6.5) 

im-^i J 

The limit 7(3;) — )■ 7 can only be taken after localizing the action. Details on this can 
be found in ref. [180] . The horizon function corresponds to a generalization of the ghost 
form factor cr(0. A) to all orders. This has explicitly been worked out up to third order 

[m]. 

The condition for being within the first Gribov region is phrased as the horizon condi- 
tion: 

dxh{x)<d^^{N'^ -l)V, (6.6) 

where d, N and V are the number of dimensions, the number of colors and the space-time 
volume, respectively. Again this condition can be enforced via a Heaviside functional. 
Details on how to get from the step function to an additional term in the action can 
be found in ref. [56]. An interesting observation is that the resulting path integral has 
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similarities with the partition function of a canonical ensemble which is equivalent to the 
microcanonical ensemble in the thermodynamic limit. Thus the Heaviside functional can 
be replaced by a Dirac delta functional here, i.e., the horizon condition becomes 

dxh{x) = d-f\N^ ~1)V. (6.7) 

This amounts to restricting the integration to configurations lying directly at the Gribov 
horizon. 

The resulting action that enforces the restriction to the first Gribov region is known as 
Gribov-Zwanziger action given by 

inon— local 



gnon-local ^ g^^ ^ g^^ ^g^^ 



W 



here 



Sfp — Sym + Sgf, (6.9) 

Sh= f dx{h{x)--fU{N^-l)), (6.10) 



SyM = ^ I dxF;,F;,, (6.11) 



Sgf = dx {i 6"(<9^A;) + c" M'' c') . (6.12) 

The trivial term '~f'^d{N'^ — 1) is introduced for later convenience. Since Sh is non-local 
the standard tools of quantum field theory can not be employed. However, it is possible 
to localize it as discussed in the next section. 

6.1.2 The local action 

In order to localize the horizon function we need two pairs of additional fields: ip°^ and 
1^"^, which are complex conjugate to each other, and uj"f and u}"f, which are also complex 
conjugate to each other but Grassmann fields. These new fields are BRST doublets: 

svf = co;\ s< = 0, (6.13) 

suf = ^f, s^f = 0. (6.14) 
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The three indices are one Lorentz and two color indices in the adjoint representation, 
i.e., /i = 1, . . . , c? and a,b = 1, . . . , A^^ — 1. The bosonic fields are used for localizing the 
non-local term Sh and the fermionic fields for canceling the determinant arising from that 
procedure. Using the Gaussian integration formulas 

(6.15) 
V[uu]e^''^''^^'''< =(det M)'^(^'-^) (6.16) 



one can show 

The local Gribov-Zwanziger action is then given by 

^'°-^ =5^p + 5gz, (6.18) 

Sgz =jdx «M"Vt^ + 7' 9 r'-A^iy,'; - ^'p - u;^Mco'; - dj' {N' - 1)) . (6.19) 

Note that a slightly different version is sometimes used that contains an additional term 

Sa,,st = jdxg r'^ul^d, [{DfC)^'^] . (6.20) 

It can be obtained by a shift of the u field, 

<(x) -> ujfix) -gjdy {M-'rix, y)r''^d, [{Df c^ {y))^^^{y)\ , (6.21) 

and has the advantage that the remaining BRST symmetry is more manifest, since almost 



all terms in eq. (6.19) are BRST exact: 



.(cD-iVf^Vi^) = ^r^^Vi' - <M»^a;J^ + g r'^ul^d, [{DtV)^';] . (6.22) 

This property is useful in the process of renormalization p3]. The additional vertex 
derived from Sabjist does not infiuence any calculations as it features the fields c and u, 
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but there is no corresponding term with c and u. Thus this vertex does not appear in 
any diagrams. 

The only term in the action 5''°'^"' + S/^g^^^r^ which is not BRST invariant is the one 
proportional to 7^: 

s (7^ g r'^A^iv"; - ^'p) = 7' 9 r'' {{-Dfcy^'; - ^'p + Ay;) . (6.23) 

This is a soft breaking, i.e., only of mass dimension twojj Hence one can show the 
renormalizability of the Gribov-Zwanziger action by embedding it into a larger BRST 
invariant action, for which the usual tools of algebraic renormalization apply [51]. The 
complete proof of renormalizability can be found in ref. |18U] . The original literature is 
refs. |33l [T82l[T83] . 

Another consequence of the BRST breaking is that physical operators can no longer be 
determined by their cohomology, because BRST exact quantities, with which a physical 
operator mixes in the renormalization process, can now influence expectation values. The 
complete form of the operator F^^ under this mixing was determined in ref. |184] for the 
Gribov-Zwanziger Lagrangian. Another approach was taken in ref. |185] . where operators 
with a physical cut structure were constructed by adding BRST exact terms. The precise 
meaning of the BRST breaking, however, remains to be clarified further. 

In terms of the new fields the horizon condition becomes 

{g r'^'A'^.i^'; - ^p) = 2d^' {N' - 1) (6.24) 



as can be inferred by differentiation of eq. (6.17) with respect to 7^. Often this is phrased 
in terms of the gap equation 

where r[0] is the vacuum energy defined by 

e-r[o] ^ f V[Acapipuuj]e-'"'" . (6.26) 

Writing the horizon condition in this form a spurious solution for the value of the Gribov 



""^The ghost fields c and c have mass dimensions zero and two, respectively 
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parameter arises, namely 7 = 0. As this corresponds to the starting action, it has to be 
discarded |97] . 

One point which has been stressed repeatedly in the literature is that the Gribov- 
Zwanziger action is tightly connected to the gap equation as was already pointed out 
by Gribov [26j. Only if 7 takes the value as determined by the horizon condition the 
Gribov-Zwanziger Lagrangian corresponds to a gauge theory. It was pointed out in ref. 
that the gap equation breaks several symmetries of the original effective action F. It 



is argued that this symmetry breaking is spontaneous and the related Goldstone particles 
are the fields c, c, u, u and parts of (f and ip. Their propagators are all found to be IR 
enhanced like l//c^. 

It is advantageous for the present task to rewrite the local Gribov-Zwanziger action 



given in eq. (6.18) by splitting the bosonic auxiliary fields into real and imaginary parts 



U and V, respectively, as done in ref. |56 



qlocal 


— Sfp + Sqz^ 


ql 


= Su + Sv + Suv - ^Muj 


Su 


= ^- fdxu;^M^''u';, 



(6.27) 
(6.28) 

(6.29) 



Sv = ^ j dx v;" ivf^^ v;^^ + 1 g -f^V2r'" f dx a;^v;^^ (6.30) 



Suv = 2^ an j dx U^Xd-^l^ (6-31) 

where U and V are defined by 

^ = ^{U + iV), ^=l^[U-iV). (6.32) 

When the Landau gauge condition d^A^ = is enforced, Cuv vanishes and the only 
mixing on the level of two-point functions is between the gluon field A and the imaginary 
part of the bosonic auxiliary field V , whereas the U field does not mix. This splitting 
simplifies calculations, because we only have to deal with a two-by-two matrix instead 
a three-by-three matrix for the mixing. Perturbative calculations for the propagators 
derived from this Lagrangian can be found in ref. [187j, where the U and V fields are 
denoted by p and ,^, respectively. 

A further simplification is achieved by combining the Faddeev-Popov ghosts c and c. 
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the fermionic auxiliary fields u and u and the real part of the bosonic auxiliary field 
U into one single field. This is possible, since all of them only interact with the gluon 
field via the Faddeev-Popov operator and are quadratic in the action. Hence they can be 
integrated out in the path integral: 



V[cc]e^^^ = detM, (6.33) 

I^Me^^'" = (detM)''(^'-^\ (6.34) 

/"p[t/]e-^^^^ = (rfetM)-i(^'-i). ^g 35) 

The different exponents of the determinant of the Faddeev-Popov operator are due to the 
different numbers of degrees of freedom. For 7 = also the V field can be integrated 
out and all determinants from auxiliary fields cancel so that the original Faddeev-Popov 
Lagrangian is recovered as required. For the present purpose we can treat all these non- 
mixing fields as new fermionic fields rj and fj with the appropriate number of degrees of 
freedom by localizing the resulting determinant again. The field V cannot be included 
due to its mixing with the gluon field and therefore the two fields V and rj can have a 



different infrared behavior. An overview of the different fields is given in tab. 6.1 The 
final action reads 



/^^U^;.^;.+^(^M^/^ 



+ Iv;- M^' V^^ + I g ^'V2r'"^AX' + \^ gr'^U^X'^^^l) ' (6-36) 

where the subscript index of the new ghost fields 77 and f] runs from 1 to ^(A^^ — 1) + 1 and 
the superscript index is the usual color index, running from 1 to A^^ — 1. Note that for odd 
dimensions and even A^ this number is half-integer and therefore this transformation is not 
directly possible. However, we can consider only integer values and perform an analytic 
continuation to half-integer values, if necessary. Alternatively one keeps the original fields 
separated and will get the appropriate numerical factors in front of diagrams. 

A further diagonalization of the Lagrangian would require a further splitting of the 
fields in color space due to the different number of color indices of the A and V fields. 
The new fields would mean a significant complication of the Lagrangian at the level of 
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Field 


Numb 


er of degrees of freedom 


Statistics 


c, c 

4>, 4> 


1 
d{N^ - 1) 
diN^ - 1) 


fermionic 

fermionic 

bosonic 


u 

V 




d/2{N^ - 1) 

d/2{N^ - 1) 

d/2{N^ - 1) + 1 


bosonic 

bosonic 

fermionic 



Table 6.1: The numbers of degrees of freedom and the statistics of the Faddeev-Popov 
ghosts c and c, the original auxiliary fields u, Co, (p and (p, the bosonic auxiliary 
fields U and V and the fermionic fields rj and f/. A factor N"^ — 1 from the 
adjoint index common to all fields is not taken into account here. 



vertices although at the two-point level it becomes simpler. Such a diagonalization was 
performed in ref. |il85j and in a slightly simpler form in ref. [183| and leads to fields 
with complex masses. Due to the resulting structure of the vertices of the diagonalized 



Lagrangian we continue with the expression given in eq. (6.36) 



6.1.3 The first Gribov region on the lattice 

The restriction to the first Gribov region works differently on the lattice. In fact, the 
usual gauge fixing algorithms automatically lead to the first Gribov region as they are 
based on the minimization of a functional corresponding to eq. (2.28). Consequently 



it is not possible to fix with such a method to another Gribov region as the first one. 
Configurations just beyond the Gribov horizon can be obtained due to limitations in the 
numerical accuracy as was demonstrated in ref. |189) . where in addition also an alternative 
approach for gauge fixing, stochastic quantization, was employed. 

Since there are several minima along a gauge orbit, it depends on the algorithm and 
its parameters which minima are chosen. This leads to the definition of several Landau 
gauges. The minimal Landau gauge takes the first minimum found by the minimization 
algorithm, see, e.g., refs. [3211190) . The absolute Landau gauge is defined by taking those 



gauge field configurations that correspond to the global minimum of eq. (2.28) and thus 



to the fundamental modular region. This is a hard computational problem similar to 
spin-glass problems and can only be realized approximately. Calculations where this was 
attempted can be found, for example, in refs. |^ IHUt 1191^ 1192] . 

Of course this is not the only possible prescription to choose a unique representative of 



93 



6.1 The Gribov-Zwanziger action 



each gauge orbit as in principle every choice is allowed. A class of recently suggested gauges 
are the Landau-B gauges [20j. Instead of using the value of the minimizing functional, 
which is equivalent to the trace of the gluon propagator, as a criterion to decide which 
gauge copy is used, one takes the value of the ghost propagator at the lowest available 
lattice momentum. This value is denoted by B. Using B instead of the trace of the 
gluon propagator is in the spirit of functional equations, where the value of the ghost 
dressing function at vanishing momentum serves as a boundary condition for the system 
of equations. Two examples of Landau-B gauges are the niin-B and the max-B gauges, 
where the lowest and highest values of B are chosen, respectively. It is possible that the 
max-B gauge corresponds to the scaling solution |2U) . but currently only data for small 
lattices in three dimensions [20] and in the strong coupling limit |193] are available. 

6.1.4 The Dyson-Schwinger equations of the Gribov-Zwanziger 
action 

A manual derivation of the Dyson-Schwinger equations becomes quite tedious because the 
mixed propagator leads to many additional terms. For example, already at the perturba- 
tive level an AAV vertex appears at one loop. For the rj field things stay simple as their 
number is conserved. This is due to invariance under the same scale transformation that 
guarantees ghost number conservation in the Faddeev- Popov theory [62] : 

r] -^rje^, (6.37) 

f] ->r/e"^ (6.38) 

Thus in a sense the r] fields behave similar as the Faddeev-Popov ghosts. However, one 
should not forget that they couple with the V field, e.g., via the (dressed) Vfjri vertex. 

The non-trivial relation between propagators and two-point functions in the Gribov- 
Zwanziger theory is not an obstacle for the derivation of the DSEs, as dressed two-point 
functions only appear on the left-hand sides of their own DSEs and the bare counterparts 
on the right-hand sides. It is again advantageous to employ the Mathematica package 
DoDSE which can also handle mixed propagators. The resulting DSEs are given dia- 



grammatically in figs. |6.1[ |6.2[ |6.3| and |6.4[ Their derivation is described in Appendix 
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Figure 6.1: The DSE of the gluon two-point function. A propagator with the exponent 
— 1 denotes the two-point function. This convention, strictly speaking being 
mathematically incorrect, is chosen for the purpose of diagrammatic represen- 
tation only. The propagators and vertices are labeled by the respective fields. 
The indices i and j denote the first and second fields of the depicted two-point 
function DSE. Internal propagators are all dressed. Bare and dressed n-point 
functions are denoted by thin and thick blobs, respectively. 
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Figure 6.2: The DSE of the V field two-point function. 
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Figure 6.3: The DSEs of the A-V"-mixed two-point function. The first is the AV DSE and 
the second one the VA DSE. 
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( _ > '' • " ■ ''— • _ 



Figure 6.4: The DSE of the rj field two-point function. 

6.2 Infrared analysis of the Gribov-Zwanziger action 

As already mentioned the matrix relation between propagators and two-point functions 
complicates the analysis as it yields non-trivial relations between the two sets of quantities. 
In this section the details of this relation are worked out. It is found that four possible 



scenarios for the IR behavior emerge which are discussed in Section |6.2.2[ Two of these 
scenarios can be discarded as they turn out to be inconsistent and the remaining two lead 
to the same qualitative IR behavior of the propagators as in the Faddeev-Popov theory. 

6.2.1 Two-point functions and propagators 

In addition to the complication due to the mixing at the two-point level the IR analysis 
is complicated by the tensor structure of the propagators/two-point functions that has 
to be considered explicitly as different dressing functions can have a different momentum 
dependence. The combination of mixing and several dressing functions leads to an even 
richer structure, as the matrix inversion counterintuitively can lead to terms where the 
determinant does not appear. 

I start by discussing the color structure. The gluon correlation function has only two 
color indices. Thus the only color tensor is 6"'''. I exclude here explicitly the occurrence 
of an antisymmetric tensor e""^ as there is no indication of such a tensor neither in lattice 
simulations nor in other calculations. The A-V mixed correlation function has three color 
indices. Excluding here the totally symmetric tensor d"'^'^, which cannot appear directly 
from any diagram, only the totally antisymmetric tensor /"^'^ remains. The VV correlation 
function poses the greatest challenge as it has four color indices. This allows in the most 
general case six different color tensors |194l 1195] : 

^ab^cd^ ^ac^bd^ ^ad^bc^ jabe jcde ^ jace jbde ^ ^abcd _ ^g gg^. 
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The last tensor is the totally antisymmetric tensor of rank four given by 

where T" is the generator of the group SU{N). Other tensors can be constructed, but 
they are not linearly independent; e.g., y^^^ey^ce jg ^gig^^gj ^q j-abej-cde ^^^ j-acejbde ^^^ ^j^g 

Jacobi identity. Similar identities exist that allow to express tensors constructed from the 
totally symmetric tensor d"''"^ with the six basis tensors given above. 

The structure in Lorentz space is simpler: As both the A and the V field have one 
Lorentz index, all three correlation functions can have a transverse and a longitudinal 
tensor. 

Finally there are the 77 and f] fields. They have one adjoint index and one index that 
runs from 1 to |(A^^ — 1) + 1. The behavior of the latter index is best explained in terms 
of the original fields c", c", w"^, w"^ and f/"^. The second color index and the Lorentz 
index of the last three fields were only introduced to obtain the correct number of degrees 
of freedom. They are in some sense static: In diagrams with external u/u/U fields and 
no 00/00/U loops they do not contribute and in u/u/U loops they only lead to additional 
numerical coefficients. Thus they can be combined to the second index of the 7] and fj fields 
together with a trivial 1 from the c and c fields. In the calculations this was most easily 
implemented by using the Faddeev-Popov ghost instead of the 7] field and supplementing 
the results afterwards with the correct coefficients. 

In the present work not the full color tensor basis was used. Instead only the tensors 
appearing in the action were taken into account for the two-point functions. Nevertheless 
this truncation results in a non-trivial propagator for the V field. The matrix of dressed 
two-point functions is defined by 

(yAA yAV ' 

where the individual two-point functions are given by 

^AA,ac ^ ^acp2^±^p2^p^^ ^ S^^^c\{p')p,p,, (6.42) 

^vv,abcd ^ s^^6'VMp')g,., (6.43) 

ri::'^' = r'^p'cAv{p')9,u■ (6.44) 
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The functions Cij{p'^) are the dressing functions and P^j^ is the transverse projector, P^j^ = 
dfiu — VtiVv/v^- The tree-level expressions are obtained by setting Cyi(p^) = 1, cyip^) = 1 
and CAviP^) = 91^lv^- 

The longitudinal part has to be added to the gluon two-point function in order to 
be able to invert the matrix. Although it may seem that this might provide a direct 
generalization of the Gribov-Zwanziger Lagrangian from the Landau gauge to general 
linear covariant gauge, this is not so. The reason is that the Faddeev-Popov operator 
is no longer hermitian and can have complex eigenvalues |196) . Consequently the usual 
definition of the Gribov region as the set of field configurations with a positive Faddeev- 
Popov operator does no longer make sense. Defining linear covariant gauges in the usual 
way via a Gaussian distribution over the gauge orbit would be the first idea. However, 
while in Landau gauge we are only dealing with a discrete set of Gribov copies inside the 
Gribov region, we have no idea what is changed by smearing out these Landau gauge 
Gribov copies over the gauge orbit. For example, if two Landau gauge Gribov copies are 
very close on the gauge orbit, the region between the two is in a sense preferred, as it 
gets contributions from both Gaussians. It can even be that there are not two maxima 
but only one, if the two Gribov copies are close enough. Such a preference does not agree 
with the original idea of a single Gaussian distribution and its significance is currently 
unknown. A study of Gribov copies for small values of the gauge fixing parameter ^ was 
done in ref. |196] . 

Normally one would expect that the inverse matrix is of the structure 

m4>4>)-i = jjH={^ ^ \ = ^ M ^ C6 45) 

^ ' \£)VA dVV r^Apyy _ (^pAy)2 \ _yVA ^aa I' v"-^^^ 

but because of the additional tensor structure in Lorentz and color space it is not that 
simple. Performing the inversion manually one finds in the employed truncation the 
following propagators: 

jjAA,ab ^^ab]^p Cyjp'^) 



p^ '"'ciip^)cvip') + 2Nclyip^y 



j^VV,abcd _ A^'^A^'^n — fabefcde}_p ^^AVKP ) ( (\ A7\ 

^l^t^ ~ ^2 ^ (^2\" " y/ii^ J / ^2^^^'' ^Lt^2\^2 t^2\ i O /\r ^2 (^2\^ (^2\' l"-^'J 



'^'^ p^cvip^) ^^^ ' ' p^ ^''c\{p^)cl{p^) + 2Nc\y{p^)cv{p^^ 

p' '-'ciip^)cvip') + 2Nclyip^y 



p,AV,abc _ „■ fabc ^ p V^CAvjp ) , . . 
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While, except for some coefficients, the AA and AV propagators indeed resemble the 
conjectured form, the VV propagator comes with two tensors. Thus the employed trun- 
cation leads to a non-trivial structure of the propagators what will be important for the 
IR analysis. 

For the fermionic ghost the standard relation is valid: 

DT' = (nrV = -5^'SJ-^. (6.49) 

The minus was factored out in order to have a positive dressing function, which is required 
for being inside the ffist Gribov region. 

6.2.2 Infrared behavior of the propagators 

In order to assess the behavior at low momenta we make as usual the ansatz that all 
dressing functions have a power law form in the infrared: 

C^AP') = d^Ap'r' ■ (6-50) 

The asymptotic form of the propagators depends on the behavior of the determinant 
ci(p^)cv'(p^) + 2N cj^yi^p^) in the IR. We can think of four possible cases: 

I: c\y > CaCv ^ Ka + Ky > 2kav 

II: caCv > c\y ^ 2kav > i^a + i^v 
III: c\y ~ CaCv <-)■ k^ + kv = ^kav^ no cancelations 

IV: c\y ~ CaCv ^ ha + i^v = 2k av, cancelations 

In case I the dressing function of the mixed two-point function dominates the determinant, 
in case II the combination of gluon and V field dressing functions and in cases III and IV 
all terms contribute equally. The difference between cases III and IV is that it is possible 
that the leading contributions cancel exactly and the determinant is less IR divergent. 
This possibility is considered as case IV. Details on this cancelation are given below in 
Section 16X21 

In the following I will derive the IR behavior of the propagators in all four cases. This 
yields relations between the IREs of the two-point functions and the propagators. With 
this information the IR analysis will be performed based on the results of Chapter |4j 
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Case I: k^ + /ty > ^kav 

By assumption the dressing function cav dominates in the numerator. Therefore the 
propagators take the following form: 

r,AA,ab _ ^ab}^p Cyjp ) 

D^J'^'^' = ^^ {5^'5'%. - r'^^'-P,.}^) , (6.52) 

j^AV,abc ^ _■ fabc]_p 1 _ .g_53^ 

P' "" V2N CAvip') 

Note that both tensors contribute to the VV propagator, since they have the same IRE. 
The relations between the IREs of propagators and two-point functions are 

Sa = i^v - '^i^Av, (6.54) 

Sv = —t^v, (6.55) 

6av = -kav- (6.56) 

These relations allow to calculate the additional term that appeared in the equations used 
in the IR analysis, Aav'- 

Aav = '^^AV — ?>A — ^v = —'^i^Av — f^V + '^I^AV + Ky = 0. (6.57) 

Consequently all equations become the same as for the non-mixing case and the analysis 
is straightforward. 

I will demonstrate now that case I does not allow a consistent solution. For this we 



need the lower bound for IREs given by eq. (A. 19), 



(^"^) (l-^E"^0-^E^^^- (6.58) 



and the inequality derived from the IR leading diagram, eq. (A. 23) 



K, + l5^5,m, -<...,^ ((i " ^""^^ AY.K'""^') ^ 0- 



(6.59) 
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V V-^ ^A A^ 

> — A — • • — V — • > — y — • 

Y Y-^ ^y Y^ 

Figure 6.5: Diagrams of the AA and W DSEs used explicitly in the text. 
The latter formula is used for the W DSE. Note that in this DSE only the bare AVV 



vertex appears, and consequently eq. (6.59) has only one possible realization: 



f^v + ^'26v-^{6A + 26v)-^(^-2]>0 => 2kv - Sa- ^ + 2 > 0. (6.60) 



Using the expression for 6a, eq. (6.54), we get 



2kv - Ky + 2kav - 2 + 2 > => 2kav > -f^v + o ~ 2- (6.61) 

This is the first condition we need for the inconsistency proof. 

Next we consider the AA DSE, to be more precise the VV loop it contains as depicted 



in fig. 6.5 on the left. In the usual way we obtain the inequality 



^ d d 

ka < 2by + - - 2 + kavv => fiA + f^v < -i^v + ^ - 2 (6.62) 

by comparing the left- and right-hand side. We also used kavv ^ from the AVV DSE. 
This inequality constitutes the second condition. 

Now we insert both conditions into the defining inequality of case I: 

d d 

-i^v + -^-2>ka + i^v> 2kav > -«v + 9 ~ 2, 

> 0. (6.63) 

As this inequality cannot be fulfilled, we conclude that case I is not consistent. 

Note that in this proof we did not need the DSE of the AV two-point function explicitly, 
i.e., we did not have to assume that any specific part in its DSE is IR leading. This point 
is important, because the analysis includes also the case that its tree-level part is IR 
leading. Thus there is no need to think about the cancelation of the tree- level part. 
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Case II: ka + k,v < S/t^y (The simple scaling solution) 

According to the definition the combination ca{p'^)cv{p'^) dominates over the mixed dress- 
ing function in this case. The propagators take then the following form: 

K^"'' = ^'''^2P.^;Th-v (6.64) 



P^ c^(p^ 



^VV,abcd I XfJcrM rabe rede p ^AV \.r J \ J^ ^ac^bd, 

'^^ ~p'Cy{p^)\^"' ^ ^ ''ci{p^)cv{p')J P'cvip^ 



(6.65) 



^AV.abc _ ■ fabc}^p V^Cav{p^ 

p2 ^^ c\{p^)cv{p^)' 



At leading order only the first part of the VV propagator contributes, as the second one 
is suppressed compared to the first due to the definition of case II. The relations between 
the IREs are 

5a = -f^A, (6.67) 

Sy = — /ty, (6.68) 

Sav = i^Av — i^A — t^v- (6.69) 

Inserting these expressions into the additional term for mixed two-point functions, Aav, 

we get 

Aav = 2(5^y — 5a — Sv = 2k Av — 2k A — 2kv + ka + kv = 2k Av — ka — f^v > 0- (6.70) 

The last step follows from the definition of case II. As a consequence all diagrams that 
feature a mixed propagator are additionally suppressed, as the coefficient of Aav is the 



number of mixed propagators, see eq. (4.41 ). This also entails that vertices are additionally 
suppressed that require a mixed propagator, as, for example, the AAV vertex. In the 
following these vertices are called V^-odd vertices, as the number of their V legs is odd. 
For the functional equations this means that we can neglect all diagrams containing AV 
propagators and \^-odd vertices in the IR. Additionally, as the V propagator reduces to 
the same form as the U and u propagators, we can subsume all graphs containing rj and 
V fields. The DSEs reduce then to the familiar DSEs of the Faddeev-Popov theory plus 
the DSEs of the AV two-point function and l^-odd vertices. A crucial point is that the 
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Faddeev-Popov DSEs do not depend on the contributions from mixed propagators and 
can thus be solved self-consistently. The mixed correlation function, on the other hand, 
depends on the gluon and V field propagators. Hence, having their solutions, one can 
also calculate them. 

First I will give the results for the diagonal correlation functions and then for the mixed 
one. As we have recovered the standard Faddeev-Popov theory, we can directly use the 
known results [TB| [2Ul \UE\ ESI E21 1147) . The propagators of the Faddeev-Popov ghosts 
and all auxiliary fields are described by the same IRE n := k^ = Kuj = i^u = i^v- In the 
following I denote all these fields as ghost field. The unique scaling relation is given by 



2ka + k = --2 (6.71) 

and the result for an n-point function with In ghost- and m gluon-legs by [721 1147] 

i^2n,m = in- m)K + (1 - n) ( - - 2 J . (6.72) 



This equation is obtained as described in Section 4.2.3 One starts with an appropriate 
two-point function, takes the IR leading diagram and plugs in as many ghost-gluon vertices 
as required. The value of k can be calculated from the leading diagrams in the gluon and 
ghost DSEs using bare vertices as 0.595 353 |29| [65 j . This value was reproduced from 
the equations of the Gribov-Zwanziger theory in order check the code employed in other 
calculations of this thesis. Everything boils down to solve the equation 

(l + /t)(2 + /t) 



12(3 + 4(-2 + fi;)/t) 



1, (6.73) 



which results in k = ^ (93 — vT20l) = 0.595 353 as the only solution with k < 1. 

Next the value of the IRE of the mixed two-point function is evaluated. The diagrams 



that require at least one AV propagator are depicted in fig. 6.6 There is one other 



diagram which has two AV propagators and is thus IR suppressed, see fig. 6.3 The 



diagrams of fig. 6.6 have the following IREs: 



^A + ^AV + i^AAA, 5a + 5y + i^AAV, ^V + ^ AV + ^AVV- (6.74) 
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Figure 6.6: Diagrams of the VA DSE used explicitly in the text. 

kavv is zero due to the scaling relation and the IREs of kaaa and kaav can be determined 
from their DSEs as 

kaaa = 35y, (6.75) 

K-AAv = Sav + 2(5v (6.76) 



Plugging these into eq. (6.74) and using eqs. (6.67) - (6.69) as well as the scaling 
relation eq. (6.71), we self-consistently get Kav for all three diagrams. In the numerical 
calculation of kav we have several unknowns: The coefficients of the dressing functions 
{dAv, dv, dA) and the IREs ny and kav- However, dAv always drops out of the equations 
and dy and dA appear in the combination g'^/dAdy, which can be calculated from the 
Faddeev-Popov part as 0.0267784. That leaves only kav- 

For a practical calculation a further truncation is required and we only take into account 
the AV-VV loop. This is in a sense a one-loop truncation and also employed for the 
Faddeev-Popov theory. There the ghost-gluon vertex is taken as bare, but in principle 
also other diagrams in its DSE contribute at leading IR order. However, if one inserted all 
those IR leading diagrams into the two-point function DSE, one would obtain two-loop 
diagrams. Due to the same argument we only consider diagrams in the truncation with 
an AVV vertex, since vertices like AAV effectively lead to two-loop diagrams. In other 
words we take only into account the order dAdy and neglect everything of order (d^dy)^ 
and higher. However, one should keep in mind that such truncations are only needed to 
calculate a numeric value for the IREs and that the qualitative features and the scaling 
relations themselves are derived for the full system without truncations. 

The equation to determine kav is obtained by projecting the AV-VV loop transversely 
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in Lorentz space and by /"^'^ in color space: 

(6.77) 

The coefficient dAv drops out of tlie equation. Tlie final equation that has to be solved 
numerically is 

_ 0.05130931(1 - fi;Ay)r(0.595353 + kav) 0.0370825r(-KAy)r(0.595353 + kav) 
~ r(1.40465-KAy)r(l + KAy) r(1.40465 - kav)T{2 + kav) 

0.0513093r(-KAV')r(l. 59535 + kav) 0.0210772r(-l - KAy)r(0.595353 + kav) 
r(0.404647 - KAy)r(2 + Kav) 1(1.40465 - KAy)r(3 + kav) 

_ 0.0247217r(-l - fi;Ay)r(1.59535 + kav) 0.0513093r(-l - fi;Ay)r(2.59535 + kav) 
r(0.404647 - KAy)r(3 + kav) r(-0.595353 - /«Ay)r(3 + kav) 

0.0237307r(-2 - /tAy)r(0.595353 + kav) 0.0210772r(-2 - fi;Ay)r(l. 59535 + kav) 

r(1.40465 - fi;Ay)r(4 + kav) r(0.404647 - KAy)r(4 + kav) 

0.0370825r(-2 - fi;Ay)r(2. 59535 + kav) 0.0513093r(-2 - /tAy)r(3.59535 + kav) 
r(-0.595353 - fi;Ay)r(4 + kav) ^ r(-l. 59535 - /€Ay)r(4 + kav) ' 

(6.78) 



The left- and the right-hand sides of this equation are plotted in fig. 6/7 As one can 
see, there are several solutions. A constraint for kav is given by the definition of case II: 
2k, Av > Ky + ka = —i^v + d/2 — 2. The smallest solutions fulfilling this constraint are 
0.0668776 and 0.981386. Further solutions always just a little bit below integer numbers 
follow. At this point we cannot say which one is correct. In our truncation all of them 
are admissible. 



An interesting observation is the fact that the AV-VV loop (see fig. 6.6 ) is automatically 
at leading order: 

Sv + Sav + K.AVV = —f^v + f^AV — Ka — Kv = Kav- (6.79) 

If the tree-level part is kept in the AV DSE and it should be subdominant, kav has to 
be lesser than —1. However, the constraint on Kav is 2kav > — ^y + d/2 — 2 and thus 
Ky > |. Since this is not true for the value determined for Ky, there has to be some 
mechanism that cancels the tree-level term in order to allow this solution. Although the 
details of such a cancelation are not known, it can be expected that it is connected to the 
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1.5 r 




kav 



Figure 6.7: The plot to determ ine solutions for kav- On the y-axis the left- and right-hand 
sides of eq. (6.78) are plotted. 



horizon condition, which is responsible for the cancelation in the Faddeev-Popov ghost 
DSE, see e.g., [56]. For the u, U and V two-point functions the same argument as for the 
Faddeev-Popov ghost can be employed without problems. 

Case III: ka + k-v = '^f^Av< "O cancelations (The strict scaling solution) 

Here both terms contribute equally in the determinant. Abbreviating it as detC = 
caCv + 2Nc\y = {p^f'^^'^detD = {p^Y^+'''' det D , we get 



D^^' =d —P, 



jjVV,abcd 



p2 ""^ (p2)KA det D 

1 1 



s'''6'"'g,, - r'^r^'p^. 



2d^y^ 
detD 



D 



AV,abc 



-^f 



abc 1 p dAV V^ 



P 



,2" ^''' (r^2\K 



(yyAv detD' 



(6.80) 
(6.81) 
(6.82) 



By definition no cancelations in the determinant occur and detD is just a momentum 
independent constant. 

The IREs of the propagators are 



5, 



-I^A, 



(6.83) 
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6v = -Kv, (6.84) 

Sav = -i^Av- (6.85) 



The additional term in the formulas for the IREs, eqs. (A. 17) and (A. 23), becomes 



^AV 



25av — Sa — Sv = —'2k,av + ka + /ty = 0. (6.86) 



Thus again the additional terms vanish, i.e., we can directly get the scaling relation from 
the vertices in the Lagrangian, if we assume that the tree-level two-point function in 
the VV DSE gets canceled in a similar way as that for the other IR enhanced fields. 
Employing the usual technique, we consider possible scaling solutions as obtained from 
the interactions in the Lagrangian. The only non-trivial solution arises from the Afji] and 
AVV vertices: 

K:= Ky = Kj^, Ka + 2k= - —2, KavV = f^Ariij = 0. (6.87) 

The IRE of the mixed two-point function can be calculated from the defining assumption 
of case III as 

'^AV = -^ + ^-l. (6.88) 

In all DSEs the diagrams with a bare Afji] or AVV vertex are IR leading and the formula 
for the IRE of a vertex is 

K2n,m = {n- m)K -^ (l - n) f - - 2 j , (6.89) 

where m is the number of gluon legs and n the number of legs of ghosts or auxiliary fields. 

Thus the qualitative behavior is the same as in case II, i.e., the propagators of the 
Faddeev-Popov ghost and the auxiliary fields are IR enhanced and the gluon propagator 
is IR suppressed, as is the mixed propagator. The difference between case II and III is 
that the mixed propagator is more IR suppressed with respect to k in case II, since there 
it holds that 5av > i^/2, whereas in case III we have 6av = i^/2. 

Having found a scaling relation does, however, not yet mean that we can be sure there 
really exists a corresponding scaling solution. A first indication of its existence is a solution 
for the IRE k. This is an intricate task as the system of equalities that has to be solved 
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involves the coefficients of the power laws, the dij, in a non-linear way. In case II only the 
IRE n appears non-linearly and the remaining equation could be rewritten into a simple 



form, see eq. (6.73). For the present system, however, this is not possible. 

For case III it is necessary that the tree-level term in the AV DSE gets canceled by 
some mechanism that may be related to the horizon condition. If this is not the case the 
AV-VV loop can only be leading for k > 2 or - if the tree-level term leads - we get k, = 2. 

Case IV: ka + ny = 2kav< cancelations 

The final case that remains to be investigated has the same condition on the IREs of 
the two-point functions as case III, i.e., both terms in the determinant scale equally. 
However, even the coefficients of the power laws are the same so that the leading terms 
of each expression cancel each other: 

dAdv = -2Nd\y. (6.90) 

The IRE of the term that takes over is undetermined, so we introduce an additional IRE 
kd which gives the correction of the IRE of the determinant: 

det C =caCv + '2Nc\y = 

= {p'Y^^^-{dAdv + 2Nd\y + doiv'T^ + ...) = 
= (p')"-^+"^(t/D(/)"^ + ...). 



(6.91) 



Consequently n£, is non-negative, Ho > 0. 

The propagators are then almost the same as in case III, except for the VV propagator, 
where now the ffist term is suppressed because there is no determinant: 



^AA,ab _ <:ab ]^ p "'V 

M- p2 ^'^{p2YAdetD 



^AAa^^^a. p^^ V_ (g_g2^ 



T^VV,abcd _ 1 ^ I A^^A^'^n _ fabercdep ^^AV \ _, ^ ^ ( _ rabe rede p ^'^AV 

^^ ~ p^ cvip^) \ ^^^ ^ ^ ^"detDj^ p^cv{p'')\ ^ ^ """detD 

(6.93) 



^AV,abe _ ■ fabe]^p "AV v2 

p2 ^'^ (p2)KAV det D 



JjAV,abe ^ _^ fabe p^^A^y^_ (g^g4) 
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The IREs of the propagators are in this case 

5a = -K,A - K.D, (6.95) 

6v = — Ky — i^D, (6.96) 

6av = -hav - f^D- (6.97) 

Hence we have again Aav = '^^av — ^a — ^v = 0. 

The appearance of the new IRE k^) does not influence the general analysis as described 



in Section 4.5 , because it only appears in the relation between the IREs of propagators 



and two-point functions, eqs. (6.95) - (6.97). So we can employ in the usual way eq. (A. 23) 



for the VV two-point function, where only a bare AVV vertex can appear: 

Kv + <5y - 2 {^A + 2(5y) >0, 

KA + 2fi;y + kd >0. (6.98) 



On the other hand, from counting the VV-AA loop in the VV DSE, see fig. 6.5 on the 
right, we get another bound: 

K.V <'^y + 8a + i^avv , 
KA + 2kv + 2kd <0, (6.99) 

where we used kavv ^ from the DSE of the AVV vertex. Note that these two inequal- 
ities look alike except that the coefficients of kd are different. Combining them yields 
therefore additional information on kd'- 

-kd < ka + 2kv < -2k£,, 

t^D > 2kd. (6.100) 

Since no is non-negative this inequality yields kd = 0. This corresponds to case III and 
thus case IV does not yield a solution. 
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6.2.3 The infrared behavior of the propagators in the 
Gribov-Zwanziger theory 

Although the analysis of the Gribov-Zwanziger action was complicated by the mixing of 
the gluon field A and the auxiliary field V, it was possible to obtain a scaling relation. 
The most important observation is certainly that the qualitative picture for the behavior 
of the gluon and ghost propagators is not changed, as the scaling relation 6a = — 2(5c 
is still valid, i.e., the gluon propagator is IR suppressed and the Faddeev-Popov ghost 
propagator IR enhanced. Furthermore, the IREs of the propagators of the auxiliary fields 
u, U and V are the same as for the Faddeev-Popov ghost propagator. The two solutions 
obtained differ mainly in the behavior of the mixed propagator: Though IR suppressed 
in both solutions, this suppression is more pronounced in case II, where 6av > '^/2. In 
case III the IRE 6av is determined by the scaling relation as 6av = '^/2. 

One of the solutions found was not accessible to a numerical solution for k, so it remains 
unclear, if it really exists. For the other solution, case II, the calculation of k, was unex- 
pectedly easy as the system of DSEs reduces in the IR to the system of the Faddeev-Popov 
theory. Due to this one can directly employ all the results known from there. The mixed 
propagator does not influence the Faddeev-Popov system and its IRE can be calculated 
with input from there. Several allowed solutions were found. 

The obtained results confirm the conjecture by Zwanziger that the cutoff of the inte- 
gration at the Gribov horizon does not influence the results from DSEs, but it has to 
be taken into account via choosing appropriate boundary conditions |29j. The boundary 
condition adopted here is derived from the horizon condition and allowed the IR enhance- 
ment of the propagators of the ghosts and the auxiliary fields. It is expected that the use 
of an alternative boundary condition will lead to a decoupling type of solution as in the 
Faddeev-Popov theory [19]. It is possible to arrive at such an alternative condition by 
adding dimension two condensates to the Gribov-Zwanziger action [961 [971 1197]. Their 
effect is that no longer the Gribov copies directly at the horizon dominate in the path 
integral but gauge field configurations lying within the Gribov region. This should be in 
direct relation to the boundary condition imposed on the DSEs of the two-point functions 
of the ghost and the auxiliary fields. 
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Conclusions 



The central topic of this thesis has been the IR behavior of correlation functions in dif- 
ferent gauges. The IR regime is directly linked to the property of confinement of the 
elementary fields. Prominent confinement scenarios are the Kugo-Ojima or the Gribov- 
Zwanziger scenarios, which are based on completely different arguments but lead to the 
same conclusions in the Landau gauge. Also the hypothesis of Abelian IR dominance, 
motivated by the dual superconductor picture of confinement, is related to IR proper- 
ties. Consequently it is of interest to determine the behavior of propagators and vertices 
at low momenta and test the different scenarios where applicable. Furthermore, know- 
ing the relations between different gauges further elucidates the overall picture of the 
non-perturbative regime of Yang-Mills theory. 

A significant obstacle in the investigation of other gauges than the Landau gauge is the 
complication of the structure of the interactions between the fields. This makes their IR 
analysis a very tedious task with the conventional techniques. In Chapter |4] I described 
an extension of these methods that allows to assess the existence of different IR solutions. 
Although the proof is rather technical, the main statement is rather simple: One can 
infer possible IR solutions directly from the interactions of the Lagrangian. This is a 
remarkable result insofar as one can derive the IR behavior of the whole infinite tower 
of functional equations. The derivation of the known Landau gauge results reduces with 
this method to a few lines. Further gauges shortly investigated in Chapter |4] are linear 
covariant and ghost anti-ghost symmetric gauges. 

Another complication of non-Landau gauges is that the mere derivation of the DSEs 
becomes very time consuming. The natural consequence was to do this task with the help 
of a computer algebra system and led to the development of the Mathematica package 
DoDSE. Its use is not restricted to the purposes of this thesis, but it can derive the DSEs 
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for arbitrary actions. The underlying algorithm is also suited for manual calculations and 
described in Chapter [3j Having the DSEs directly available on the computer alleviates 
many task as, for example, the calculation of IREs. The length of the resulting expres- 
sions makes clear that manual calculations are extremely difficult when going beyond the 
Landau gauge. 

The maximally Abelian gauge was investigated in Chapter [5j Due to the complexity 
of this gauge the improved method for the assessment of the IR behavior of correlation 
functions was of great help here. The main results were an IR enhanced diagonal gluon 
propagator and IR suppressed off-diagonal propagators. This very nicely supports the 
hypothesis of Abelian IR dominance according to which the off-diagonal degrees of freedom 
should be suppressed compared to the diagonal ones in the IR. It has to be mentioned, 
however, that this result is not in direct agreement with current lattice simulations. The 
reason for this disagreement is most probably the same as in the Landau gauge, i.e., 
two solutions exist which depend on the boundary conditions imposed on the system of 
functional equations. For the maximally Abelian gauge this would be implemented by 
the renormalization condition for the diagonal gluon propagator. The scaling relation 
of the maximally Abelian gauge constitutes the first consistent solution of the system of 
infrared exponents that was found beyond the Landau gauge. Furthermore, the connection 
between ghost dominance in the Landau gauge and Abelian dominance in the maximally 
Abelian gauge is further elucidated by these results. Also on the infiuence of the shape of 
the Gribov region, which is different for the Landau and the maximally Abelian gauges, on 
the IR behavior of Green functions was speculated. Finally, it should be mentioned that 
this work was the first one to investigate the maximally Abelian gauge for the physical 
gauge group SU{3). As the action is different in SU(2) and SU{N > 2), it is rather 
non-trivial that the obtained IR behavior of propagators and vertices is the same in both 
cases. 

Chapter |6] concerned the non-perturbative gauge fixing in Landau gauge. A common 
way to overcome difficulties due to Gribov copies to a certain extent is to use the Gribov- 
Zwanziger action. Its IR behavior was determined with the methods developed in Chapter 
|4j The semi-perturbative results of an IR suppressed gluon and an IR enhanced ghost 
propagator were confirmed with this non-perturbative analysis. At the same time the 
qualitative behavior is the same as obtained from functional equations when the standard 
Landau gauge fixing is employed. This corroborates the conjecture by Zwanziger that 
only the proper choice of the boundary condition is relevant for the functional equations. 
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as one can formally cut the integration at the Gribov horizon. 

The methods developed in this thesis allow first steps in the investigation of Yang- 
Mills theory beyond the Landau gauge. Due to the inherent complications the following 
two tools are of great help: An automated derivation of the functional equations and a 
method for a qualitative determination of the IR behavior capable of dealing with many 
interaction terms. Their formulations are not specific to the problems investigated in this 
thesis but have an even greater area of applicability. Thus they form a sound basis for 
further investigations to improve our understanding of the strong interaction. 
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Appendix A 

Details on the formulae for the 
scaling analysis 

I present in this appendix the derivations of several equations of Chapter |4} All calcu- 
lations are done in d dimensions. Another generalization is that mixed propagators are 
taken into account wherever necessary, but the case of a diagonal two-point matrix can 
be directly inferred from the results. 

A.l Inequalities from DSEs and FRGEs in d 
dimensions 



For completeness I shortly summarize the derivation of inequality (4.16) in d space-time 



dimensions which was already given in four dimensions in Section |4.2.1 

We start with three-point functions. From the FRGE of a generic three-point function, 
as depicted in fig. |A.1[ one can extract the inequality 

(A.l) 



For four-point functions we get from fig. A. 2 the following two inequalities: 



-5a-6b-[--2\<kaabb. (A.2) 

i^AABB < 2kabcd + 5c + 5d + --2. (A.3) 
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Figure A.l: One specific diagram in the FRGE of a generic tliree-point function. Internal 
lines represent dressed propagators, black blobs dressed vertices. The grey 
blob is a regulator insertion. 



Combining them yields 



KABCD + ^{Sa + Sb + 5c + Sd) + 2 I-- 1 ] > 0. 



The inequalities for the three- and four-point functions fulfill 

(^ -l)(r -2) + /.,...,, + ^^A;---5,>0, 



(A.4) 



(A.5) 



where r = 3, 4 denotes the number of legs the n-point function has. Higher n-point 
function also obey this inequality as can be shown by induction. 



For the proof we need the following two inequalities which can be inferred from fig. |A.3 



I^ABCD- < I^ABCDEE- + 0_E + ;^ — 2, 
K-AABBCC- < 2kaBCDE- + 5d + 5e + - — 2. 



(A.6) 
(A.7) 



a) 









b) 



Figure A. 2: Parts of the FRGEs of generic four-point functions. 
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k 



dk 





A / \ D A / \ D 

a) B C B C 



b) 





Figure A. 3: Parts of the FRGEs of generic n-point functions. 



The dots represent further legs as indicated in fig. A. 3 Note that these inequahties are 



generahzations of eqs. (A. 2) and (A. 3). The first inequahty can be used to write down 



an equation for the vertex AABBCC as appearing in the second inequahty: 



I^AABB- 



[l^AABB- 



^C 



d \ d 

- — 2 I < KaABBCC- < ^I^ABCDE- + (JD + ^S + ^ 



^C 



5d - 5e) 



d 



— 2 < kabcde- ■ 



(A.8) 
(A.9) 



This inequahty connects an n-point function with a (2n-6)-point function. The goal is to 
rewrite the equation such that only propagator IREs and the IRE of the n-point function 



remain. For this one can successively use eq. (A. 6) to replace the IRE of the remaining 
other vertex and arrives at 



-^ i^A + Sb + Sc + 6d + Se + . ■ ■) - [ - - I ] (r - 2) < kabcde 



(A.IO) 



where r gives the number of legs of the n-point function. The coefficient of (i/4 — 1 is 



obtained as follows: One gets 2{d/A — 1) in eq. (A.9) and d/A — 1 from every further 
application of eq. (A. 6). The latter has to be done until 2ra — 6 — 2a; = 2, where x is 
the number of iterations. Hence the total coefficient of d/A — lis2-|-x = n — 2. This 



establishes eq. (A. 5) for all n-point functions. 



For the second group of inequalities the non-positivity of the IREs of interactions ap- 



pearing in the Lagrangian, eq. (4.20), is used in eq. (A.5): 



I 1 I (r — 2) H — 2_, k^/"'^''Si > V primitively divergent vertices. (A. 11) 
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A. 2 Infrared exponent for an arbitrary diagram 

It is possible to express the IRE k^ of an arbitrary diagram v by counting the IREs of 
all its propagators and dressed vertices. For d dimensions also the canonical dimensions 
have to be taken into account: 

i vertices ,r> 3 

+ Yl <...vCn...v-c.. (A. 12) 

vertices,r>3 

Hi are the number of internal propagators with IR exponents (5j, whereas the numbers of 
vertices (pi^- ■ ■ (pi^ are ni^,,,i^. Superscripts d and h stand for dressed and bare, respectively. 
In case none is given, I refer to both. The sums X]i;ertjces r>3 extend over all vertices with 
r legs. TJii is the number of external legs of field type 0j and / is the number of loops. The 
canonical dimensions of the vertex 04i ■ ■ ■ 0j,, are given by Ci^,,,i^. 

Without mixed propagators it is possible to use topological relations to get completely 
rid of the internal propagators. In the more general case, however, when mixed propaga- 
tors appear, the dependence on them will remain. To avoid too cumbersome notation we 
restrict ourselves to the case of two fields only. We denote them explicitly by A and V , 
as these are the fields for which the mixing occurs in the Gribov-Zwanziger Lagrangian; 



see Section 6.1.2 for details. The case of no mixing fields can be obtained by setting the 
number of mixing propagators uav to zero and extending the sums over all fields in the 
formulae to come. Moreover, we have 5i = — Kj in this case. 

The first topological relation we need gives the number of loops in terms of the numbers 
of propagators and vertices: 

1= X^ ^i + iT-AV + 1 - X] ni^...ir- (A. 13) 

i=A,V vertices,r>3 

Furthermore, we express the number of internal A- and V^-propagators by 

_ I Y^ k^-'^^Ui,,,,,^ -m.-nAvV i = A,V, (A. 14) 

\vertices ,r>3 / 



2 



where /c*^'"*'' denotes the number of times the field 0j appears in the vertex (pi^ 
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Plugging these expressions into eq. (A.12), we get 

5Z 2 ( 5Z ^r^'^'^n-v -mi-UAvj +'i- + nAv- Yl "'^i-^'- ) 2^ 

i=A,V \vertices,r>3 / vertices,r>3 J 

+ nAv{^Av - 1) + XI 2 ( ^ ^i'"'^ni^...ir - ^i - hav J {5i - 1) + 

i=A,V \vertices,r>3 / 

vertices,r>3 vertices,r>3 i=A,V 

^ -^ \ i=A,V / i=Ay vertices,r>3 ^ -^ 

+ E <...v (-^ + ^n...v + 2 - + n^^ ( (5^v - 1) - ^ E (^^ - 1) ) + 

vertices,r>3 ^ \ i=A,V / 

+ E ^ f E ^r-"^n...v(^ + '^^ - 1) ) , (A.15) 



2 , 

i=A,V \vertices ,r>3 



vertices 



where it was used that the canonical dimension of a vertex is given by (4 — r)/2 with r 
being the number of external legs. Reordering terms yields 

K„ = (--2J [1-2 E ^A ~ 2 ^ ^'^' ^ 2'^av{'2Sav -Sa- Sv)+ 

^ ^ \ i=A,V / i=A,V 

E <...v f (i - 1) (^ - 2) + '^n...v + \ E ^"•■■"^^) + 

:ices,r>3 \ ^ ^ J=A,y / 

+ E <...v f (i - 1) (- - 2) + 1 5^ ^r-^^^.) • (A.16) 

vertices,r>3 \ ^ '^ i=A,V / 

This is the formula for the IRE of the diagram v with rrii external legs of the field (/)j and 
the number and type of vertices given by the n'' and n'^. The dependence on the internal 
propagators is given via the term UAvi'^^AV — ^a ~ '^v)/2- The maximally IR divergent 
solution in the d-dimensional case is 

K.v,max = (2~2j (^"o E"^M~2 E ^^^' + ^nAvC^SAV - Sa - Sy). (A. 17) 
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It is derived from eq. (A. 16) and the non-negativity of the last two terms, see eqs. (A. 5) 



and (A.ll) 



The corresponding formulae in the case of no mixing fields are obtained by setting 
uav = 0: 



and 



vertices,r>3 \ i / 

vertices,r>3 \ i / 



(A.18) 



(A.19) 



for the maximally IR divergent solution. 

A. 3 Inequality from the leading diagram 

Under the assumption that a certain diagram is leading in a two-point DSE one can 
derive an additional inequality constraining the system. For this one starts with the 



general expression of a two-point function IRE, obtained from eq. (A. 16): 



r\jn 



+ E <...v-n...v+ E -n...v((^-l)(r-2) + i5^fc--5.] 

dressed all \ j / 



(A.20) 
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A. 3 Inequality from the leading diagram 



Here Kj can be any two-point IRE and 6j is restricted to IREs of non-mixing propagators. 
Using V . mj = 2 in the case of propagators, this can be written as 

3 
all \ j / dressed 



(A.21) 



vertices 



vertices 



We can use the lower bound for the IREs of the vertices on the right-hand side given by 



the maximally IR divergent solution, eq. (A. 17), to get a new inequality: 



'^^ + 9 E ^^'^^ ~ t^'^avC^Sav 



— Sa — Si 



> 



vertices 

E<..,((^^)^45:*H- 

ressed \ ^ ^ \ j / 

irtices 



(A.22) 



Here fc^y *'' indicates the number of times a mixed propagator is contained in the diagram 
that determines the IRE of the vertex (pi^- ■ ■ (pi^. Note that fc^^*'' only is different from 
zero for vertices that necessarily contain an AV-propagator like the AAV"-vertex. The 
right-hand side depends on dressed vertices only, indicated by the superscript d of n. On 
the other hand, the left-hand side sums over dressed and bare vertices, so that in total 
only the bare vertex remains in the sums over vertices: 

'^^ + ^ E ^^■"^^- - <...v ( (^ - 1) (^ - 2) + ^ E ^r •■"^^) - 



M 



AV 



>V 



Uav 



+ 2^ ''^il...ir^AV 



\ 



dressed 
vertices 



\ 



J 



>0. 



(A.23) 
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Here J2i ^f *'^ — ^ ^^^ used. 

Again the case of a Lagrangian with only diagonal two-point functions is obtained by 
setting uav and consequently also fc^y*'' to zero. Additionally we have ki = —6i so that 
the first two terms cancel: 



Combining it with eq. (A. 11) we obtain the equation 

<.... ((i-l) (r-2) + ^5^fc--5,') =0, (A.25) 

where the hat indicates that it is only valid for a specific n-point function. 
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Grassmann fields 



In many quantum field theories not only commuting, but also anti-commuting fields ap- 
pear, be they fundamental particles like electrons or quarks or mathematical constructs 
like ghosts. These Grassmann fields entail additional complications as one has to care- 
fully take into account their anti-commutativity. For their treatment one has to choose 
one of several conventions concerning the direction of derivatives and the ordering of the 
fields. It is also possible to introduce a metric, see, for example, ref. |126| . This ap- 
pendix summarizes the convention employed in this thesis, which has also been adopted 
for the program DoDSE. In the following Grassmann fields are denoted by tp and their 
anti-fields, here called anti-Grassmann fields, by ip. The corresponding sources are fj and 
r], respectively. 

Derivatives with respect to anti-commuting fields or sources are defined as acting from 
the right and left for Grassmann and anti-Grassmann quantities, respectively: 

A .= 1 A .= 1 m n 

Hence by definition all derivatives are written at the left side of an expression but act from 
the correct side. Consequently quantities should always be ordered such that Grassmann 
fields are right of anti-Grassmann fields. This is valid for the products of fields and sources 
appearing in the path integral as well as for derivatives. So the path integral for a theory 
with the two Grassmann fields ip and ip reads 

Z{r],f]]= /"r'[V'^]e-^[^'^]+*''^+'^'' (B.2) 
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and a quartic Grassmann interaction has the form 



ijkl 



di/jidipjdipkdiJi 



- ■ (B.3) 

ip=4)=o 



For easier readabihty (and also in correspondence with DoDSE) the indices of Tij^i do not 
reflect the order of how the derivatives are performed, but rather have the order in which 
the derivatives appear, i.e., the differentiation with respect to ipj has to be performed 
before that with respect to ipi, but ipk comes before t/^;. 

The replacements of the fields to obtain the generating DSE for IPI functions, see 



eq. (3.14), are given by 



i'r^^r + Dtf"^, (B.4) 

i,^ ^ i,^ + Dif" ^^ (B.5) 



for anti-commuting fields. Furthermore, eqs. (3.16) have to be amended by 

_r,J - r,J pJ,Vi**r)J m fip^ 

X 7 jk -^jm^ imn -^nki yi->.\ja) 

^Dj, = D^^Tilf^Di,, (B.6b) 

^^i-.. = r/...„., (B.6d) 

where the fields were added as superscripts for easier identification. Care has to be taken 
for the internal super fields $, which can be commuting or anti-commuting. As long as 
it is not determined if they are Grassmann or anti-Grassmann, i.e., if the derivatives are 
right- or left-derivatives, respectively, they should be considered as floating. Only at the 
end, when the sources are set to zero and the super fields become physical fields, one 
places them at the left or the right side, but one still has to obey the order in which the 
derivatives have been applied. This is an important point as the required ordering yields 
the signs expected normally for Feynman diagrams with fermion loops. 

An example of how to order the fields should help to illustrate this point. We consider 
the quark respectively ghost loop in the gluon DSE of Landau gauge where ip = {q,c}. 
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After performing the first derivative with respect to the gluon field Ai, we have to replace 



the fields as indicated in eq. (3.14). For the ant i- Grassmann field this is ipr ^ i^r + 

-Sti'^Ai^s ^ -Sti'' {Ai^s + Dif^) . (B.7) 

Differentiating once more with respect to Aj yields 

5T 






+ gluonic terms, (B.8) 

j=o 



where the external sources have already been set to zero. Ordering the derivatives to get 
the canonical order changes the sign of the expression and leads to the expected relative 
minus sign of closed fermion loops: 

rf = Stt^'DitDtXi:^, + gluonic terms. (B.9) 

Finally, we consider the expansion of the action when Grassmann fields are involved. 
First, by definition all ant i- Grassmann fields have to be left of the Grassmann fields. 
Second, the expansion coefficients are antisymmetric in the indices belonging to anti- 
commuting fields. This entails that we can differentiate with respect to these fields as 
usual, e.g., 

^siit^'AAMu =stt:!'^sMu - sS:!:^rAA = 

=stsi:^HsAA + sti:^H^^,^^ = 2strHsAA. (b.io) 

The additional rules for Grassmann fields described in this section allow their inclusion 
in the derivation of DSEs, but also for FRGEs this convention is adequate. 
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Appendix C 

Feynman rules of the maximally 
Abelian gauge 



I collect here the Feynman rules of the MAG. The action is given in eq. (5.31) with the 
interpolating gauge fixing parameter chosen as C = 1 and the gauge fixing parameters a 
and ^ unfixed. The indices i, j, k are diagonal and a, 6, c, c?, e, off-diagonal; r stands for 
both. The momenta are always chosen as ingoing, except for anti-ghosts. The momentum 
convention for the Fourier transformation to momentum space is defined with a positive 
sign for ingoing momenta, e.g., 

rAt!?(P3,Pi,P2) = JdxdyT{z,x,y)e'^^P-^+yP^-^P^^ (C.l) 

for a ghost-gluon vertex, where p2 and ps are incoming and pi is outgoing. The order of 
indices corresponds to the order of fields given as subscripts indicating the vertex. 
The bare propagators are 

Df'\p') = 5^^-^ (^,. - (1 - 0^) , (C.2) 

Df^'\p') = S^"^^ (^g,^ - (1 - a)^) , (C.3) 

Df '(0)(p2) = _<5«fel. (C.4) 

The bare three-point vertices are 
^abbLMp^'P^) = 9^r'\2ny6{p^+p,+ps)x 
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X f ^M^(P2 - Pl)p + gfip{pi - P3)u + gupiPS - P2)t, + -idupPl^ - 9ppP2^ 

'(C.5) 
X {g^,v{p2 - Pi)p + gpp{pi - P3)u + gupips - ^2)^) , (C.6) 

^7SiP3,PuP2) = ^^r'^(2vr)^5(-pi +P2 +P3)(P2 +Pi)m, (C.7) 



rst!^P3,Pi,P2) = gtr'%2ny6{-p^+p2+P3) (^^3+^2) . 



(C.^ 



X 



Finally, the bare four-point vertices are 

^fABB,,upAP3,P^,Pl,P2) = g\27rf6{p,+p2+P3+p,)nf=' 

X ( 2g^^gp^ - ( 1 - - J gpag^p - ( 1 - - ) d^ipg^a j , (C.9) 

^ABBB,pupaiP4,Pl,P2,P3) = /(27r)^5(pi + ^2 + Ps + P4) X 

X {f-'^f'^'^ {g^^g,, - g^^g.p) + /^-/«- {g^,g^, - g^^g,,) + /-^^^ {g^^g.^ - g^,g^^)) , 

(CIO) 

^BBBB,pupaiPl^P2,P3,P4) = /(27r)^5(pi + ^2 + Ps + P4) X 

X (/'^'V^'^ {gppgua - gp.g.p) + /'^7'^*' {g,ugp. - g^pg..) + nT' {gp.g.p - gp.gp.)) , 

(C.ll) 
1 

rSct(P3,P4,Pl,P2) = /(27r)45(-pi +P2 +P3 +P4)2(?^.r7'^^ (CIS) 

r?S,t(p3,P4,Pi,P2) = -g\2T,)%-p,+p,+p,+p,)g,, {r'^r'^ + r'f'') , 

(C.14) 

^ltcc:pJ{Pl,P2,P3,PA) = -/(27r)^(5(-pi -P2 +P3 +P4)X 
y _ f A eabk rcdk < c\ rabe rede < cade rbee cace cbde\ //"< 1 r \ 



^TBiUp4,P3,Pi,P2) = g\2n)%-p,+p,+p,+p,)g^,- (2r^7«- - r'^r^) , (c.12) 
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Calculation of the sunset diagram 



The solution of the sunset diagram is 

]d„ ,]d. 



=(4vr)-'^(p2) 



d. 2,d+a+6+er(a + rf/2)r(6 + rf/2)r(c + d/2)Ti-a -h-c-d) 



r(-a)r(-6)r(-c)r(a + h + c + 3d/2) 

(D.l) 

It can be derived, for example, with NDIM or by using the solution for the one-loop 



two-point function, eq. (4.3). Raising the number of loops by adding further propagators 
between the two vertices one gets the so-called water melon diagram for which the solution 
is also known, see e.g., ref. |198] for the massive case. Introductions to NDIM can be 
found, for example, in refs. ^M Iffl 1X99] . 

The main complication is that we need not only this integral but also the case where 
factors oi pq, pr and q r appear. In contrast to the one-loop case it is here not possible 
to express the scalar products by the invariants appearing in the integral. Following an 
idea by Suzuki and Schmidt |135j . who calculated the integral when arbitrary powers of 
(2pq) appear, one can use NDIM for the calculation of 

Issia, b, c, e, /; /) := J ^^(9')1(0?[(P " Q " r)2]'=(2pg)^(2pr)/. (D.2) 

This indeed suffices as the scalar product between the two loop momenta can be expressed 
via the other scalar products: 

qr = - ((p — q — t)^ — p^ — q^ — r"^ + 2pq + 2pr) . (D.3) 
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Following the NDIM procedure one obtains 276 solutions in terms of five-dimensional 
hypergeometric series, which should all be equivalent. We exploit this huge number by 
choosing one solution where the hypergeometric series terminates: 

Issia,b,c,e,f;p'^) = 

^ .. .^a+6+c+.+e+/oe+f (-«> 2a + 6 + d/2){-b, 2b + f + d/2){-C, 2c + d/2) ^ 

^^ ' -a-b-c-d,2a + 2b + 2c + e + f + 5/2) 

"y^^^ (_l)-"4-«82»6(c^^/2,n3 + n7) 

^^ nsln^lriQlrijlnslil — a — d/2 — e,n% + n^ — ng) 

^ P(e, ng + n4, nQ)P{f, rij + ng, tiq) 

(1 — 6 — d/2 — f,nj + ns — n4) (1 + a + b + c + d, n^ + n^ + uq + rij + ng) 

The symbol P is defined as 

P{a, 6, c) := (-a/2, 6 + c/2)(l/2 - a/2, 6 + c/2) (D.5) 

and 

r(a) 

is the Pochhammer symbol. Although a five-fold series is normally quite a nasty object 
to deal with, this is not a real problem here. The reason is that the series is truncated 
as we can see from the definition of the symbol P. Depending on if the first argument is 
even or odd, P vanishes from certain values of the summation variables on, for example, 

P(e, ns + n4, ne) = for e < 2n3 + 2n4 + uq if e/2 G I A ne/2 G I. (D.7) 

This can most easily be taken into account by restricting the sum to max[e, /]. As e and 
/ are typically quite small this is sufficiently fast in a numerical implementation. 
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Appendix E 

Deriving Dyson-Schwinger equations 
with DoDSE 



In the course of this thesis a Mathematica package [llOj was programmed to automate 
the derivation of DSEs. As we turn to more and more comphcated gauges such a tool 
is indispensable, because the numbers of diagrams become very large. DoDSE, what is 
short for "Derivation of Dyson-Schwinger equations", can handle complicated Lagrangians 
easily by resorting to a symbolic notation. This alleviates the process of the derivation 
considerably. After the DSEs have been derived in this symbolic notation the individual 
terms can be transformed into the proper algebraic expressions. 

Available documentation for DoDSE include ref. |lU8j and a Mathematica notebook 
distributed together with the package, which, among other places, can be obtained from 



http://cpc.cs.qub.ac.uk/summaries/AECT_vl_0.html Furthermore, the syntax and ex- 



amples for all commands are listed via the Mathematica command ?, e.g., ?doDSE. The 
present Appendix describes some general aspects of the program and provides examples 
used in the calculations of Chapters [5] and |6| The final section contains a list with all 
functions of DoDSE. 

E.l Overview over DoDSE 

The main challenge in doing functional calculations with a symbolic programming lan- 
guage is to implement the properties of fields, propagators and vertices properly. A very 
convenient way for this is provided by the index notation introduced in Sec. |3j Thereby all 
indices, e.g., color, Lorentz, Dirac, and the space or momentum dependence are described 
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by only one index. This suffices as long as all quantities like color or momentum "flow" 
through the diagram together. The evident advantage is the reduction of redundancy in 
the notation. 

Instead of using the standard multiplication operator DoDSE employs its own function 
called op. The main reason is that the order of quantities is not changed unpredictably 
by Mathematica within a calculation and also the anticommutativity property of Grass- 
niann flelds can be taken into account. The arguments of an op function can be flelds, 
propagators or vertices. Fields are defined as a list with two entries: The first determines 
the name of the field and the second gives its index, e.g., {A,i}-. If a field appears as 
an argument of an op function it is an external field, n-point functions exist in bare and 
dressed form, S and V, respectively. Dressed propagators are denoted by P. Note that in 
the derivation of DSEs there is no need for a bare propagator but only for a bare two- 
point function, which is denoted by S. The arguments of V, P and S are fields. A simple 
example containing an external field, dressed propagators and bare and dressed vertices 
is the expression 



op[{A,i}, 
P[{A,jl},{A,j2>] , P[{A,kl>,{A,k2>] , 
S[{A,i>,{A,jl>,{A,kl>,{A,a}] , 
V[{A,j2>,{A,k2>,{A,b}]] 



It is depicted in fig. |E.1[ All internal indices are summed and the free indices a and b 
correspond to the external legs of the graph. 

For using DoDSE the package file DoDSE. m has to be loaded. The standard way is 

«DoDSE' 

if the package resides in the subdirectory DoDSE of $UserAddOnsDirectory. 

The next step is the definition of an action. For the simplest cases DoDSE only needs 
its basic structure given by lists of fields. The action of Yang-Mills theory, for example, 
is defined as 

actionYM = {{A, A>, {c, cb>, {A, A, A>, {A, cb, c>, {A, A, A, A»; 

From this definition DoDSE will automatically infer that the ghost and anti-ghost fields 
c and cb, respectively, are Grassmann fields as they appear in a pair. If this assump- 
tion is unwanted, because the fields are not anti-commuting, one can provide the option 
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Figure E.l: The graphical representation of the example provided in the text. The circle 
denotes the external field. 



specif icFieldDefinitions to declare fermions and bosons specifically. An example to 



this is given in Section |E.3.1[ A DSE is then derived from the action with the function 
doDSE. Its two required arguments are the action and the list of fields with respect to 
which the derivatives are performed. For example, the two-point DSEs in Landau gauge 
are obtained by 

AADSE = doDSE[actionYM, {A, A}]; 
ccDSE = doDSEEactionYM, {cb, c}] ; 

The result of doDSE is a (sum of) op function(s) each representing a single diagram. 
There are basically two ways of getting this into a useful representation. One can us the 
function shortExpression (or sE) to get the result in a shorthand notation where A 
represents a propagator, F a dressed vertex and S a bare vertex. For the gluon two-point 
DSE this looks like 

sE [AADSE] 



qAA_WcAAAAaAA\ _WcAAApAAA aAAaAA\_ 
^il i2 2 V il i2 rl sl^rl si;! 9 v il ^l ^^ i2 tl ul^rl tl^sl ul) 



/oAAAApAAAA 
Wil rl r2 si"*- i2 s2 t2 u2 



aAA aAA aAA \ _i_ qA cb c pA cb c a c cb a c cb 
^rl s2^r2 t2^sl u2) "r "-"il rl si"*- i2 tl ul^sl tl^ul rl' 



/ cA AAApAAApAAA aAAaAAaAA aAA\ 
Wil rl r2 si""- 12 s2 tl""- ul v2 wl^rl s2^r2 v2^sl wl^ul tlj 
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The superscripts denote the fields and the subscripts the corresponding indices. 

It is also possible to use DoDSE for drawing Feynman graphs. This is done with the 
function DSEPlot: 

DSEPlotEAADSE, actionYM, {{A, Red>, {c, Green»] 

As second argument one has to give the action. The third argument is optional to give 
some graphics primitives for the different fields. 

I would like to make a short comment on the graphical representation of DSEs using 
DoDSE. DSEPlot employs the Mathematica function GraphPlot that originally is intended 
for drawing graphs in graph theory. However, it is possible to abuse it also for drawing 
Feynman diagrams, but with a few limitations: It is not possible to use wiggly lines and the 
alignment of some graphs may be arbitrary. And when it comes to non-planar diagrams 
the usefulness of GraphPlot comes to an end, as the result cannot be distinguished from 
a planar diagram. Thus the ability to represent DSEs graphically is limited, but with a 
little bit of effort it can produce presentable figures. 

In order to get the algebraic expressions with all the indices and the integrals one uses 
the function getAlg. Before it can be used one has to specify the Feynman rules. These 
are defined globally by overloading the propagator and vertex functions. For example, the 
gluon propagator, the ghost propagator and the ghost-gluon vertex in the Landau gauge 
are defined as 

S[c[cl_, pl_] , cb[c2_, p2_] , explicit -> True] := 

-l/SPD[pl] SDEcl, c2]; 
P[c[cl_, pl_] , cb[c2_, p2_] , explicit -> True] := 

-SPD [pi] ~(\ [Delta] c - 1) SD[cl, c2] ; 
S[A[cl_, \[Mu]_, pl_] , A[c2_, \[Nu]_, p2_] , explicit -> True] : = 

SD[cl, c2] (MTD[\[Mu], \ [Nu] ] - 
FVD[pl, \[Mu]] FVD[pl, \[Nu]]/SPD[pl])/SPD[pl]; 
P[A[cl_, \[Mu]_, pl_] , A[c2_, \[Nu]_, p2_] , explicit -> True] : = 

SD[cl, c2] (MTD[\[Mu], \ [Nu] ] - 
FVD[pl, \[Mu]] FVD[pl, \[Nu]]/SPD[pl]) SPD [pi] ~(\ [Delta] A - 1); 
S[A[cl_, \[Mu]_, pl_] , cb[c2_, p2_] , c[c3_, p3_] , explicit -> True] : = 

1 g SUNF[cl, c2, c3] FVD[p2, \ [Mu] ] ; 
V[A[cl_, \[Mu]_, pl_] , cb[c2_, p2_] , c[c3_, p3_] , explicit -> True] := 

S[A[clA[Mu],pl], cb[c2, p2], c[c3,p3]]; 
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Here SPD, SD[a,b] and FVD represent scalar products, the color tensor 5°"^ and momenta. 
The choice of convention is up to the user and depends on the further processing of the 
expressions. The propagators were dressed with a power law and the dressed ghost-gluon 
vertex was taken bare. These functions have to be defined as above so DoDSE can use 
them. Their arguments are fields and the arguments of those are their indices and their 
momentum argument. The task of getAlg is to replace all the propagators and vertices 
in the result of a doDSE calculation by expressions with the correct indices and the correct 
flow of the momenta through the integral. After this it sets the option explicit to True 
and the algebraic expressions are plugged in. 

Let me illustrate this with the diagrams of the two-point functions in the Landau gauge 
containing ghosts. For the ghost two-point function this is 

IGhDSELoop = getAlg [ccDSE[ [2] ] , 
{{A, adj , lor>, {c, adj», 
{A, {c, cb», 

{pl, -pl>, 

{LorentzContract, SUNContract}-] ; 

and for the gluon two-point function 

IGluonDSELoop = getAlg [AADSE [ [4] ] , 
{{A, adj, lor>, {c, adj», 
{A, {c, cb>}, 

{pl, -pl}, 

{LorentzContract, SUNContract}]; 

The first argument of getAlg is the expression for the integrals we got from doDSE. The 
second tells which indices the fields have and the third is a list of the fields. The fourth list 
contains the external momenta of the integral. The final argument is a list of functions 
applied to the output. LorentzContract and SUNContract are user-defined functions for 
the contraction of the indices. 

With the resulting expressions one can calculate the IRE k of Landau gauge |29t 165] : 
The gluon DSE is projected transversely and the integrals are performed with the help of 



eq. (4.3). The two equations are then combined to solve for the IRE k. 



136 



E.2 DSEs of the maximally Abelian gauge 



E.2 DSEs of the maximally Abelian gauge 

The derivation of the two-point DSEs in the MAG is presented in this section and the 
way of discarding unphysical vertices is explained. 

E.2.1 Derivation of the DSEs in the maximally Abelian gauge 



First we define the action of the MAG based on eq. (5.31): 



IAMAG={{A, A}, {B, B>, {cb, c>, {db, d>, 
{A, B, B>, {A, cb, c>, {A, A, B, B>, {A, A, cb, c}, 
{B, B, B, B>, {B, B, cb, c>, {cb, cb, c, c>,{B, B, B>, 
{B, cb, c>, {A, B, cb, c>, {A, B, B, B}}; 

The fields A, B, c and d are the diagonal gluon, the off-diagonal gluon, the off-diagonal 
ghost and the diagonal ghost, respectively, cb and db are the corresponding anti-fields. 
DoDSE does not know anything about diagonal or off-diagonal indices and allows all 
combinations of fields. However, some structure functions vanish and only certain vertices 



exist, see Section 5.2.1 Hence we have to devise some test function that discards non- 
existent vertices. The function vertexTest only allows three-point functions with no or 
one diagonal field and four-point functions with at most two: 

Clear@vertexTest ; 

vertexTest [a_V] : = 
Not@(Length(§a == 3 kk MatchQ [LengthOCases [a. A, \ [Infinity]] , 2 I 3]) && 
Not@(LengthOa == 4 && MatchQ [LengthOCases [a. A, \ [Infinity]] , 3 I 4]) 

The derivation of the two-point DSEs can now be done. It is important to provide the 
function vertexTest as argument to discard unphysical vertices. Note that the indices 
are only given together with the fields to respect the convention of the MAG described in 
Section [5.2.1 Otherwise DoDSE would automatically take il and i2 as indices: 



AA = doDSE[IAMAG, {{A, i}, {A, j», vertexTest]; 
BB = Delete [doDSE[IAMAG, {{B, a>, {B, b}>, vertexTest], 23]; 
cc = Delete [doDSE[IAMAG, {{c, a>, {cb, b», vertexTest], 15]; 
dd = doDSE[IAMAG, {{d, i}, {db, j», vertexTest]; 
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With hindsight those diagrams are discarded manually in the off-diagonal equations which 



vanish due to their color structure, see Section [5. 2. 3 

For plotting the DSEs we can define some graphics primitives. The black and white def- 
inition as used in this thesis and a color definition which is better suited for presentations 
are 

fieldPlotRulesBW = {{A, Thickness [0. 01] >, {B, Thickness [0.01] , 
Dashing [0. 05] >, {c, Dotted}, {d, Dotted, Thickness [0.03] »; 

fieldPlotRulesColor = {{A, Red, Thickness [0.01] }, {B, Purple, 
Thickness [0.01]}, {c, DarkerSDarkerOGreen, Thickness [0.01] }, 
{d, LightGreen, Thickness [0.01]}}; 

When plotting the DSEs we use the argument f actorStyle to adjust the size and style 
of the numeric coefficients and set the number of diagrams per row to four: 

DSEPlot[#, IAMAGSU2, fieldPlotRulesBW, 4, 
f actorStyle :> {FontSize :> 20, FontWeight :> Bold}]& /0{AA, BB, cc, dd} 



This yields the two-point DSEs of the MAG as shown in figs. 5.1, 5.2 and 5.3 



E.2. 2 Obtaining and calculating the integrals 

The integrals of the two-point DSEs are obtained with the function get Alg. The Feynman 
rules of bare vertices are defined according to Appendix [C| Full propagators are dressed 
with a power law. Required arguments are a list with the fields and their indices, the list 
of fields and a list with the external momenta of the expression. For example, the sunset 
of the diagonal gluon two-point DSE with a bare A ABB vertex is obtained by 

getAlg[AA[[6]], 
{{A, adj , lor}, {B, adj , lor}, {c, adj}}, 
{A, B, {c, cb}}, 

-Cp, -p}, 

{LorentzContract , SUNMagContract}] ; 

The integral can be calculated with the solution for the sunset diagram given in Appendix 



D eq. (D.4), after a projection in Lorentz space to get a scalar expression. 
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E.3 DSEs of the Gribov-Zwanziger action 



E.3 DSEs of the Gribov-Zwanziger action 

The derivation of the two-point DSEs of the Gribov-Zwanziger action is presented in 
this Section. As in this action fields mix at the two-point level one needs the option 
specif icFieldDefinitions. Furthermore, this mixing requires an additional step to 
identify equal diagrams at the end. 

E.3.1 Deriving the DSEs of the Gribov-Zwanziger action 



The Gribov-Zwanziger action, see eq. (6.36), is defined for DoDSE as follows: 



lAFP = {{A, A}, {A, A, A}, {A, A, A, A}}; 

lAGZ = {{n, nb>, {A, nb, n>, {A, W, W>, {A, W}, {W, W}}; 

lA = JoinElAGZ, lAFP] ; 

The action contains the gluon field A, the pair of Grassmann fields n and nb and the 
auxiliary field W. They correspond to the fields A, i], ff and V of Chapter [6| respectivelyjj 
From the action we derive the individual propagator equations: 

AADSE = doDSEElA, {A, A>, {{A, A>, {n, nb>, {W, W>, {A, W>}, 

specif icFieldDefinitions -> {A, W, {n, nb}}] ; 
nnbDSE = doDSE[IA, {n, nb}, {{A, A}, {n, nb}, {W, W}, {A, W}}, 

specif icFieldDefinitions -> {A, W, {n, nb}}] ; 
WWDSE = doDSE[IA, {W, W}, {{A, A}, {n, nb}, {W, W}, {A, W}}, 

specif icFieldDefinitions -> {A, W, {n, nb}}] ; 
AWDSE = doDSEElA, {A, W}, {{A, A}, {n, nb}, {W, W}, {A, W}}, 

specif icFieldDefinitions -> {A, W, {n, nb}}] ; 
WADSE = doDSEElA, {W, A}, {{A, A}, {n, nb}, {W, W}, {A, W}}, 

specif icFieldDefinitions -> {A, W, {n, nb}}] ; 

Note that the option specif icFieldDefinitions is required to define W as a bosonic 
field. Otherwise DoDSE would assume automatically from the definition of the action 
that it is a fermion field. 

The number of terms in the DSEs can be determined with countTerms: 



""^The symbol W was chosen, as V represents in DoDSE a dressed vertex. 
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E.4 Tables of DoDSE functions 



countTerms /O {AADSE, nnbDSE, WWDSE, AWDSE, WADSE} 

{51, 3, 5, 50, 5} 

However, this is not the final number. Because of the mixed propagators not all diagrams 
were yet identified properly and this has to be done in a separate step: 

{AADSEId, nnbDSEId, WWDSEId, AWDSEId, WADSEId> = 

identifyGraphs [#, compareFunction :> compareGraphs2] & /@ 
{AADSE, nnbDSE, WWDSE, AWDSE, WADSE>; 

Now the number of terms decreases considerably for DSEs, where the first derivative was 
done with respect to the gluon field: 

countTerms /@ {AADSEId, nnbDSEId, WWDSEId, AWDSEId, WADSEId} 

{37, 3, 5, 36, 5} 

Finally, we can plot the diagrams 

pAADSE = DSEPlot [AADSEId, lA, 4, factorStyle -> {FontSize -> 20>] 
pnnbDSE = DSEPlot [nnbDSEId, lA, 4, factorStyle -> {FontSize -> 20}] 
pWWDSE = DSEPlot [WWDSEId, lA, 4, factorStyle -> {FontSize -> 20}] 
pAWDSE = DSEPlot [AWDSEId, lA, 4, factorStyle -> {FontSize -> 20}] 
pWADSE = DSEPlot [WADSEId, lA, 4, factorStyle -> {FontSize -> 20}] 



to obtain figs. 6J^, 6^, 6^ and 6A As we do not give any graphics primitives for drawing 
the propagators, DoDSE tags them by their field content. 

Finally, the full expressions are obtained with the function getAlg. After suitable 



projections in color and Lorentz space, the integrals can be performed using eq. (4.3). 



E.4 Tables of DoDSE functions 

In the following I provide lists with all public functions of DoDSE. 
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E.4 Tables of DoDSE functions 



Main functions 



Command 



Description 



doDSEEilist, clist] 

doDSEEilist, clist [, 
props, vertexTest, opts]] 



shortExpression[expr, opts] 
sE[expr, opts] 



DSEPlotEexpr, 

ilist [,f Rules, len, opts]] 



getAlgEexp, ilist, flist, 
mlist, funclist 



Derives tlie DSE for the correlation function 
clist for a tlieory with interactions ilist. 
vertexTest is a function for determining if 
a vertex respects the symmetries of the La- 
grangian. props is a hst of allowed propaga- 
tors given in the form {{fieldla, fieldlb}-, 
{field2a, field2b>, ...>. doDSE accepts 
the options specif icFieldDefinitions and 
sourcesZero (prevents the replacement of 
super-field propagators and vertices when set to 
False). 

Rewrites a DoDSE expression into a 
shorter form using $bareVertexSymbol, 
$vertexSymbol and $propagatorSymbol 
for representation. Options of the internal 
Mathematica function Style can be given. 
Plots the full DSE in graphical form, expr is 
an expression containing op functions, ilist 
the list of interactions and fRules a list of op- 
tions for plotting individual fields, len deter- 
mines how many graphs are shown in one line. 
If fRules is not given, the lines are named 
according to the fields. Possible options are: 
output ->Li St, to get the result in list form, 
and indexStyle and f actorStyle to change the 
style of the indices and the prefactors (e.g. font 
size or color). 

Derives the algebraic expression from exp. 
ilist contains the lists of fields and all their 
indices, flist the list of fields, mlist the list of 
external momenta and funclist a list of further 
functions to be applied on the result. 
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E.4 Tables of DoDSE functions 



Functions for the individual computation steps 



Command 



Description 



generateAction[ilist [,f list]] Generates the action in internal representation 

from the interactions of the theory given in 
ilist. For mixed propagators f list specifies 
explicitly the type of fields in the form {bosonl , 
boson2, ..., {fermionl, antif ermionl}-, 
{fermion2, antif ermion2}-, . . .}. 
Differentiate expr with respect to the fields in 
dlists. 

Replaces the fields in expr by the corresponding 
expressions after the first differentiation is done 
to change from full to IPI Green functions. 
Adds up equivalent graphs in expr. cfunc 
can be compareGraphs (standard) or 
compareGraphs2, the latter being necessary for 
mixed propagators but taking longer. 
Sets the external fields in flist to zero, i.e. 
only physical propagators and vertices are left. 
vertexTest is a function for determining if 
a vertex respects the symmetries of the La- 
grangian. props is a list of allowed propaga- 
tors given in the form {{fieldla, fieldlb}-, 
{field2a, field2b>, . . .}. 
Orders derivatives with respect to Grassmann 
fields such that the anti-fields are left of the fields 
thereby possibly giving a minus sign, expr is an 
op-function or a sum of those. Bare vertices are 
not affected by the ordering. 



deriv [expr, dlists] 
replaceFields [expr] 



identif yCraphs [expr [ , 
compareGraphs -> cfunc] ] 



setSourcesZero[expr, flist 
[, props, vertexTest]] 



orderFermions [expr] 
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E.4 Tables of DoDSE functions 



Functions for checks and tools 


Command 


Description 


countTerms [expr] 


Counts tlie number of terms appearing in the 




expression. 


fieldQEf] 


Determines if expression f is defined as a field. 


bosonQEf ] 


Determines if expression f is defined as a bosonic 




field. 


f ermionQ [f ] 


Determines if expression f is defined as a 




fermionic field. 


antiFermionQ [f ] 


Determines if expression f is defined as an anti- 




field to a fermionic field. 


checkFields [expr] 


Checks if all fields in the expression are defined 


checklndices [expr] 


Checks if an index appears more often than 




twice. 


checkSyntax [expr] 


Checks if expr has the correct syntax, i.e. op 




functions only contain propagators, vertices and 




fields. 


checkAction [expr] 


Checks if all indices appear exactly twice, the 




syntax is ok and all fields are defined. 


checkAll [expr] 


Performs a series of checks on expr 




(checklndices, checkSyntaix, checkFields). 


def ineFields [f list] 


Defines the fields of the action that are given in 




f list as single entries for bosons and grouped 




by braces for fermions. 


$vertexSymbol 


Symbol representing a vertex in 




shortExpression. Standard value: F. 


$bareVertexSymbol 


Symbol representing a bare vertex in 




shortExpression. Standard value: S. 


$PropagatorSymbol 


Symbol representing a propagator in 




shortExpression. Standard value: A. 
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